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ABSTRACT 

This is the second paper in a series where we study the influence of transport processes on the chemical 
evolution of protoplanetary disks. Our analysis is based on a 1+1D flared a-model of a ~ 5 Myr DM Tau-like 
system, coupled to a large gas-grain chemical network. To account for production of complex molecules, the 
chemical network is supplied with an extended set of surface reactions and photo-processes in ice mantles. Our 
chemo-dynamical disk model covers a wide range of radii, 10-800 AU (from a Jovian planet-forming zone 
to the outer disk edge). Turbulent transport of gases and ices is implicitly modeled in full 2D along with the 
time-dependent chemistry, using the mixing-length approximation. Two regimes are considered, with high and 
low efficiency of turbulent mixing. The results of the chemical model with suppressed turbulent diffusion are 
close to those from the laminar model, but not completely. A simple analysis for the laminar chemical model to 
highlight potential sensitivity of a molecule to transport processes is performed. It is shown that the higher the 
ratio of the characteristic chemical timescale to the turbulent transport timescale for a given molecule, the higher 
the probability that its column density will be affected by diffusion. We find that turbulent transport enhances 
abundances and column densities of many gas-phase species and ices, particularly, complex ones. For such 
species a chemical steady-state is not reached due to long timescales associated with evaporation and surface 
photoprocessing and recombination (t > 10 5 years). When a grain with an icy mantle is transported from a cold 
disk midplane into a warm intermediate/inner region, heavy radicals become mobile on the surface, enriching 
the mantle with complex ices, which are eventually released into the gas phase. The influence of turbulent 
mixing on disk chemistry is more pronounced in the inner, planet-forming disk region where gradients of 
temperature and high-energy radiation intensities are steeper than in the outer region. In contrast, simple 
radicals and molecular ions, which chemical evolution is fast and proceeds solely in the gas phase, are not 
much affected by dynamics. All molecules are divided into three groups according to the sensitivity of their 
column densities to the turbulent diffusion. The molecules that are unresponsive to transport include, e.g., 
C 2 H, C + , CH 4 , CN, CO, HCN, HNC, H 2 CO, OH, as well as water and ammonia ice. Their column densities 
computed with the laminar and 2D-mixing model differ by a factor of < 2-5 ("steadfast" species). Molecules 
which vertical column densities in the laminar and dynamical models differ by up to 2 order of magnitude 
include, e.g., C 2 H 2 , some carbon chains, CS, H 2 CS, H 2 0, HCO + , HCOOH, HNCO, N 2 H + , NH 3 , CO ice, 
H 2 CO ice, CH3OH ice, and electrons ("sensitive" species). Molecules which column densities are altered by 
diffusion by more than 2 orders of magnitude include, e.g., C 2 S, C3S, CeH6, C0 2 , 2 , SiO, SO, S0 2 , long 
carbon chain ices, CH3CHO ice, HCOOH ice, 2 ice, and OCN ice ("hypersensitive" species). The chemical 
evolution of assorted molecules in the laminar and turbulent models is thouroughly analyzed and compared 
with previous studies. We find that column densities of observed gas-phase molecules in the DM Tau disk 
are well reproduced by both the laminar and the chemo-dynamical disk models. The observed abundances of 
many reduced and oxidized cometary ices are also successfully reproduced by the both models. We indicate 
several observable or potentially detectable tracers of transport processes in protoplanetary disks and the Solar 
nebula, e.g., elevated concentrations of heavy hydrocarbon ices, complex organics, C0 2 , 2 , SO, S0 2 , C 2 S, 
C3S compared to CO and the water ice. A combination of UV photodesorption, grain growth, and turbulent 
mixing leads to non-negligible amount of molecular gases in the cold disk midplane. 

Subject headings: accretion, accretion disks — astrochemistry — molecular processes — protoplanetary disks 
— stars: DM Tau — turbulence 



1. INTRODUCTION 

One of the most exciting questions in astrophysics is the 
genesis of prebiotic molecules served as life-building blocks 
in the Solar system and preserved in meteorites and comets. 
Nowadays protoplanetary disks are believed to be birth places 
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of planetary systems, so we may expect prebiotic molecules 
to be the eventual outcome of the disk evolution. How- 
ever, our understanding of the chemical composition and 
evolution of protoplanetary disks is far from being com- 
plete. Apart from CO and its isotopologues, and occa- 
sionally HCO + , DCO + , CN, HCN, DCN, CCH, H 2 CO, and 
CS, the molecular co ntent of protoplanetary disks remains 
largely unknown (e . g..lDutrey et al.ll997tlKastner et aljl997t 
lAikawa et alJl2003llThi et alJ2004llPietu et alJ2007t lOietal] 
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l2008tlHenning et alfe OlO). Molecular line data are limited in 
their sensitivity and resolution. This means that spatial distri- 
bution of m olecular abundances in disks is still poorly deter- 
mined (e.g., Pietu et al. 2005; Dutrey et al. 2007b; Panic et al. 
2009), hampering a detailed comparison with existing chem- 
ical models, which is often based on global data (e.g. inte- 
grated line profiles). 

Multi-molecule, multi-transition interferometric observa- 
tions, coupled to line radiative transfer and chemical mod- 
eling, allowed to constrain disk sizes, kinematics, distribution 
of temperature, surfac e density, and molecular c olumn den- 
sities (s ee reviews by Bergin et al.l (120071) and iDutrev et alj 
(2007a)). The measured line intensities a re indicative of ver- 
tical temper ature gradients in disks (e.g.. iDartois et al.ll2.003b 
lOi et al.ll2006l) . though several disks with large inner cavities 
do not sho w evidence for such a gradient (e.g., G M Aur and 
LkCa 15; IDutrev et"aill200l iHughes et alj |2009). A signif- 
icant reservoir of very cold CO, HCO + , CN and HCN gases 
has been found in the disk of DM Tau at temperatures < 6- 
17 K, which cannot be explained by conventional chemical 
models without invoking a non-thermal desorption or trans- 
port mechanism (e .g., iSemenov et al.l I2006t lAikawal 120071: 
iHersant et al.l 12009). Non- thermal broadening of emission 
lines of ~ 0.1 kms" 1 has been reported (e. g., [Bergin et al. 
l2007tlDutrev et al.l2007atlHughes et alJ20ll . which is likely 
due to sub sonic turbulence driven by the magnetorotational 
instability dBalbus & Hawlevll 19911) . 

Recently, with space-borne (Spitzer) and ground-based 
(Keck, VLT, Subaru) infrared telescopes, molecules have 
been detected in very inner zones of planet-forming sys- 
tems, at r < 1-10 AU. Rotational-vibrational emission 
lines from CO, C0 2 , C 2 H 2 , HCN, OH, H 2 imply a rich 
chemistry driven by endothermic reactio ns or reactions with 
activation barriers an d photoproc esses (lLahuis et al.l 120061; 
Carr & Naiital 120081: fSalyk et all 120081: IPontoppidan et ail 
2008al; iPascucci et al.l 120091; Ivan der Plas et all 120091: 
Salvk et al.l 1201 11) . Through ISO and Spitzer infrared spec- 
troscopy abundant ices in cold disk regions consisting of 
water ice and substantial amounts (~ 1 - 30%) of volatile 
materials like CO, CO ?, NH 3 , CH 4 , H 2 CO, and HCOOH have 
been detected (e.g.. IPon toppidan et alj 12005: Terada et al. 
120071; IZasowski et al JI2009I ). 

The conditions of planets formation in the early Solar sys- 
tem have been revealed by a detailed analysis of chemi- 
cal and mineralogical composi tion of meteo ritic samples and 
cometary dust particles (e.g., Bradley 2005). The recent Star- 
dust and Genesis space missions have returned first samples of 
pristine materials, likely of cometary origin, showing a com- 
plex structure of high-temperature cryst alline silicates embed- 
ded in low-temp e rature condensates (|Brownlee et al.l 120041 
iFlvnn et alj|2006t iBrownlee et all [2008). Comets have been 
assembled around or beyond Neptune and expelled gravita- 
tionally outward, however the presence of Mg-rich crystalline 
silicates in cometary dust indicates annealing of amorphous 
preso l ar gra ins at temperatures above 800 K dWooden et al.l 
[19991 120071) . The presence of crystalline silicates in outer 
regions of protoplanetary disks has also been revealed (e.g., 
Ivan Boekel et al.ll2004Uuhasz et al.ll20loh . 

Recently, an anti-correlation between the age of a disk and 
the X-ray hardness/luminosity of a central star, and the ob- 
served crystallinity fraction has been inferred , making the 
overa ll picture even more complicated (e.g., Glaus er et al.l 
2009). 



An isotopic analysis of refractory condensates in un- 
altered chondritic meteorites shows strong evidence that 
the inner part of the Solar Nebula has been almost com- 
pletel y mixed during the first several Myr of evolution 
(e.g., IBossI 12004 iCieslal 120091) . This mixing (either ad- 
vective or turbulent), along with high-energy irradiation, 
could have also been important for fractionation ratios in 
both gas-phas e and solid compounds in the nebu l a (e.g. , 
[ Clayton! 119931; IClayton & Maved ai 1 19961; ILee et all Il998bt 
Hersa ntet al.H200U iLvons & Y oung 2005|). The rich variety 
of organic compounds in meteorites, including amino acids, 
suggest that these complex species have formed just prior or 
during the formation of planets in heavily irradiated, warm 
regions of the Solar Nebula (e.g., Ehrenfreund & Charnley 
2000; Bus emann et alj|2006l) . Combustion and pyrolysis of 
hydrocarbons at high temperatures, coupled with outward 
transport, has been inferred to explain the omni-presence 
of kerogene-like (mainly aromatic) carbonaceous mat erial in 
meteoritic and cometary samples (Morg arTet al.lll99Tb . 

These intriguing findings are partly understood 
in modern astrochemical model s of protop la netary 
disks dwTllacy & Langeri l2000t lAikawa et all 120021 
i Markwick et al.N2002t Ivan Zadelhoff et al.H2003l: lllgner et al l 
20041: iKamp & Dullemondl 12004 ISemenov et alj 1 2005 1: 
Aikawa & Nomura! 120061: iTschamuter & Gaill 120071 
Agun dez et alj 120081: IWoods & Willacvl 120091; IVisser et all 
2009b; Wa lsh et al] |2010). The major result of the chemical 
modeling is that disks have a layered chemical structure due 
to heavy freeze-out of gas-phase molecules in the cold mid- 
plane and their photodissociation in the atmosphere. Vertical 
column densities of CO, HCO + , N 2 H + , CN, HCN, HNC, 
CS, etc. are reproduced with modern chemical models (e.g., 
Aikawa et alJl2002tlSemenov et al.ll200~5l; I Willacy et alj|2006t 
Dutrey et aLl I2007bt iSchrever et all 120081 iHenning et alj 
2010)'. 

While most of the chemical studies are still based on lam- 
inar disk models, evidences for mixing, mentioned above, 
call for a more sophisticated treatment. A few such dy- 
namical studies have been presented. Models of the early 
Solar nebula with ra dial transport by advective flow s have 
been developed (e . g..lMorfill & VoM19"84tlGailll2001ll2002l; 
Wehrstedt & Gail 120021: IBossI 12004 iKeller & Gaill 120041) . 
Ilgner et al.1 (120041) for the first time simultaneously modeled 
the influence of turbulent diffusion in the vertical direction 
and advection flows in the radial direction on the chemical 
composition of the inner disk region. They found that dy- 
namical processes signific antly affect the chemi cal evolution 
of sulfur-bearing species. Willac y et al.l (120061) have shown 
that ID vertical mixing modifies chemical composition of 
the outer disk re gion and that the mixin g results bette r agree 
to obs ervations. ISemenov et al.l (|2006, Paper I) and|Aikawa 
(2007) have found that turbulent transport allows explain- 
ing the presence of a large amount of cold (< 15 K ) CO 
gas in the disk of DM Tau. Tscha rnuter & Gaill (120071) have 
used a 2D disk chemo-hydrodynamical model and showed 
that in the disk midplane matter moves outward, carrying 
out the angular momentum, while the accretion flows to- 
ward the star are located at elevated altitudes. Consequently, 
gas-phase species produced by warm chemistry in the in- 
ner nebula can be st eadily transported into the cold outer re- 
gion and freeze out. Tscharnu ter & GaiH (12007b have claimed 
that global radial advective flows dominate over diffusive 
mixing for the disk chemical evoluti on. A radial advec- 
tion model has also been utilized by iNomura et alj (120091) . 
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who have demonstrated that inward radial transport enhances 
abundances of organic molec ules (produced mainly o n dust 
surfaces in cold outer regions). Woods & Willacy (2009) have 
elaborated a disk chemical model with improved heating and 
cooling balance and accurate modeling of the UV radiation 
field and found that the I2 C/ I3 C fractionation in the Solar sys- 
tem comets can be explained by the re processing of preso- 
lar materials in a warm nebular region. lHersantetal1(B009l) 
have studied various mechanisms to retain gas-phase CO in 
very cold disk regions. They concluded that efficient pho- 
todesorption in moderately obscured disk regions (Ay < 5 m ) 
greatly enhances gas-phase CO concentration s, while the role 
of vert ical mixing is less important. Finally. iHeinzeller et al.1 
d201 ll) have investigated the disk chemical evolution with ra- 
dial advection, vertical mixing, and vertical wind transport 
processes. They have found that the disk wind has a negli- 
gible effect on disk chemistry, whereas the radial accretion 
alters the molecular abundances in the cold midplane, and the 
vertical turbulent mixing affects the chemistry in the warm 
molecular layer. The abundances of NH3, CH3OH, C2H2 and 
sulfur-containing species are the most enhanced by the trans- 
port. 

In this paper we continue our detailed study of chemo- 
dynamical interactions in protoplanetary disk started in Pa- 
per I with an intent to find out if chemistry can be used as a 
diagnostic of dynamical processes in a protoplanetary accre- 
tion disk. For the first time we utilize a large-scale disk physi- 
cal model along with an extended gas-grain chemical network 
coupled to 2D turbulent transport. A wide range of tempera- 
tures, densities, and X-ray/ UV radiation intensities encoun- 
tered in this dynamical model allows us to follow formation 
and destruction of various molecules, possibly detectable with 
ALMA and Herschel. The primary aim of the present study 
is to characterize the importance of turbulent diffusion for the 
chemical evolution in various disk domains and for various 
chemical families. We argue that even though the overall effi- 
ciency of the diffusive transport in the outer disk is dominated 
by vertical mixing, one has to consider vertical and radial 
mixing simultaneously in the planet-forming region. Many 
complex ices and their gas-phase counterparts are enhanced 
by turbulent diffusion, in particular sulfur-bearing and other 
heavy (complex) species. 

The organization of the paper is the following. In Sec- 
tion |2] we describe the adop ted d isk physical model and the 
chemical network. In Sect. 13.11 basic chemical and dynam- 
ical timescales in protoplanetary disks are outlined and dis- 
cussed. A general scheme to estimate possible sensitivity 
of a given molecule to transport processes is presented in 
Sect. 13. 21 Influence of the 2D-turbulent diffusion on the chem- 
ical evolution of dominant ions as well as C-, O-, N-, S- 
bearing species and com plex organic molecules is studied in 
detail in Sects. I4.1H4.6I Detailed comparison with the pre- 
vious studies and future impro vem ents of chemo-dynamical 
models are discussed in Sect. 15.11 We discuss importance 
of turbulent diffus ion f or the presence of cold gases in disk 
midplanes in Sect. 15.21 We verify feasibility of our chemical 
and physical disk models by comparing the calculated and ob- 
served column densities in the DM Tau dis k an d abundances 
of cometary ices in the Solar system (Sect. 15.3b . Finally, de- 
tected or potentially detectable molecular tracers of trans port 
processes in protoplanetary disks are summarized in Sect. 15.41 
Summary and conclusions follow. 

2. MODEL 



2.1. Disk structure 

We focus on the DM Tau system because it is one of the 
most observationally studied and molecularly-rich protoplan- 
etary disk around a T Tauri star. The adopted flaring disk 
structure is based on a 1+1D steady-state a-model similar to 
iD'Alessio etail dl999t) model. The DM Tau is a single M0.5 
dwarf ( T s ff = 3720 K), with a mass of 0.65M f?), and a radius of 
\.2Rq (lMazzitellilll989l:ISimon et"aTll2000l) . The non-thermal 
FUV rad iation fi e ld from DM Tau is represented by the scaled 
ISRF of Draine dl978l) . with t he un-attenuated in tensity at 
100 AU of x*(100) = 410 (e.g. JBergin etafl 120031) . For the 
X-ray luminosity of the star we adopt a value of 10 30 erg s" 1 , 
which is constrained by recent measurements with Chandra 
and XMM in the range of 0.3-10 keV (M. Guedel, priv. 
comm.). 

The disk has an inner radius ro = 0.03 AU (dust subli- 
mation front, T w 1 500 K), an outer radius r\ = 800 AU, 
an accretion rate M = 410~ 9 Mq yr , a viscosity parameter 
a = 0.01, and a mas s of M = O.O66M (iDutrey et al.ll2007bt 
lHenning etapl2010i) . The DM Tau disk age is about 5-7 Myr 
dSimon et akteOOoHh so we adopted 5 Myr as a limiting time in 
our che mical simulations. According to the Spitzer IRS obser- 
vations dCalvet et al.l l2005). the inner DM Tau disk is cleared 
of small dust (< 3 AU) and is in a pre-transitional phase 
already. Therefore, in the chemical simulations a disk re- 
gion beyond 10 AU is considered, where dust grain evolution 
seem s to be slow and grain growth is moderate (Birnstiel et al. 
2010). In chemical modeling, the dust grains are assumed to 
be uniform 0.1 /jm amorphous olivine particles (with density 
of pd = 3 gem" 3 ). This is the size representing a mean radius 
in the dust size distribution. Equal gas and dust temperatures 
are assumed. Gas becomes hotter than dust only in an up- 
per, tenuous and heavily irradiated disk layer, which often has 
negligible contribution to molecular column densities. 

The turbulence in disks is likely driven by the magnetorota- 
tional instability (MRI), which is operative even in a weakly 
ionized medium, and is essentially a 3-D phenomenon (e.g., 
Balbus & Hawley H"l991l) . This turbulence causes anomalous 
viscosity that enables efficient redistribution of the angular 
momentum and accretion of disk matter onto the star. How- 
ever, inner disk midplane can be almost neutral and decou- 
pled from magnetic fields, formin g a region with redu ced, in - 
herited turb ulence (" dead zone"). iGammiel (see, e.g., 1996); 
Sano et all (se e, e.g.. I2000T): ISemenov et al.1 (see, e.g., 2004); 
Wiinsch et alJ (see, e.g.. 120061). With modern computational 



facilities coupled chemo-MHD 3-D models are only manage- 
able for extremely limited che mical networks and a restricted 
disk domain (see, for example, Turner et al. 200 7]). We hav e 
followed the parametrization of Shakura & Suny aevl (119731) . 
where turbulent viscosity v is related to local disk proper- 
ties such as the characteristic (vertical) spatial scale H(r), the 
sound speed c s (r,z), and the dimensionless parameter a: 



v(r,z) = ac s (r,z)H(r). 



(D 



From observational cons t raints a is ~ 0.0 01 - 0.1 
dAndrews & Williams! 120071; iGuilloteau et al.1 1201 lb , similar 
to values obtained from MHD modeling of the MRI (e.g., 
iDzyurkevich et al.1 120 101) . However, the magnitude of MHD 
viscous stresses changes throughout the disk, and thus in this 
simplistic parametrization the a-parameter should also be 
variable. Unfortunately, without detailed MHD studies it is 
hard to characterize a, so we adopt the constant value of 0.01 . 
Since only a disk region beyond 10 AU is studied, a "dead 
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Figure 1. (Left to right) Distributions of the temperature, particle density (log 10 scale), diffusion coefficient in cm 2 s ' (log 10 scale), and pressure scale height in 
AU (logjQ scale) in the DM Tau disk model. In the first 3 panels the Y-axis is given in units of the pressure scale height. 



zone" where effective a is very low, < 10 , is avoided. 
Consequently, the diffusion coefficient is calculated as 

AurbO, z) = v(r, z)/Sc, (2) 

where Sc is the Schmidt num ber describing efficiency 
of turbulent diffusivity (se e e.g. [Shakura & Sunyaevl 119731: 
Schra pler & Hen ning 2004). In our simulations we assume 
that gas-phase species and dust grains are well mixed, and 
transported with the same diffusion coefficient. We treat dif- 
fusion of mantle materials similarly to gas-phase molecules, 
without relating it to individual grain dynamics. Hence 
two chemo-dynamical models are considered: (1) the high- 
efficiency mixing model with Sc = 1 and (2) the low-efficiency 
mixing model with Sc = 100. In the Sc = 1 model the diffusion 
coefficient in outer d isk regions is ~ 10 18 cm 2 s" 1 , similar to 
IWillacv et al.l (12006b . The second model represents a hypo- 
thetical case when mixing of molecules occurs much slower 
than turbulent eddy turn-over speeds yet faster than in the pure 
laminar gas. The temperature and density structure, diffusion 
coefficient D tur \,, and vertical pressure scale height H(r) of the 
DM Tau-like disk model are shown in Fig.Q] 

2.2. Chemical network 

The adopted gas-grain chemical model is described in our 
recent papers on ben chmarking of disk chemical models 
dSemenov et al.ll2010bl) a nd observations of CC H in DM Tau, 
LkCa 15, and MWC 480 dHenning et alj|2010h . A brief sum- 
mary is provided below. The chemical network is based on 
the osu.2007 ratefile with recent updates to reaction ratesQ. 
(Note that in Paper I we used UMIST 95 ratefile.) A new class 
of X-ray-driven reactions leading to production of O" 1 " 1 ", C ++ , 
N ++ , S ++ , Fe" 1 " 1 ", Si" 1 " 1 " is added. Their neutralization reactions 
by electrons and charge transfer re actions with molecules are 
adopted from Stauber et al. (20 051). The phot oionization cross 
sections are taken fro m IVerner et all ([1993 ), as described in 
iMalonev et alj d!996f) . Seco ndary electron impact ion ization 
cross sections are taken from Meijerink & Spaans (2005). 

To calculate UV ionization and dissociation rates, the mean 
FUV intensity at a given disk location is obtained by adding 
the stellar x*( r ) = 410(r,AU)/(100 AU) 2 and interstellar xo 
components that are scaled down by the visual extinction 
in the vertical direction and in the direction to the central 
star (ID plane-parallel approximation). Several tens of pho- 
toreaction rates are upda ted using the new calculations of 
Ivan Dishoeck et afl (l2006). which are publicly available^. The 
self-shielding of H2 from photodissociation is calculated by 

2 See:|http : / /www .physics ■ ohio- state ■ edu/$\sim$eric7 
3 |http :77wwwT^trwTTeidenuniv7nT7~ewin^ 



Eq. (37) from lDraine & Bertoldil (11996b . The shielding of CO 
by dust grains, H2, a nd its self-shield ing is calculated using 
precomputed table of lLee et ail (119961 Table 11). 

We model the attenuation of cosmic rays (CRP) by Eq. (3) 
from ISemenov et al.1 (120041) . using the standard CRP ion- 
ization rate £cr = 1.3 10~ 17 s" 1 . Ionization due to the de- 
cay of short-living radionuclides is taken into account, £rn = 
6.5 • 10" 19 s" 1 (iFinocchi et alJll997l) . The ste llar X-ray radia- 
tion is modeled using observational results of Glassgold et al. 
(2005) and the approximate expressions (7-9) from the 2D 
Monte Carlo simulations of Glass goTd et al.l dl997alfbT) . The 
typical X-ray photon energy is 3 keV, and the X-ray emitting 
source is located at 12 stellar radii above the midplane. The 
X-ray ionization rates exceed that of the CRPs in the disk re- 
gions above the midplane. 

The gas-grain interactions include sticking of neutral 
species and electrons to uniformly-sized dust grains with 
100% probability, release of frozen molecules by thermal, 
CRP-, and UV-induced desorption, dissociative recombina- 
tion and radiative neutralization of ions on charged grains, 
and grain re-charging. We do not allow H2 to stick to grains 
because th e binding energy of H2 to pure H2 mantle is low, 
~ 100 K dLedll972l) . and it freezes out in substantial quan- 
tities only at temperatures below ss 4 K. Chemisorption of 
surface molecules is not considered. We considered vari- 
ous UV photode sorption yields betw een 10~ 5 and 10~ 3 (e.g., 
iGreenbergll 1 9731: IQberg et"aDl2009bllah and found that in this 
range the exact yield value has a negligible impact on the 
modeling results. To allow synthesis of complex (organic) 
molecules, an extended list of surface reactions together with 
desorption energies and a list of p hotodissociation reaction s 
of surface species is adopted from iGarrod & Herb st (2006). 
Desorption energies for assorted molecules are listed in Ta- 
ble Q] We assume that each 0.1/.tm spherical olivine grain 
provides ps 2 10 6 surface sites, and that surface recombination 
proceeds solely through the Langmuir-Hinshelwood forma- 
tion mechanism. Upon a surface recombination, there is a 5% 
chance for the products to leave the grain. Following inter- 
pretations of experimental resu lts on the formati on of molec- 
ular hydrogen on dust grains (Katz et aD 1 19991) . we employ 
the standard rate equation approach to the surface chemistry 
without H and H2 tunneling either through the potential walls 
of the surfa ce sites or through re action barriers. As it has been 
shown bv lVasvunin et alj d2009l) . when surface recombination 
rates are slow, the stochastic effects are of no importance for 
the surface chemistry. 

Overall, the disk chemical network consists of 657 species 
reseSBadetoJml3 elements, and 7306 reactions. In contrast to the 
Paper I, we do not apply our 'Automatic Reduction Tech- 
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Table 1 

Desorption energies 



Table 2 

Initial abundances 



Species 



Energy, K 



C 

Cl 

&H 

C,H 2 
C 2 S 

c 3 

C 3 H 2 

C 6 H 6 

C 8 H 2 

CHoOH 

OfcCO 

CH 3 CHO 

CH^OH 

CH 4 

CN 

CO 

COt 

cs 

H 

H 2 

H 2 S 

Hi CO 

H,CS 

H,0 

HCN 

HCOOH 

HNC 

HNCO 

HNO 

N 

N 2 

NH 

NH, 
NH 3 
NO 
O 

O, 

OCN 

OH 

S 

SO 

so 2 



800 

1600 

2140 

2590 

2700 

2400 

3390 

7590 

7390 

5080 

2200 

2870 

5530 

1300 

1600 

1150 

2580 

1900 

624 

552 

2740 

2050 

2700 

5700 

2050 

5570 

2050 

2850 

2050 

800 

1000 

2380 

3960 

5530 

1600 

800 

1000 

2400 

2850 

1100 

2600 

3400 



nique" (ART) to reduce the size of this network (thanks to 
increased performance of our che mo-dynam i cal cod e). The 
"low metals" initial abundances of lLee etail (11998a!) are uti- 
lized (Table |2). The choice of initial abundances does not 
affect the resulting molecular abundances and column densi- 
ties due to relatively high densities in and long evolutionary 
timescales of protopla netary disks, which essentially reset a 
chemical "clock" (e.g.. Willa cy et al.ll2006l) . 

2.3. Modeling chemistry with transport 

The disk physical structure and the chemical model de- 
scribed above are used to solve chemical kinetics equations 
together with turbulent transport terms. Since a while such 
kind of models have been employed in atmospheric chemistry, 
planetary atmosp here chemistry, and chemistry of molec- 
ular clouds (e.g., [G arcia & Solomon 1983; [Xie et al.l 11995b 
Yate & Millar 2003). We utilize the mixing approach of 
Xie et alj d 1995b . which is based on the Fickian's diffusion 
laws: 



dni 
dt 



(r, z) = Fi(r, z)-Li(r, z)- V ■ Aufcfo z)n H (r, z) V 



nj(r,z) 
nu(f,z) 



Species 


Relative abundance 


H 2 


0.499 




2 00(— 3) 


He 


9.75(-2) 


C 


7.86(-5) 


N 


2.47(-5) 


O 


1.80(-4) 


S 


9.14(-8) 


Si 


9.74(-9) 


Na 


2.25(-9) 


Mg 


1.09(-8) 


Fe 


2.74(-9) 


P 


2.16(-10) 


CI 


1.00(-9) 



Table 3 

Characteristic Timescales: Inner Disk (10 AU) 



Processes 


Midplane 


Warm layer* 


Atmosphere* 




[yr] 


[yr] 


[yr] 


Mixing 


3.4 (3) 


3.4 (3) 


1.3(3) 


Gas-phase 


1.4 (-5) 


1.3 (-4) 


1.0 (-2) 


UV 


>1.0 (7) 


3.3 (4) 


5.9 (2) 


Accretion 


1.2 (-2) 


1.1 (-1) 


5.4 (0) 


Desorption 


5.8 (-7) 


6.0 (-7) 


<1.0 (-7) 


Surface 


>1.0(7) 


>1.0(7) 


<1.0(-7) 


* The warm 


layer and atmosphere are located at the z/H t 


= 0.8 and 1.75, 


respectively. 












Table 4 






Characteristic Timescales: Outer Disk (250 AU) 


Processes 


Midplane 


Warm layer* 


Atmosphere* 




[yr] 


[yr] 


[yr] 


Mixing 


1.0 (6) 


2.5 (5) 


1.4 (5) 


Gas-phase 


2.0 (-2) 


1.8 (-1) 


2.9 (0) 


UV 


>1.0 (7) 


1.2(6) 


3.1(1) 


Accretion 


2.7 (1) 


1.8 (2) 


2.3 (3) 


Desorption 


1.0 (6) 


4.3 (0) 


<1.0 (-7) 


Surface 


>1.0(7) 


>1.0(7) 


1.4 (5) 



The warm layer and atmosphere are located at the z/H t 
respectively. 



0.8 and 1.75, 



Here n, is concentrations of the i-th species (cm -3 ), Fj and 
Li are formation and destruction (loss) terms. In the follow- 
ing we will use relative abundances for considered species 
X/ = n,/«H where «h is the total hydrogen nucleus number 
density. In this formalism the turbulence mixing rate for a 
certain species depends on its chemical gradient. 

The formation and destruction of molecules are governed 
by the chemical kinetics: 



dn; 
dt 



1. 111 



(4) 



dn 
dt 



(3) 



where n\ is the surface concentration of the 2-th species 
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Figure 2. (Top to bottom) Distributions of timescales to reach a chemical steady-state during the 5 Myr of evolution for C + , CO, CO?, CH3OH ice, NH3, and 
SO2. For relative abundances below 1CT 25 this timescale is assumed to be 1 year. 



(cm ), ki m and ki are the gas-phase reaction rates (in units 
of s" 1 for the first-order kinetics and cm 3 s _1 for the second- 
order kinetics), kf cc and kf es denote the accretion and des- 
orption rates (s _I ), and k\ m and k\ are surface reaction rates 
(cm 3 s" 1 ). 

The stiff equations of chemical kinetics are integrated si- 
multaneously with the diffusion terms in the Eulerian descrip- 
tion, using a fully implicit 2D integration scheme. As bound- 
ary conditions for mixing, we assume that there is no inward 
and outward diffusion across boundaries of the disk domain, 
and that there is no flux through the midplane. We do not 
employ approximate, operator-splitting integration schemes, 
in which transport and chemical processes are treated sep- 
arately, and instead integrate PD E system directly. O ur 
"ALCHEMIC" FORTRAN77 code dSemenov et aT]|2010bl) is 
based on the Double-precision Variable-coefficient Ordinary 
Differential equation solver with the Preconditioned Krylov 
(DVODPK) method GMRES for the solution of linear sys- 
temfl The approximate Jacobi matrix is generated automat- 
ically from the supplied chemical network (without transport 
terms) and serves as a lefthand preconditioner. For astrochem- 
ical models dominated by hydrogen reactions the Jacobi ma- 
trix is sparse, with < 1% of non-zero elements. The corre- 
sponding linearized system of algebraic equations is solved 
using a high-performance sparse unsymmetric MA48 solver 
from the Harwell Mathematical Software LibraryB All the 
equations are solved on a non-uniform staggered grid con- 
sisting of 41 radial points (from 10 to 800 AU) and 91 verti- 



http : / /www . netlib ■ org/ ode/vodpk . f| 



httpTT7ww^7hsTTTTT^cTuk7[ 



cal points. This resolution is found to be optimal for mixing 
problems in a protoplanetary accretion disk, keeping the com- 
putation tractable and still providing enough accuracy for the 
analysis. With a typical 10~ 6 relative and 10~ 15 absolute er- 
rors, the 2D-mixing model with high mixing efficiency has 
about 15 million non-zero Jacobi matrix elements and takes 
about 48 hours of CPU time (Xeon 3 .0 GHz, 4 Gb RAM, gfor- 
tran 4.4-x64) to calculate the disk chemical structure within 
5 Myr. 

Our main set of chemical simulations consists of three runs: 
(1) the laminar disk chemical model (no transport processes 
are taken into account), (2) the fast 2D-mixing model (Sc =1), 
and (3) the slow 2D-mixing transport model (Sc = 100). The 
chemical evolutionary time span is 5 Myr. 

3. WHEN TURBULENCE AFFECTS DISK CHEMISTRY 

Before drilling into complex numerics, it is reasonable to 
analyze general conditions at which turbulent diffusion may 
affect the chemical evolution in a protoplanetary disk. 

3.1. Chemical and dynamical timescales 

In computational fluid dynamics with reacting flows 
Damkohler number Da is often used as a measure of the influ- 
ence of dynamical processes on the chemical evolution. This 
number is simply the ratio of a physical timescale to a chemi- 
cal timescale: 

Da = Tphy S /T c hem- (6) 

When Da < 1, chemical evolution of a molecule is slow and 
therefore may be sensitive to changes in physical conditions 
due to the medium flow. In contrast, when Da 1 , chemical 
evolution is fast and is not affected by transport processes. 
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Figure 3. (From left to right) Distributions of the characteristic timescales: turbulent diffusion, ion-molecule chemistry, photochemistry, freeze-out, desorption, 
and surface chemistry (log 10 scale). 



The physical (or dynamical) timescale is often written as 
Tphy S = L/V, where V is a characteristic value of velocity fiuc- 
tuatio ns set by turbulence and L is their correlation length 
(e.g.. lCant & Mastor akis 2008). For our a-disk model, char- 
acteristic physical timescale is a turbulent mixing timescale: 

T phys (r 1 z) = H(r) 2 /D tmh (r 1 z). (7) 

The distribution of r p h ys (r,z) in the adopted disk model is 
shown in Fig.[3](lst panel). In the outer disk region, r ~ 100— 
800 AU Tphys exceeds 10 5 years, while in the Jovian planet- 
forming zone the mixing timescale is faster, < 10 4 years. Due 
to the vertical temperature gradient this timescale is slightly 
shorter in the superheated disk upper region (see Fig. Q] 
1st panel). Further we compare the mixing timescale with 
timescales of key chemical processes. 

The typical timescale of a first-order kinetics reaction (e.g., 
photodissociation) with a rate k is r c h e m ~ k~ l . For a second- 
order reaction involving species A and B with a rate &ab the 
reaction timescale for the species A is T c h em ~ 1 Aab«b- Lets 
assess the characteristic timescale of ion-molecule chemistry 
in disks, using HCO + as an example. The ion-molecule re- 
actions are rapid even at very cold tempera t ures and usu- 
ally have no barriers ( e.g., ISmith et ail l2004t IWoodall et al.l 
I2007t IWakelamll2009h . We assume that the HCO + evolu- 
tion is governed by the following formation and destruction 
reactions, CO + -> HCO + + H 2 (h = 1.61 10" 9 cm 3 s" 1 ) 
and dissociative recombination HCO + + e" — > CO + H (£2 = 
2.4 10- 7 (r/300)-° 69 cm 3 s" 1 ). Then the corresponding HCO + 
ion-molecule (IM) chemistry timescale can be estimated as 

tim ~ (fcinconH+/nHCO + -^2nHCO + r 1 - (8) 

The IM timescale distribution in the disk is shown in Fig. [3] 
(2nd panel), where we take nco = 6 10~ 5 hh, «h+ = 10~ 10 «h, and 
"hco+ = 10~ 9 «h- As can be clearly seen, ion-molecule chem- 
istry is very rapid, with a typical timescale of < 10 - 10 3 years 
even in low density disk regions. This is also true for neutral- 
neutral reactions without barriers or with small barriers, in- 
volving radicals and open-shell species. They ha ve compara- 
ble tim escales even in the outer cold disk region (Smi th et al.l 
120041) . In the warm (T > 50-100 K) inner disk region 
other neutral-neutral reactions with considerable barriers of > 
1 000 become competitive. Overall, ion-molecule and neutral- 
neutral reactions without large barriers proceed faster than the 
turbulent transport. 

Chemical evolution in upper disk layers is determined by 
photochemistry, which is driven by intense high-energy stel- 
lar (UV, X-rays) and interstellar (UV, CRP) radiation. The 
corresponding timescale is primarily set by rates of the UV 



dissociation and X-rays ionization processes, 

For a CO-like molecule and without shielding, the UV and 
X-ray photorates are = 2 10" 1() (exp" L7A v ^ +exp" L7A v ^ 

and fcpj = 3((x + CcrpX respectively. The calculated photo- 
chemistry timescale is presented in Fig. [3] (3rd panel). Pho- 
tochemistry is fast (< 1 year) in the disk atmosphere (faster 
than the turbulent transport) and becomes slow (> 10 6 years) 
in dense dark disk regions close to the midplane. 

The next important process in disk chemistry is freeze-out 
of neutral species onto dust grain surfaces in cold and dense 
region (T < 20- 120 K). The inferred substantial depletions of 
observed gas-phase molecules in disks compared to the ISM 
are generally interpreted as a combined action of photodis- 
sociation and freez e -out processes (e.g.. iBergin et al.l 12007b 
iDutrev et all l2007at ISemenov etail l2010al) . The adsorption 
(AD) timescale is 

tad ~ k A D = 1 / (7rag r Vthnd) , (10) 

where a gI is the grain radius (cm), Vth is the kinetic veloc- 
ity of molecules (cms -1 ), and «d is the grain concentration 
(cm -3 ). The freeze-out timescale for CO in our disk model 
is depicted in Fig. [3] (4th panel). This value is mostly de- 
termined by density (assuming homogeneous dust and gas 
mixture) and only slightly by temperature and mass of a 
molecule. It varies between < 1 and 10 3 years around mid- 
plane. At higher, less dense disk regions tad is longer, up to 
~ 10 s — 10 6 years, though at such conditions evaporation rate 
will be much shorter. The adsorption timescale is in general 
shorter than the mixing timescale in the disk regions favorable 
for freeze-out. 

A process, competitive to adsorption, is evaporation of icy 
mantles in warm and/or irradiated disk regions. The evap- 
oration timescale in disks is a combination of thermal des- 
orption, CRP-induced desorption, photodesorption and pos- 
sibly other non-ther mal desorption mechanism s (e.g., explo- 
sive desorption; Shalabiea & Greenberg 1994). Thermal des- 
orption will prevail in warm viscously-heated midplane re- 
gion (r ~ 1-5 AU) and across intermediate molecular layer, 
whereas CRP-desorption operates in the coldest, dark outer 
disk midplane, and photodesorption becomes competitive in 
upper disk layer. Thus, evaporation timescale can be written 
as: 

where 

C = lO" 3 ™^. (exp" 2A v x , +exp- 2A v ) , (12) 
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4f = 2.41(r 2 CcRP^oexp-^, 



and 



4= s = ^>exp 



-7/r 



(13) 



(14) 



Here, 7crp = (4.36 1 5 + r 3 ) 1 / 3 is a peak grain temperature 
( K) reached upon a hit b y relativistic iron nucleus (see Eq. 6 
in iSemenov et alj 120041) . i> {) is characteristic vibrational fre- 
quency (s _1 ) of a molecule, and 7 its desorption energy (K). 
In Fig. [3] (5th panel) we show distribution of the evaporation 
timescale for CO in the DM Tau disk model. In cold and 
dark outer midplane, where temperatures are below « 20 K, 
desorption of CO is too slow, and Td es exceeds about 1 Myr, 
whereas in upper disk layer and in the warm planet-forming 
zone desorption is a rapid process. 

Finally, surface chemistry timescale in the absence of tun- 
neling is controlled by thermal hopping rates of the reactants, 
and th e reaction barrier (see Eqs. 12, 14 in ISemenov et al.l 
2010b). As an example, we consider a slow first step in sur- 
face hydrogenation sequence of CO into CH3OH, namely, H 
(ice) + CO (ice) — > HCO (ice), which has a barrier of about 
2500 K caused by bond restructuring. In Fig. [3] (last panel) 
the timescale of the CO surface chemistry is presented. As 
can be clearly seen, surface chemistry has a long timescale of 
~ 10 6 years around midplane and in cold outer disk region, 
which is similar to the turbulent mixing timescale. There- 
fore, chemical species produced mainly via surface processes 
(e.g., ices and complex organics) shall be sensitive to turbu- 
lent diffusion mixing in disks. The aforementioned character- 
istic timescales in the inner (10 AU) and outer (250 AU) disk 
regions are also compared in Tables |3H 

In reality chemical evolution of a large multi-component 
mixture is controlled by both fast and slow processes, so it is 
hard to obtain a single characteristic evolutionary timescale 
as discussed above. Usually detailed analysis of eigenval- 
ues and eigenvectors of Jacobi matrix or its diagonal and off- 
diagonal terms is used to isolate fast an d slow evolving sub - 
systems in the chemical network (e.g.. lLovrics et all 12006). 
However, for a tightly coupled, extended set of chemical re- 
actions, like a typical hydrogen- and carbon-dominated astro- 
chemical network, such an analysis is hard to perform. We 
elaborate a different approach and use the results of our non- 
mixing disk model to derive individual chemical timescales 
as the time needed to reach a chemical quasi steady-state for 
a given species. The steady-state time is the evolutionary 
moment when molecular abundances in the subsequent time- 
steps change by a factor of < 3 or assumed 5 Myr otherwise. 
In Fig. [2] the chemical steady-state times for a few assorted 
species in the DM Tau disk model are presented. 

The steady-state timescale distributions exhibit a complex 
pattern that is not easily related to characteristic chemical 
timescales shown in Fig. [3] Noteworthy, this pattern is simi- 
lar for many species in the network, with slow surface pro- 
cesses, r c hem ~ 10 s — 10 6 years, determining the chemical 
timescales in the midplane, z/H T < 0.8 (here H r is the pres- 
sure scale height). In the warm inner midplane, r < 20 AU, 
the characteristic chemical timescales are large because not 
only H and H2 but also heavy radicals become mobile at 
T ~ 30-40 K, lea ding to active surface chemistry (see, e.g., 
iGarrod et all 120081) . Within ~ 30-50 AU temperatures in 
the midplane are between 15-25 K, and hydrogen evaporates 
too rapidly, while other radicals are too heavy to activate sur- 
face chemistry. This results in a dearth of rapid surface re- 
actions. Therefore, chemical timescales are usually shorter, 



Tchem ~ 1 4 — 1 5 years, in this region, and are governed by 
gas-grain interaction rates. Above the cold and dark mid- 
plane, at w 1 pressure scale height, where a warm molecu- 
lar layer is located, the characteristic chemical timescales are 
shorter, < 10 5 years, since thermal desorption and photodes- 
orption start to prevail there over surface chemistry and ac- 
cretion. Within the inner ~ 200 AU of the molecular layer 
lukewarm temperatures and elevated X-ray ionization enable 
slow photodestruction of well-bound ices. In the upper disk 
region, > 1.5 pressure scale height, where molecules are ef- 
fectively destroyed by high-energy UV and X-ray photons, 
photochemistry coupled to slow neutral-neutral reactions set 
chemical timescales that increase outward from ~ 10 4 years 
till > 1 Myr. The outward decrease of the gas density in this 
region makes recombination rates slower (as these scale down 
as density squared), whereas the stellar UV and X-ray fluxes 
also decrease with the radial distance from the star. Note that 
chemical timescales are shorter in the CO molecular layer 
and, in particular, for C + , which is abundant in the disk atmo- 
sphere at z/H T < 1.8-2, while r c h e m is long and comparable 
to the dynamical timescales in regions where CO2, CH3OH 
ice, and SO2 are produced (see Fig. |2). As we shall see later, 
indeed chemical evolution of these (and many others) species 
is affected by the turbulent diffusion. 

3.2. Mixing Importance Measure 

Considerations presented in the previous subsection pro- 
vide some clues on the mixing sensitivity for various species. 
However, it would be instructive to find some general quan- 
titative measure of this sensitivity. We suggest a new Mixing 
Importance Measure (MIM) to find necessary conditions for 
sensitivity of column densities of a given molecule to the tur- 
bulent mixing, based on results of the laminar model. This 
quantity comes in three varieties, namely, a local value 



MIM,,, : 
a vertically integrated value 



1 n{r,z) 
Da(r,z) N(r) ' 



MIM r = 



,(r) 



and a global value 



MIM: 



1 



(n - r ) 



MTM rz dz, 



MIMrdr, 



(15) 



(16) 



(17) 



where Da(r,z) is the Damkohler number (Eq. |6), N(r) is the 
total vertical column density (cm -2 ) of a given species at the 
radius r, n(r,z) (cm -3 ) is its number density at the disk loca- 
tion (r,z), ZmaxM is the disk height at the radius r, and ro, r\ 
are the disk inner and outer radii, respectively. The charac- 
teristic physical timescale is the diffusion timescale given by 
Eq.|7] and the characteristic chemical timescale of a molecule 
is its quasi steady-state timescale defined above (see Fig. [2]). 
By definition, the MIM for a molecule allows localizing those 
disk regions where the chemical evolution is slow and which 
contribute most to the vertical column density (a typical ob- 
servationally inferred quantity). Thus, the larger the MIM, 
the stronger possible changes in vertical column densities of 
a considered species due to turbulent transport. However, the 
straightforward expression ([TBI is not a sufficient criterion for 
making reliable estimates on a magnitude of such changes. 
For example, in the absence of relative abundance gradients 
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Figure 4. (Top to bottom) Distributions of the Mixing Importance Measure (MIM; left panel) and vertically-integrated MIM (right panel) for C + , CO, CCh, 
CH3OH ice, NH3, and SO9. The higher the MIM, the more sensitive can be the chemical evolution of the given molecule to dynamical processes. Also given are 
globally integrated MIM values. 



even very slow chemistry (resulting in large MIM) leads to 
the same column densities as the non-mixing (laminar) disk 
model. On the other hand, if MIM is small everywhere in 
the disk, the mixing will not alter resulting column densi- 
ties. The local (Eq. [15] ), vertically-integrated (Eq. [Tol l, and 
disk-averaged (Eq.[T7li MIMs for several species are shown in 
Fig.ffl(log 10 scale). 

The MIM distribution for C + peaks at elevated disk heights, 
~ l-2H r , where the C + concentration is high, and photo- 
chemistry is relatively slow yet efficient, especially on ices 
(cf. Fig- [3j. The importance of turbulent mixing to its chem- 
ical evolution is low (lg(M/M) = 0.23) compared to the other 
species, particularly, in outer disk region at r > 300 AU. 
Two other species with similarly low global MIM values of 
\g(MIM) = 0.22 and lg(M/M) = -0.28 are CO and NH 3 , re- 
spectively. The MIM distribution for CO has a distinct peak 
around inner disk midplane at r < 30 AU, where tempera- 
tures are appropriate to its surface conversion to CO2, and 
where the gas-phase CO concentration is high. Another re- 
gion where the CO chemistry may be affected by the turbu- 
lent processes is located at the upper CO molecular layer, at 
~ 2H T , in the inner disk zone. There the UV photodissocia- 
tion of CO is slow due to self-shielding and mutual-shielding 
by H2 and dust, and X-ray ionization of He leads to slow de- 
struction of CO by He + at late times, t ~ 10 5 — 10 6 years. 

In contrast to CO and C + , the MIM for the C0 2 , NH 3 , and 
SO2 chemistry exhibit a somewhat similar pattern (Fig. |4}. 
The most likely dynamically-sensitive regions are the warm 
molecular layer, where photoprocessing is activated, and in- 
ner part of the disk midplane, where surface heavy C-, O- 
and N-bearing radicals become mobile. This is also probably 
true for S-bearing radicals as well, but our surface network in- 
cludes only a very limited sulfur chemistry, therefore the MIM 
for SO2 does not reach its maximum in the inner midplane, 
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Figure 5. Dependence of the ratio of the molecular column densities in- 
tegrated over the disk radius and computed with the 2D-mixing and laminar 
models (1 = 5 Myr) versus MIM. Shown are the molecules with the maximum 
column densities in the laminar case that exceed 10" cirT 2 . The correspond- 
ing linear correlation coefficient is also shown (top left corner). 



and peaks within the SO2 molecular layer. Finally, the MIM 
for CH3OH ice is the largest in the disk midplane where solid 
methanol is produced by surface hydrogenation of CO and by 
CH3 reacting with OH. The disk-averaged MIM values are 
relatively high for CO2 (lg(MIM) = 0.56) and especially for 
S0 2 (lg(M/M) = 0.74) and solid CH 3 OH (\g(MIM) = 0.87). 
Below we will reveal that their column densities are enhanced 
by the turbulent diffusion by more than an order of magnitude, 
unlike those of C + , CO, and NH 3 . 



4. RESULTS OF NUMERICAL MODELING 
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Table S 

Species steadfast to turbulent mixing 



lVTrilfPiilp 






lf>W raax 
'g^'stat 


c + 


0.23 


0.48 


16.79 


c 2 


0.25 


0.68 


14.41 


C 2 H 


0.40 


0.64 


14.29 


C 2 H 2 + 


0.67 


0.68 


12.36 


CH 


0.30 


0.37 


12.91 


CH 2 


0.33 


0.45 


12.73 


CH 3 + 


0.69 


0.57 


12.22 


CH 5 + 


0.34 


0.66 


13.48 


CL 


-0.17 


0.21 


13.95 


CN 


0.12 


0.56 


14.41 


CO 


0.22 


0.43 


18.67 


Fe 


0.31 


0.45 


11.31 


Fe + 


-0.62 


0.19 


13.69 


FeH 


0.22 


0.68 


12.16 


grain(O) 


-0.21 


0.35 


12.29 


H 2 


-1.08 


0.02 


24.30 


H 2 CN+ 


0.13 


0.49 


13.06 


H 2 CO 


0.32 


0.44 


13.91 


H 2 CO+ 


0.78 


0.44 


11.23 


H 2 S 
H 3 + 


0.11 


0.44 


13.23 


-0.68 


0.52 


14.81 


H 3 CO+ 


0.58 


0.37 


11.62 


H 3 0+ 


0.12 


0.64 


13.43 


HCN 


-0.00 


0.65 


14.50 


HCS 


0.27 


0.61 


11.28 


He+ 


0.07 


0.51 


15.07 


HNC 


-0.04 


0.66 


14.24 


Mg 

Mg+ 


-1.67 


0.12 


12.06 


-1.38 


0.10 


13.91 


MgH 2 


-1.38 


0.20 


13.27 


N 2 


0.19 


0.43 


17.61 


Na 


-1.62 


0.10 


11.84 


NH 2 


-0.59 


0.35 


14.18 


NH 4 + 


-0.00 


0.53 


12.17 


OH 


0.09 


0.32 


14.92 


Si 


0.69 


0.42 


13.06 


SiH 


0.49 


0.35 


11.21 


SiH 3 


-0.59 


0.38 


11.05 


C ice 


1.03 


0.50 


16.00 


C 2 H 2 ice 


0.18 


0.62 


13.91 


CH 5 N ice 


1.26 


0.67 


15.54 


FeH ice 


-1.47 


0.02 


16.04 


H ice 


-1.00 


0.37 


18.55 


H 2 ice 


1.15 


0.43 


17.08 


H 2 CS ice 


0.34 


0.52 


16.29 


H 2 ice 


0.06 


0.09 


20.84 


H 2 S ice 


-0.98 


0.05 


17.54 


HCN ice 


0.27 


0.65 


19.43 


HCSi ice 


-0.00 


0.54 


12.50 


HS ice 


1.01 


0.38 


15.54 


MgH 2 ice 


-1.29 


0.01 


16.64 


NaH ice 


-1.78 


0.01 


15.96 


NH 3 ice 


-0.08 


0.11 


19.84 


S ice 


0.96 


0.57 


15.50 


S 2 ice 


0.74 


0.42 


11.45 


SiC ice 


0.88 


0.39 


13.85 


SiCH 2 ice 


0.94 


0.23 


12.95 


SiH 4 ice 


-0.11 


0.08 


16.59 



Note. — Col. 2 The Mixing Importance Measure (MIM) integrated over 
the disk (in log scale). Col. 3 The ratio of the vertical column density at 
t = 5 Myr computed with the 2D-mixing and the laminar chemical models, 
averaged over the radius (in log scale). Col. 4 The maximal vertical column 
density for a given molecule at 5 Myr in the laminar model. 

* Listed are the molecules which vertical column densities in the disk exceed 
10 11 cm" 2 . 




10 100 
R, AU 

Figure 6. Major disk regions are designated by analyzing the distribution 
patterns of the C + and CO abundances. The atmosphere is the disk region 
dominated by C + in contrast to CO, whereas in the molecular layer carbon 
is mainly locked in the CO gas-phase molecules. The midplane is the disk 
region where depletion of gas-phase molecules occurs. 

Table 6 

Species sensitive to turbulent mixing 



Molecule* 


lg(Se) 


lg(CDR) 




c 

V 


O 7Q 


1 80 


17 89 




-0.07 


1 01 


14 35 


rv lj + 


3^ 


1 1 s 
1.10 


1 1 43 


r,H, 
V2rL4 


U.7J 


1 7S 
1 . / J 


1 1 82 


C\N 


0.76 


1.83 


12.53 


*-3 


0.52 


1.12 


14.93 


C3H 


0.64 


1.18 


13.42 




n id 
v. it 


A Q1 


1 ? 




u.ou 


1 OA 


1 J.yj 




U.J7 


1 1 1 
1.11 


1 9 (IS 
1Z.UJ 


p u + 


U.50 


1 f\A 


1 1 CI 
i 1 .63 


L-3IN 


u.oo 


1 Ai\ 
1 .4U 


1 9 rn 

1Z.U1 


C3O 


0.96 


1.25 


11.17 


c 4 


0.73 


1.17 


13.13 


C 4 H 


0.56 


1.11 


13.74 


C 4 H 2 


0.32 


1.09 


13.89 


C 4 H 3 
C 4 H 3 + 


0.56 


1.30 


13.81 


0.59 


1.13 


11.81 


C 4 H 4 


1.04 


1.61 


12.15 


C4H4" 1 " 


0.86 


0.84 


11.59 


C 4 N 


0.84 


1.41 


12.71 


c 5 


0.55 


1.21 


14.18 


C 5 H 


0.72 


1.40 


12.98 


C 5 H+ 


0.96 


1.06 


11.15 


C 5 H 2 


0.69 


1.23 


13.29 


C 5 H 2 + 


1.10 


1.62 


11.07 


C 5 H 2 N+ 


1.14 


1.48 


11.19 


C 5 N 


0.81 


1.63 


12.39 


c 6 


0.62 


1.51 


13.85 


C 6 H 


0.95 


1.59 


12.71 


C 6 H 2+ 


0.91 


1.07 


13.25 


C 6 H 2 + 


1.14 


1.84 


11.10 


C 6 H 6 


1.18 


1.59 


11.39 


c 7 


0.74 


1.39 


13.44 


C 7 H 


1.12 


1.37 


12.22 


C 7 H 2 


1.08 


1.31 


12.36 


c 8 


1.03 


1.47 


13.41 


C 8 H 


1.11 


1.30 


12.31 


C 8 H 2 


1.02 


0.96 


12.93 


C, 


1.09 


1.41 


13.33 


C 9 H 


1.09 


1.11 


11.99 


C 9 H 2 


1.11 


1.46 


12.49 


C 9 N 


1.13 


1.98 


11.05 


CH 2 CN 


0.75 


1.12 


11.80 


CH 2 CO 


1.05 


1.78 


11.29 


CH 2 NH 


0.53 


1.05 


11.56 


CH 3 


-0.35 


0.98 


15.34 


CH 3 CN 


0.76 


1.04 


11.33 


CH 4 


0.60 


0.91 


17.69 
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Table 6 — Continued 



Table 6 — Continued 



1V1U1CL Lilt 






'gV v stat > 


CH 5 N 


0.68 


1.07 


11.50 


CNC+ 


0.33 


1.20 


11.38 


CS 


0.43 


1.91 


12.15 


e~ 


0.42 


0.89 


17.25 


H 


0.93 


1.81 


22.69 


H+ 


0.88 


1.66 


17.07 


H 2 CS 


0.28 


0.93 


12.33 


H 2 NO+ 


0.19 


1.52 


11.20 


H 2 


0.37 


0.95 


15.52 


H 2 S 2 


0.98 


1.67 


11.21 


H 2 SiO 


0.65 


1.68 


11.31 


HC 3 N 


1.03 


1.17 


11.97 


HC 5 N 


1.12 


1.56 


12.15 


HCCN 


0.82 


1.52 


11.25 


HCL 


0.83 


1.62 


11.16 


HCO 


0.66 


1.25 


11.28 


HCO+ 


0.48 


1.35 


13.98 


HNCO 


0.95 


1.96 


12.14 


HNO 


0.18 


1.48 


13.67 


HNSi 


0.47 


1.99 


11.29 


HS 


0.79 


1.37 


11.60 


HS 2 


0.97 


1.67 


11.21 


N 


0.52 


0.84 


17.78 


N 2 H + 


0.47 


1.47 


12.39 


NH 


0.09 


1.13 


13.36 


NH 2 CN 


0.62 


1.07 


12.40 


NH 3 


-0.28 


0.75 


14.74 


NH 3 + 


0.34 


0.89 


11.43 


NO 


0.36 


1.43 


14.57 


NO + 


0.66 


0.89 


11.23 


O 


0.08 


0.73 


18.65 


o + 


0.80 


1.52 


12.24 


o 2 


0.58 


1.98 


16.91 


o 2 + 


0.86 


1.10 


13.41 


OCN 


0.18 


1.18 


13.69 


OCS 


0.49 


1.75 


11.98 


P 


0.26 


0.73 


12.53 


P + 


0.73 


1.56 


12.89 


s 


0.60 


1.23 


14.73 


s + 


0.68 


1.13 


15.52 


Si + 


0.69 


1.25 


14.12 


SiH 4 


0.32 


0.87 


12.07 


so + 


0.59 


1.11 


12.00 


Cjo ice 


1.04 


1.19 


15.00 


C 2 ice 


0.96 


1.91 


15.28 


C 2 H ice 


0.94 


1.91 


15.22 


C 2 H 3 ice 


1.34 


1.89 


12.79 


C 3 ice 


0.96 


1.80 


15.81 


C 3 H ice 


0.78 


1.77 


16.39 


C 3 H 2 ice 


0.70 


1.29 


19.66 


C 3 S ice 


0.92 


1.51 


15.12 


C4 ice 


1.30 


1.70 


12.84 


C4H ice 


1.30 


1.86 


13.18 


C 4 H 2 ice 


0.84 


1.74 


14.43 


C4H 3 ice 


1.28 


1.85 


12.74 


C4H4 ice 


0.79 


1.89 


18.06 


C4N ice 


0.67 


1.75 


14.78 


C4S ice 


1.16 


1.01 


12.84 


C5H 2 ice 


0.77 


1.54 


18.73 


CgH 2 ice 


0.61 


1.34 


18.59 


C6H6 ice 


1.06 


1.99 


13.81 


C7H2 ice 


0.67 


1.62 


17.95 


CgH 2 ice 


0.70 


1.82 


18.44 


CjH 2 ice 


0.58 


1.04 


18.06 


CCL ice 


1.06 


0.98 


11.76 


CH ice 


0.93 


1.20 


16.16 


CH 2 ice 


0.79 


1.37 


17.73 


CH 2 CO ice 


1.36 


1.37 


15.26 


CH 2 NH ice 


1.37 


1.45 


13.37 


CH 2 NH 2 ice 


1.37 


1.76 


12.95 


CH 3 ice 


0.40 


1.34 


17.77 


CH 3 CN ice 


0.71 


1.78 


14.65 


CH 3 NH ice 


1.38 


1.70 


12.95 


CH 3 OH ice 


0.87 


1.70 


16.08 


CH 4 ice 


0.26 


1.51 


19.95 


CL ice 


0.63 


0.88 


15.20 



1 * 

Molecule 






l(rrA7 max ^ 
'g^stat > 


CNice 


0.58 


1.00 


15.53 


CO ice 


0.92 


0.85 


19.51 


CSice 


0.66 


1.87 


14.93 


H 2 CN ice 


1.33 


1.71 


13.25 


H 2 CO ice 


1.03 


0.93 


16.73 


H 2 S 2 ice 


0.63 


1.16 


12.84 


H 2 SiO ice 


1.20 


1.46 


12.60 


HC 2 NC ice 


0.93 


1.47 


12.99 


HC 3 N ice 


1.31 


1.95 


15.57 


HC 5 N ice 


0.71 


1.47 


14.82 


HC 7 N ice 


0.99 


0.94 


14.49 


HC9N ice 


1.04 


0.72 


14.42 


HCS ice 


1.29 


0.72 


11.06 


HNC ice 


0.60 


1.41 


18.32 


HNCO ice 


0.24 


1.10 


14.80 


HNSi ice 


0.81 


1.14 


12.52 


HS 2 ice 


0.64 


1.17 


12.84 


N ice 


1.20 


1.27 


16.13 


N 2 ice 


1.30 


1.70 


14.13 


NaOH ice 


1.39 


0.87 


11.13 


NHice 


0.46 


1.01 


17.17 


NH 2 ice 


0.07 


1.47 


17.59 


NH 2 CHO ice 


1.12 


1.44 


13.23 


NH 2 CN ice 


0.45 


1.26 


13.54 


NO ice 


1.04 


1.56 


15.55 


O ice 


1.04 


1.53 


16.49 


OCN ice 


1.35 


1.63 


12.55 


OH ice 


-0.59 


0.98 


18.05 


P ice 


0.11 


0.94 


14.65 


PH ice 


1.02 


1.93 


13.69 


PH 2 ice 


0.86 


1.65 


13.75 


Si ice 


1.48 


1.49 


12.29 


SiC 2 ice 


1.25 


0.91 


12.76 


SiC 2 H ice 


1.19 


0.93 


12.22 


SiH 2 ice 


1.24 


1.57 


14.06 


SiH 3 ice 


1.24 


1.05 


14.03 


SiN ice 


0.83 


0.74 


13.40 


SiO ice 


1.28 


0.78 


14.99 


SO ice 


0.72 


1.72 


14.86 



Note. — Col.2 The Mixing Importance Measure (MIM) integrated over 
the disk (in log scale). Col. 3 The ratio of the vertical column density at 
t = 5 Myr computed with the 2D-mixing and the laminar chemical models, 
averaged over the radius (in log scale). Col.4 The maximal vertical column 
density for a given molecule at 5 Myr in the laminar model. 
* 

Listed are the molecules which vertical column densities in the disk exceed 
10" cm 4 . 



Table 7 

Species hypersensitive to turbulent mixing 



Molecule* 


lg(Se) 


lg(CDR) 




C10 


1.18 


2.01 


12.09 


C 2 S 


0.71 


2.27 


11.80 


C 3 H 3 


1.12 


2.64 


11.39 


C 3 H 4 


1.03 


3.31 


11.89 


C 3 S 


1.00 


2.31 


11.35 


C 7 N 


1.12 


2.57 


11.58 


CH3C4H 


1.10 


2.18 


11.21 


co 2 


0.56 


2.23 


14.94 


H 2 0+ 


0.95 


2.09 


12.53 


H 2 2 


0.91 


2.53 


13.54 


HC 7 N 


1.13 


2.34 


11.34 


HCOOH 


0.81 


2.07 


12.36 


HSiO+ 


0.75 


2.75 


11.38 


N + 


0.67 


2.00 


11.51 


N 2 


0.44 


2.92 


12.30 


N0 2 


0.85 


2.26 


12.06 


o 3 


0.05 


2.02 


15.06 


OH+ 


1.00 


2.68 


12.86 


SiO 


0.76 


3.33 


13.17 


SO 


0.74 


3.23 


13.33 


so 2 


0.81 


3.83 


12.45 
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Table 7 — Continued 



Molecule* lg(Se) lg(CZ)ff) lg(AO 



C2H4 ice 


0.99 


2.61 


16.71 


C2H6 ice 


0.86 


6.37 


18.21 


C2S ice 


0.95 


2.79 


12.86 


C3H3N ice 


1.14 


2.04 


11.92 


C3H4 ice 


1.00 


5.65 


19.05 


C5 ice 


1.28 


5.24 


12.63 


C 5 H ice 


1.00 


3.67 


13.38 


C5H4 ice 


1.57 


7.39 


14.81 


Cg ice 


1.33 


4.78 


12.51 


C 6 H ice 


0.91 


3.77 


13.08 


CfiHa ice 


1.10 


2.21 


14.51 


C7 ice 


1.39 


5.34 


12.22 


C 7 H ice 


1.11 


3.93 


12.65 


C7H4 ice 


1.53 


8.73 


13.81 


Cg ice 


1.32 


4.65 


12.72 


C 8 H ice 


1.01 


2.83 


13.24 


C 8 H 4 ice 


1.53 


9.24 


14.19 


Cg ice 


1.35 


4.83 


12.79 


C9H ice 


1.14 


3.52 


13.60 


C9H4 ice 


1.54 


10.09 


13.83 


CCP ice 


1.30 


2.74 


11.47 


CH 2 PH ice 


1.54 


6.99 


12.26 


CH3C3N ice 


1.37 


2.18 


13.23 


CH3C4H ice 


1.04 


2.01 


14.79 


CH3C 6 H ice 


1.14 


2.13 


13.66 


CH3CHO ice 


1.17 


3.26 


13.15 


C0 2 ice 


1.03 


3.60 


19.62 


CP ice 


1.12 


2.66 


14.68 


Fe ice 


1.27 


3.14 


13.35 


H 2 C 3 ice 


1.34 


2.96 


15.94 


H2O2 ice 


1.40 


2.20 


16.97 


H 5 C 3 N ice 


1.55 


7.22 


13.16 


HCCP ice 


1.41 


2.85 


11.08 


HCL ice 


0.64 


2.86 


15.60 


HCOOH ice 


1.32 


3.47 


16.71 


HCP ice 


1.21 


3.21 


14.22 


HNC3 ice 


0.91 


2.76 


11.09 


HNO ice 


0.81 


2.38 


17.90 


HPO ice 


1.37 


5.45 


13.28 


Mg ice 


1.66 


7.32 


12.74 


MgH ice 


1.25 


4.37 


14.39 


N2 ice 


0.86 


4.24 


18.37 


N 2 H 2 ice 


0.76 


2.62 


16.65 


Na ice 


1.35 


5.44 


12.99 


N0 2 ice 


1.25 


2.84 


13.20 


NS ice 


1.19 


2.37 


14.27 


O2 ice 


0.99 


2.25 


17.31 


2 H ice 


1.35 


2.79 


11.28 


O3 ice 


0.91 


2.58 


18.07 


PN ice 


1.30 


2.25 


13.69 


PO ice 


1.09 


3.05 


14.17 


SiH ice 


1.44 


2.24 


12.99 


SiS ice 


1.30 


2.14 


11.63 


C3H3 ice 


1.40 


6.40 


11.31 



Note. — Col. 2 The Mixing Importance Measure (MIM) integrated over 
the disk (in log scale). Col. 3 The ratio of the vertical column density at 
t = 5 Myr computed with the 2D-mixing and the laminar chemical models, 
averaged over the radius (in log scale). Col. 4 The maximal vertical column 
density for a given molecule at 5 Myr in the laminar model. 

Listed are the molecules which vertical column densities in the disk exceed 
10" cm- 2 . 

As we have already mentioned, we present results for three 
models. In the laminar model no diffusion is taken into ac- 
count. In the fast diffusion model Sc = 1 is assumed, while in 
the slow diffusion model we use Sc = 100. In Fig.[5]the disk- 
averaged MIMs (Eqs. ITSHTTb are plotted versus the Column 
Density Ratios (CDRs) for fast mixing and laminar cases and 
for all species potentially observable with ALMA (with maxi- 
mum column densities above 10 11 cm" 2 ). The CDR is defined 
as: 

\g(CDR) = max (| lg(^ 2D (r)-lg(^laminar(r)|) (18) 



As we have mentioned in Sect. 13.21 high MIM value is only 
the necessary condition for molecular concentrations to be 
altered by the turbulence. This is clearly demonstrated in 
the Figure [5] where for species with the low MIM values of 
< 0.25 the maximum ratio between the vertical column den- 
sities calculated with the efficient 2D-mixing and the lami- 
nar models does not exceed an order of magnitude. For ex- 
ample, C + , CO, and NH3 belong to this group of species. 
For larger MIMs the scatter in CDR values rapidly increases, 
reaching 10 orders of magnitude. The CH3OH ice, SO2, and 
CO2 falls into this MIM range, with the CDRs in the range 
of 1.7-3.8 dex. The most extreme group of species populat- 
ing the right top corner in Fig. [5] are carbon chains sitting in 
dust grain mantles (e.g., C9H4 ice). The linear Pearson cor- 
relation coefficient calculated for the MIM-CDR distribution 
is « 0.5 - a medium/large confidence correlation. Thus our 
straightforwardly-defined MIM, which is easy to calculate as 
it only requires a non-mixing chemical model, can be useful 
in a priori analysis of a sensitivity of a given chemical species 
to the effects of transport. 

To facilitate further interpretation of our results for the 
reade r, we divide the disk into 3 chemically dist inct regions 
(e.g., lAikawa et alJ 120021 : ISemenov et alJ [2004), using the 
abundance distributions of C + and CO as a proxy (see Fig.[6]l. 
The hot and UV-irradiated atmosphere at z > 1.5— 2H t is 
molecularly poor, with ionized carbon being the primal C- 
bearing species. Deeper in the disk, in the warm molecu- 
lar layer, C + is converted into CO and a multitude of other 
molecules that remain in the gas phase (located between ^0.7 
and ~ 2 pressure scale heights). The dense and cold disk mid- 
plane region is marked by the onset of molecular depletion 
from the gas phase (< 0.5-0.9 scale heights). Also, in what 
follows we refer to the regions located at radii < 100 AU as 
to the "inner" disk regions, and the "outer" disk regions oth- 
erwise. 

The individual families of primal ions, C-, O-, N-, S- 
bearing and complex (organic) species are depicted in Figs. [Tr- 
ill] and discussed in detail below. The major influence of tur- 
bulent transport on the disk chemical evolution can be sum- 
marized as follows. First, two-dimensional mixing behaves 
as a combination of vertical and radial mixing processes, that 
have to be considered simultaneously. Vertical mixing is more 
important as it affects the evolution of gas-phase and surface 
species of any kind, whereas the effect of radial mixing is 
pronounced mostly for the evolution of ices (e.g. the CO2 ice 
from the OH and CO ices). The reason is that radial tem- 
perature gradient is weaker, and thus is only relevant for the 
evolution of polyatomic ices formed via surface reactions of 
heavy radicals. On the other hand, steep vertical gradients of 
temperature and high-energy radiation intensity cause much 
sharper transition from the ice-dominated chemistry in the 
disk midplane to the oasis of the gas-phase chemistry in the 
warm molecular layer. 

Second, the inefficient diffusion in the slow mixing model 
leads to molecular abundances and column densities that al- 
most coincide with those calculated with the laminar (non- 
mixing) model, but not always. The large Schmidt number 
Sc = 100 assumed for this model makes its characteristic dy- 
namical timescale too long, ~ 1 - 100 Myr, to be competitive 
with most of the slowest chemical timescales deterimed by 
surface processes. 

Third, diffusion in the fast mixing model enhances abun- 
dances and column densities of many species by several or- 
ders of magnitude (Tables [5}Q- We divide all the consid- 
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ered species into three groups. Species with fast chemistry 
(faster than the diffusion timescale of ~ 10 5 - 10 6 years) are 
rather insensitive to turbulent transport. Further we call them 
"steadfast" species (Table [5]). These include simple radicals 
and ions (e.g., C + , Mg + , CO, OH, C 2 H, H 3 0\ HCN, N 2 ) 
and few abundant ices (e.g., water and ammonia ices). Their 
column densities calculated with the laminar model and the 
fast mixing model differ by no more than a factor of 3-5, 
which is comparable to intrinsic uncertainties in molecular 
conce ntrations caused by uncertainties in the reaction rates 
fe.g.. lVasyunin et alJl2008h . 

On the other hand, turbulent diffusion alters abundances 
of many polyatomic species, in particular complex (organic) 
molecules and their ices produced (at least partially) on dust 
grain surfaces (e.g., C0 2 , HCOOH, CH 3 CHO ice, carbon 
chain and cyanopolyyne ices, O2 ice, SO, SO2). These species 
are included in "sensitive" and "hypersensitive" groups. Col- 



umn densities of "sensitive" species (Tables [6]) are altered by 
diffusion by up to 2 orders of magnitude, and even stronger 
for a "hypersensitive" group (Table|7}. Mixing steadily trans- 
ports ice-coated grains in warmer regions, allowing more ef- 
ficient surface processing due to enhanced hopping rates of 
heavy radicals. In warm intermediate layer these "rich" ices 
eventually evaporate, and in the inner disk they can also be 
photodissociated by CRP/X-ray-induced UV photons. The 
importance of mixing is higher in an inner, planet-forming 
disk zone, where thermal, density, and high-energy radiation 
gradients are stronger than in the outer region. Finally, even 
efficient 2D transport in the fast mixing model cannot com- 
pletely erase the layered chemical structure in the disk, leav- 
ing the midplane a gas-phase molecular "desert". 

4. 1 . Major atomic and molecular ions and the ionization 



Table 8 

Key chemical processes: ionization degree 



Reaction 


a 


13 


7 


tin in 


mti ax 




[(cm 3 )s-'] 




[K] 


[yr] 


[yr] 


+ 

M2 + X-rays — ^ e + H2 


u.yjc,x 




u 


1 on 
1 .uu 


1 KM A 

r .zu ^4 


rl2 + ^.l\..r\ — t e + rl2 




u 


u 


1 on 
1 .uu 


D .UU 


H + X-rays — y H + + e 


U.40C,X 




u 


1 nn 
1 .uu 


J.UU^O 


C + + X-rays — > C ++ + e~ 


5i"v 

v. J^x 





() 


1.00 


5.00(6 


S + + X-rays -> S ++ + e~ 


Cx 








1.00 


5.00(6 


H 2 + X-rays — > e~ + H + H + 


0.02Cx 








1.00 


5.00(6 


C + UV -> e- + C+ 


0.31 (-9) 





3.33 


1.00 


5.00(6 


H 2 CO + UV — > HCO+ + H + e~ 


0.48 (-9) 





3.21 


4.18(2) 


1.20(4 


H 2 + + UV ->• H+ + H 


0.57 (-9) 





2.37 


1.00 


55.90 


H 3 + + UV -> H + + H 2 


0.20 (-7) 





1.80 


1.00 


5.00(6 


H+ + O -» 0+ + H 


0.70 (-9) 





2.32(2) 


1.00 


5.00(6 


H+ + 2 -> 2 + + H 
H + + S — > S + + H 


0.12 (-8) 








1.00 


5.00(6 


0.1 3 (-8) 








1.00 


5.00(6 


He+ + CO -> C+ + + He 


0.1 6 (-8) 








1.00 


5.00(6 


C+ + H 2 -> CH 2 + 


0.40(-15) 


-0.20 





1.00 


5.00(6 


S+ + H 2 -> H 2 S + 


1.00(-17) 


-0.20 





1.00 


5.00(6 


H 2 + + H 2 -> H 3 + + H 


0.21 (-8) 








1.00 


5.00(6 


H 2 + + H2-s>H 3 + + H 
N 2 + + H 2 -> N 2 H+ + H 


0.61 (-9) 








1.00 


5.00(6 


0.17 (-8) 








1.00 


5.00(6 


NH+ + H 2 -> H 3 + + N 


0.23 (-9) 








1.96 


5.00(6 


C + + S ->• S+ + C 


0.15 (-8) 








1.00 


5.00(6 


H 3 + + CO ->• HCO+ + H 2 


0.1 6 (-8) 








1.00 


5.00(6 


H 3 + + No -> N 2 H+ + Ho 
HS + + H — > S + + H 2 


0.17 (-8) 








1.00 


5.00(6 


0.11 (-9) 








1.00 


5.00(6 


N 2 H+ + CO -> HCO+ + Ni 
NH 2 + + N ->• NiH+ + H 


0.88 (-9) 








1.00 


5.00(6 


0.91 (-10) 








1.00 


5.00(6 


0+ + H -> H+ + O 


0.70 (-9) 








7.48 


5.00(6 


0H+ + N 2 -> N 2 H+ + O 


0.36 (-9) 








2.34(4) 


5.00(6 


S+ + 2 -> S0+ + 


0.20(-10) 








1.00 


5.00(6 


S+ + OH -» SO+ + H 


0.46 (-8) 


-0.50 





1.00 


5.00(6 


S + + NH 3 -> NH 3 + + S 


0.1 6 (-8) 


-0.50 





1.00 


5.00(6 


S + + e~ -> S 


0.39 (-11) 


-0.63 





4.18(2) 


5.00(6 


HCO+ + e " -> CO + H 


0.24 (-6) 


-0.69 





1.00 


5.00(6 


N 2 H+ + e- -> NH + N 


0.64 (-7) 


-0.51 





1.00 


5.00(6 


H 3 + + e->OH + H + H 


0.26 (-6) 


-0.50 





1.00 


5.00(6 


e~ + grain(O) — > grain(-) 


1.00 








1.00 


5.00(6 


HCO + + grain(-) ->• CO + H + grain(O) 


1.00 








1.00 


5.00(6 



The magnetorotational instability has been shown to be the 
most promising mechanism in driving turbulence in weakly 
ionized PPDs (Balb us & Hawlevl[T9 91). However, even high- 
energy cosmic ray particles cannot penetrate in cold and 
dense disk interiors (Umebayashi & Nakanol [l981l) . and the 
ionization degree drops so low that MRI may not be oper- 



ational at all ("dead zone", e.g.. [Gammi el [19961: ISano etalJ 
I2000t lArmitage et al.ll2003b lllgner & Nelsonll2008l) . At these 
conditions slow radiative recombination of metallic ions be- 
comes important for the evolution of ioniz ation fraction, 
if these are still present i n the gas (e.g., Fromang et al. 
12001 ISemenov et aT]l2004l lllgner & Nelsonl l2006al) . ~On 
the other hand, the viscous stresses can be transported to 
the "dead zone" from above MRI-active accretion layers by 




Figure 7. Abundances and column densities of electrons and selected ions in the DM Tau disk at 5 Myr. (Top to bottom) The log of relative abundances (with 
respect to the total number of hydrogen nuclei) and vertical column densities of electrons, H + , C + , HCO + , H30 + , N2H + , Hj, and S + . (Left to right) Results for 
the three disk models are shown: (1) laminar chemistry, 2) the 2D-mixing chemistry (5c = 100), and (3) the 2D-mixing chemistry with Sc = 1. 



diffusive mixin g, maintaining a low level of accretion in 
this region (e.g.. Fleming & StondfcOOl IWunsch et aLll2006t 
Turner et al. 2007). In this Section we analyze in detail chem- 
ical processes responsible for the evolution of the ionization 
state of protoplanetary disks and how these are altered by tur- 
bulence. 

Turbulent diffusion does not affect column densities of 
many atomic ions and a few simple molecular ions, e.g. C + , 
Mg\ Fe + , He+, H+, CH3 + , NHJ (by less than a factor of 3, 
see Tabled. The charged species sensitive to mixing include, 
e.g., hydrocarbons, electrons, H + , O" 1 ", S + , N2H" 1 ", and HCO + 
(their column densities are altered by factors of ~ 3-50; Ta- 
ble^. Only a few ions are hypersensitive to mixing, e.g. N + , 
OH\ H z O + (CDRs are - 100; Table E). This is a manifes- 
tation of the fact that the global evolution of the ionization 
degree is partly affected by slo w surface chemistry (recombi- 
nation, dissociation) (see also lSemenov et al. 20QJ). 

In Fig.|7]the distributions of molecular abundances and col- 
umn densities at 5 Myr of several major ions and the disk 
ionization degree calculated with the laminar and the mixing 
models are shown. The global ionization structure of the disk 
shows a layered structure similar to that of Photon-Dominated 
Regions (see also, e.g. , ISemenov et al.l 12004; Bergi n et al.1 
I2007t iRollig et al.ll2007h : (1) heavily irradiated and ionized, 
hot atmosphere where the key ions are C + and H + , (2) partly 
UV-shielded, warm molecular layer where carbon is locked 
in CO and major charged species are X-ray-produced H + and 
polyatomic ions like HCO + and H3, and (3) dark, dense and 
cold midplane where most of molecules are frozen out onto 
dust grains and the most abundant charged species are dust 
grains and H|. 

Turbulent diffusion lowers abundances of electrons and 
atomic ions such as C + , S + , and H + by up to 3 orders of 



magnitude, mainly in the inner molecular layer (r < 100 AU, 
~ 1 -2H r ), see Fig. [7] Consequently, their column densities 
within ^100 AU are decreased by factors 3 - > 100, and 
much less in the outer disk. The effect is exactly opposite for 
H3O 1 ", HCO + , and N2H + , which column densities and abun- 
dances are enhanced by the mixing in this region, albeit only 
up to an order of magnitude (Fig. |7}. The mixing remark- 
ably expands their inner molecular layers that are rather con- 
fined in the laminar case (particularly, for N2H" 1 "). The abun- 
dance distribution of H3 shows similar expansion due to the 
mixing as the other polyatomic ions, while its column den- 
sity is slightly decreased at r < 30 AU. Overall, turbulent stir- 
ring makes chemical gradients of key ionized species stronger 
compared to the laminar disk model. On the other hand, the 
slow mixing model does not differ significantly from the lam- 
inar model (compare 1st and 2nd panels in Fig. [7}. 

To better understand these results, we performed detailed 
analysis of the chemical evolution of the dominant ions and 
the ionization degree shown in Fig. [7] in two disk vertical 
slices at r = 10 and 250 AU (the laminar chemical model). The 
most important reactions responsible for the time-dependent 
net change of the electron, H\ C + , HCO\ H 3 + , N 2 H\ HJ, 
and S + concentrations in the midplane, molecular layer, and 
atmosphere are presented in Table [8] both for the inner and 
outer disk regions. Time interval when these processes are 
active is also shown. Note that our list contains only top 15 re- 
actions for the selected species per region (midplane, molecu- 
lar layer, atmosphere) for the entire 5 Myr time span, with all 
repetitions removed. 

The energy deposition into the disk matter leads to ioniza- 
tion of molecular and atomic hydrogen by high-energy cos- 
mic ray particles and the stellar X-ray radiation. Typical X- 
ray ionization rates in the inner and outer disk atmosphere 
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are (x ~ 10~ 12 s" 1 and w 2 10~ 13 s" 1 , respectively. The H£ is 
rapidly converted to H3 by collisions with H2, while in the 
disk atmosphere both H| and H3 can be dissociated by UV, 
producing H + . The line UV-photodissociation of molecular 
hydrogen is hindered by strong self-shielding of this mole cule 
( Solomon & Wickramasinghe 1969; van Dishoeck 1988). In 
the atmosphere, the X-rays ionize other neutral atoms (e.g., 
He), and atomic ions such as C + and S + , leading to abun- 
dant doubl y-charged atomic ions in the very i nner region 
(see, e.g. lMeijerink et al.l2008l;lNajita et al.l2010l) . The heavy 
atomic ions are also produced by charge transfer reactions 
with H + and C + . The atomic ions further slowly radiatively 
recombine with electrons (e.g., S + + e" — > S) or converted 
to polyatomic ions by radiative recombination reactions (e.g., 
C + + H2 — > CHJ) or via fast ion-molecule pathways involv- 
ing simple radicals (e.g., S + + OH — > SO + + H ). As a re- 
sult, the disk fractional ionization in the atmosphere is high, 
X(e~) ss 1CT 4 , and dominant charged species are C + (with the 
abundance < 10" 4 ) and H + (~ 10" 5 - 10" 3 ). The pace of their 
chemical evolution is restricted to sl ow radiative recombina - 
tion and association processes (e.g., Wakel am et al.l l2010b), 
Tchem ~ 10 3 - 10 4 yrs. These timescales are comparable to the 
vertical mixing timescales (r p h ys ~ 10 3 yrs) only in the inner 
disk region (see Figs. [3] [2] and Tabled. 

In molecular layer located beneath the disk atmosphere the 
ionization chemistry is mostly regulated by fast ion-molecule 
reactions (see Fig. [3] 2nd panel), e.g. H3 + N2 — > N2H + 
+ H2. Typical X-ray ionization rates in the inner and outer 
molecular layer are £x w 10~ 12 s" 1 and w 310~ 14 s , re- 
spectively. The partly absorbed X-rays ionize H2, produc- 
ing Hj and H + . In the inner disk molecular layer, the X- 
ray-driven chemistry destroys 99% of molecular hydrogen 
(r < 100 AU). There H2 is converted to atomic hydrogen via 
fast ion-molecule reactions involving OH + and H2O" 1 " (e.g., 
H 2 + + H 2 H 3 + +H; Table S). In turn, these ions be- 
come abundant due to steady release of elemental oxygen 
from water, which begins from photodissociation of the wa- 
ter ice. The H£ is rapidly converted to H3 that protonates 
abundant radicals (e.g., CO, N2), forming various polyatomic 
ions (e.g., HCO + , N2H + ). Thus, the abundance distributions 
of the polyatomic ions are peaked along a rather narrow stripe, 
0.5 H r , where both H3 and progenitor molecules are abun- 
dant. The molecular ion layers are vertically expanded in the 
outer disk at r > 200 AU, where the surface density drops 
so low that the disk becomes partly transparent to the inter- 
stellar UV radiation. The N2H + is destroyed by reactions 
with abundant CO and, thence, the N2H 1 " layer is particu- 
larly thin (see Fig. |7). The N 2 H + , HCO + , and H 3 + lay- 
ers are rather sharp-edged from below where their parental 
molecules ar e locked in dust icy mantles (e.g., Oberg et al. 
120071 l2009bllal) . whereas the H3 layer ex tends more toward the 
cold region as H2 does not deplete (e.g..lSandford et alJlT993t 
iBuch & Devlinlll994t iKristensen et alj|201 lh . The formation 
of polyatomi c ions is balanced out by dissociative recombina- 
tion with e" (Florescu-Mitchell & Mitchell 2006), and within 
« 10 4 years (inner disk) and w 10 s years (outer disk) a chem- 
ical quasi steady-state is reached. This equilibrium timescale 
is mainly set by slow X-ray-driven destruction of H2 that is 
needed to produce H3, rx ~ 2 10 3 - 10 s years. Thus the char- 
acteristic chemical timescales exceed those due to the mixing, 



making the ionization state of the inner molecular layer to be 
vulnerable to the turbulent transport. The dominant charged 
species in the warm molecular layer are H + as well HCO + , 
H3O 1 " and C + , with abundances of 10~ 10 - 10~ 5 . The resulting 
ionization fraction is Xe~ ~ 10~ 8 - 10~ s , which is s ufficient for 
MRI to be operational (e.g., Fromang et al. 2002). 

Finally, in the cold, dark midplane the X-ray ionization 
rates further decrease (£x ~ 10~ 19 - 10~ 15 s" 1 ), and the cos- 
mic ray ionization becomes important ((cr = 1.3 10~ 17 s" 1 ). 
Except of H2, majority of molecules severely deplete within 
< 10 3 years, leaving only a tiny fraction in the gas (see Fig. [3] 
4th panel, and Fig. [8). Electron sticking to grains is efficient, 
and negatively charg ed grains start playing important role as 
charge carriers (e.g.. lUmebavashi & Nakanolll980l lOkuzumil 
2009; Perez-Becker & Chiang 201 1 ) and sites of dissociative 
recombination (e.g., HCO + + grain(-) — > CO + H + grain(O); 
see Table [8). Consequently, the dominant charged species in 
the midplane are dust grains and Ht, with H + being important 
in the upper midplane layer (Fig. |7). 

The fractional ionization drops to as low values as < 
10~ 12 -10~ 9 , and a "dea d" zone for the MRI-driven accre- 
tion may develop (e.g., iGammid [19961 ; iTurner etafl 120071: 
iDzvurkevich et al] 12010). Characteristic timescales of the 
evolution of the fraction ionization in the disk midplane are ~ 
5 10 3 - 10 5 years, after which a quasi steady-state is reached. 
These values are shorter than the mixing timescale in the mid- 
plane (Fig. [3}, and influence of turbulent transport on ioniza- 
tion chemistry shall be negligible. 

Indeed, as can be clearly seen in Fig. [7] (compare 1st and 
3rd panels), turbulent diffusion most strongly affects abun- 
dances of ions in the inner warm molecular layer, at r < 
200 AU, H r w 1-2, where the characteristic diffusive mix- 
ing timescale (~ 10 3 years) is shorther than the chemical 
timescale (~ 10 4 years) regulated by the slow X-ray-driven 
destruction of H2 and radiative recombination and associa- 
tion processes involving C + . The diffusion timescale is rela- 
tively short in the inner disk because of stronger gradients of 
density and ionization rates, elevated temperatures, and short 
distances for transport (see Figs. Q~|and[3). The vertical up- 
take of molecular hydrogen from the inner midplane to the 
inner intermediate layer allows to produce more H3, while 
simultaneously lowering ionization fraction and abundances 
of H + and C + due to their enhanced neutralization in via ra- 
diative association, charge exchange, and ion-molecule reac- 
tions (e.g., C + + H2 —> CHJ; Table [8}. The increased abun- 
dances of H3 in the inner molecular layer leads to more ef- 
ficient production of polyatomic ions through protonation re- 
actions, which lead to remarkable peaks in the relative abun- 
dances of HCO + , N2H + , H30 + , and other molecular ions. The 
chemistry of OH + , H20 + , and N + ions is most sensitive to the 
X-ray-drive n, partly surface, p rocesses (similar results were 
obtained by lAresu et aTll201 II) . Naturally, OH + , H 2 + , and 
N + are most sensitive to the mixing in our model. 

To summarize, turbulent mixing lowers the disk ionization 
fraction and the abundances of atomic ions in the inner disk 
region at intermediate heights (r < 200 AU, ~ 1 - 2H r ) by 
up to 3 orders of magnitude, whereas the concentrations of 
polyatomic ions are enhanced there. 

4.2. Carbon-containing molecules 



16 



Semenov & Wiebe 




10 100 10 100 10 100 10 100 10 100 10 100 10 100 10 100 

R, AU R, AU 



Figure 8. The same as in Fig.[7]but for the C-containing species. Results are shown for CH3 + , CH4, C2H, C2H2, Qjlrb, C2H2 ice, CgHj ice, and C8H2 ice. 
Abundance distribution patterns of many of the carbon chains look similar to that of C 8 H2- 

Table 9 

Key chemical processes: C-bearing species 



Reaction 


a 


P 


7 


tmin 


knax 




[(cm 3 )s-'] 




[K] 


[yr] 


[yr] 


CgH6 ice + hi'cRP — > C6H4 ice + tb ice 


3.00(3) 








1.00 


5.00(6 


C 8 H2 ice + hfcRP — > C 8 H ice + H ice 


1.75(3) 








1.00 


5.00(6 


CgtLt ice + hi/cRp — > CgH2 ice + tb ice 


7.50(3) 








1.00 


5.00(6 


H 2 C 3 ice + hi/ C Rp — > C 2 H 2 ice + CO ice 


1.80(3) 








1.00 


8.18(2 


C2H2 ice + UV — > C2 ice + H ice + H ice 


0.66 (-10) 








1.00 


5.00(6 


CgH6 ice + UV — > C 6 H 4 ice + H2 ice 


1.00 (-9) 





1.70 


1.00 


5.00(6 


C 8 H 2 ice + UV -> C 8 H ice + H ice 


1.00 (-9) 





1.70 


1.00 


5.00(6 


CH 3 + UV — > CH 3 + + e" 


1.00 (-10) 





2.10 


1.00 


5.00(6 


C 2 H + UV — > C? + H 


0.52 (-9) 





2.30 


1.00 


5.00(6 


CoH 2 + UV -> C 2 H + H 


0.33 (-8) 





2.27 


1.00 


5.00(6 


C 2 H 3 + UV -> C 2 H 2 + H 


1.00 (-9) 





1.70 


1.00 


1.75(5 


CoH 4 + UV -> C 2 H 2 + H 2 


0.30 (-8) 





2.10 


2.14(2) 


1.75(5 


C 8 H 2 + UV -> C 8 H + H 


1.00 (-9) 





1.70 


1.00 


5.00(6 


CH 4 + UV — > CH 2 + H 2 


0.84 (-10) 





2.59 


1.00 


5.00(6 


C 2 H + grain — > C 2 H ice 


1.00 








1.00 


5.00(6 


C6H6 + grain — > CgHg ice 


1.00 








1.00 


5.00(6 


C2H2 + grain — > C2H2 ice 


1.00 








1.00 


5.00(6 


C 8 H2 + grain — > C 8 H2 ice 


1.00 








1.00 


5.00(6 


CH4 + grain — > CH4 ice 


1.00 








1.00 


5.00(6 


C 2 H 2 ice -> C 2 H 2 


1.00 





2.59(3) 


1.00 


5.00(6 


C 8 H 2 ice -> C 8 H 2 


1.00 





7.39(3) 


1.00 


5.00(6 


CH 4 ice -» CH4 


1.00 





1.30(3) 


1.00 


5.00(6 


H ice + CH 3 ice — Y CH4 ice 


1.00 








1.00 


5.00(6 


H ice + CH 3 ice — > CH 4 


1.00 








1.00 


5.00(6 


H ice + C2 ice — > C 2 H 


1.00 








1.00 


5.00(6 


H ice + C2H ice — ¥ C2H2 ice 


1.00 








1.00 


5.00(6 


H ice + C2H ice -> C2H2 


1.00 








1.00 


5.00(6 


H ice + C 8 H ice — > C 8 H 2 ice 


1.00 








1.00 


5.00(6 


H ice + C 8 H ice C 8 H2 


1.00 








1.00 


5.00(6 


C 8 H2 ice + H ice — ¥ C 8 H 3 ice 


1.00 





1.21(3) 


1.00 


5.00(6 


CH 3 + + H 2 -> CH 5 + 


0.13 (-13) 


-1.00 





1.00 


5.00(6 


CH 2 + + H 2 -»• CH 3 + + H 


0.12(-8) 








1.00 


5.00(6 


CH 5 + + C -> CH 4 + CH+ 
C + + C 2 H 6 -> CH 4 + C 2 H 2 + 
CH 5 + + CO ->• CH 4 + HCO+ 
CH 4 + H+ -> CH 3 + + H 2 


1.00 (-9) 








1.00 


5.00(6 


0.17 (-9) 








4.18(2) 


5.00(6 


0.32 (-9) 


-0.50 





1.00 


5.00(6 


0.23 (-8) 








1.00 


5.00(6 


+ CH 3 -> CH 4 + e - 


1.00 (-9) 








1.09(2) 


5.00(6 


C + C 8 H. -> C 9 H + H 


0.90 (-9) 








1.00 


5.00(6 
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Table 9 

Key chemical processes: C-bearing species 



f~1 , /it t , in TT i T T 

C + CH3 — > C2H2 + H 


1 An / 1 n\ 

1.00 (— 10) 








1.00 


5.00(6) 


f 1 /"'"LJ T-J . 1 T_| 1 O 

L. 4- L-ri^L-gri — ^ ^8^2 + rb 


0.74 (—9) 








1 .00 


1.31 (6) 


C 2 H + -» CO + CH 


U. 1 / (— 1U) 


U 


ft 

u 


1 nn 

1 .UU 


^ nn (& \ 


f\ r I , / 1 . ( ' [ i . (V| 

l~3xl + \J — ? L-2^l ' *~ 


n 1 7 f 1 n\ 
U. 1 / (— 1U) 


u 


u 


1 nn 


J.oZ 


L^3rl2 ■+■ — r ^2^2 ^-AJ 


n ^ n f 1 n 1 ! 


n ^n 

U. jU 


u 


1 nn 
1 .uu 


J.UU (0 ) 


C ' II 1 ' ( • I 1 . /" 1 - 

L-2rl2 + L — ? ^2^2 ^ *- 


U. 1 1 (— 0) 


U 


U 


1 nn 
1 .uu 


^ nn /^.t 


('II / | . c " it ni'n 1 

^9x13 T W — r L- 8 rl2 + rll^U 


n on f o\ 


u 


u 


1 nn 
1 .uu 




CiH2 + + e~ — > C2H + H 


0.29 (-6) 


-0.50 





1.00 


5.00(6) 


C 2 H 3 + + e" -> C 2 H 2 + H 


0.14(-6) 


-0.84 





1.00 


5.00(6) 


C,H 3 + + e ^C 2 H + H + H 


0.30 (-6) 


-0.84 





1.00 


5.00(6) 


C 3 H 2 + + e" -> C 2 H 2 + C 


0.30 (-7) 


-0.50 





1.00 


1.20(4) 


C3H 2 N+ + e" C 2 H 2 + CN 


0.23 (-6) 


-0.50 





1.00 


5.00(6) 


C 3 H 2 N + + e~ -> C 2 H + HNC 
C 8 H 3 + + e-^C 8 H2+H 


0.75 (-7) 


-0.50 





55.90 


5.00(6) 


1.00 (-6) 


-0.30 





1.00 


5.00(6) 


C 8 H4 + + e-^C 8 H2 + H 2 


1.00 (-6) 


-0.30 





1.00 


5.00(6) 


CH9CO+ + e" -» C9H2 + 


0.20 (-6) 


-0.50 





1.00 


5.00(6) 


CH 5 + + e- -> CH 4 + H 


0.14(-7) 


-0.52 





1.00 


5.00(6) 


C 2 H 3 + + grain(-) -> C 2 H + H + H + grain(O) 


0.59 








7.48 


5.00(6) 



In this subsection we analyze in detail chemical and mixing 
processes responsible for the evolution of the hydrocarbons in 
protoplanetary disks. 

There are only several steadfast light hydrocarbons, includ- 
ing 4 neutral diatomic and triatomic species (e.g., CH, CH2, 
C 2 H), 4 ions (e.g., CH^, CU + 5 ), and the C 2 H 2 ice. Their col- 
umn densities in the laminar model and the fast mixing model 
differ by a factorof < 3. In contrast, majority of hydrocarbons 
are sensitive to the turbulent transport that alters their column 
densities by up to 2 orders of magnitude. The turbulence- 
sensitive hydrocarbons include 28 neutral molecules (e.g., 
CH 3 , CH 4 , C 2 H 2 , C 3 H 2 , C 6 H 6 ,..., C 9 H 2 ), 8 ions (e.g.,, C 2 H+, 
C 3 H,, C 6 H£), and 22 ices (CH 4 , C 2 , C 2 H, C 2 H 3 , C 9 H 2 ; 
58 species in total). Diffusive mixing modifies their column 
densities by up to 2 orders of magnitude. Finally, for 25 hy- 
persensitive hydrocarbons the turbulent diffusion changes the 
column densities by more than 2 orders of magnitude (up to 
a factor of 10 10 for the C9H4 ice). These species include no 
ions, 4 neutral hydrocarbons (e.g., C10, C3H3, CH3C4H), and 
21 ices (C 2 H 6 , C3H4, C 5 , C5H4,..., C9H, C9H4). A rough trend 
is that heavier and more saturated hydrocarbons are stronger 
affected by the turbulent transport than lighter species, and 
that abundances and column densities of hydrocarbon ices 
are stronger altered compared to their gas-phase counterparts. 
Overall, frozen hydrocarbons are among the most sensitive 
species to the mixing in our chemical model. Unlike the 
chemistry of the primal ions and the ionization degree dis- 
cussed in the previous subsection, the chemical evolution of 
hydrocarbons should be dominated by slow surface hydro- 
genation and radical -radical growth processes, and slow evap- 
oration. 

In Fig. [8] distributions of relative abundances and column 
densities at 5 Myr of CH}, CH 4 , C 2 H, C 2 H 2 , C 8 2H 2 , C 2 H 2 
ice, CeH6 ice and C 8 H 2 ice calculated with the laminar and 
mixing models are presented. The hydrocarbon abundances 
show a 3-layered structure similar to that of the polyatomic 
ions, with rather narrow molecular layers of 0.2-0 .5 H r at 
z ~ 1 H r and typical values of ~ 10~ 10 - 10~ 7 . Note that rela- 
tive abundances of hydrocarbon ices are also high in the warm 
molecular layer (X ~ 10~ 10 - 10~ 7 ). The hydrocarbons are 
fragile to the UV and X-ray irradiation and thus are absent in 
the disk atmosphere, apart from the photostable ethynyl rad- 



ical (C 2 H) and CH3 ion that are abundant in the lower disk 
atmosphere, at z ~ 2H r . In the midplane the gas-phase hy- 
drocarbons are depleted, whereas their ices are moderately 
abundant. Turbulent mixing expands and enhances molecu- 
lar layers of gas-phase hydrocarbons at almost all radii. For 
the ices there is no such a clear trend. As a rule, abundances 
of carbon chains that serve as intermediate products to sur- 
face hydrogenation are lowered by the turbulent transport, e.g. 
the C 8 H 2 ice, whereas more saturated hydrocarbons show in- 
creased abundances, e.g. the C^He ice (Fig. [HJ. Note that 
the turbulent transport affects gas-phase and solid abundances 
of heavy hydrocarbons (like C 8 H 2 ) in an opposite way, en- 
hancing concentration of a gas-phase species and reducing its 
solid-state abundances. We attribute such a behavior to partic- 
ularly slow evaporation of frozen heavy carbon chains that is 
comparable or longer than the dynamical timescale, < 1 Myr, 
even in the warm molecular layer. 

To better understand these results, we performed detailed 
analysis of the chemical evolution of the hydrocarbons, shown 
in Fig. [8] in two disk vertical slices at r = 10 and 250 AU 
(the laminar chemical model). The most important reactions 
responsible for the time-dependent evolution of their abun- 
dances in the midplane, the molecular layer, and the atmo- 
sphere are presented in Table both for the inner and outer 
disk regions. The final list contains only top 25 reactions for 
CH+, CH 4 , C 2 H, C 2 H 2 , C 8 H 2 , C 2 H 2 ice, C 6 H 6 ice and C 8 H 2 
ice per region (midplane, molecular layer, atmosphere) for the 
entire 5 Myr time span, with all repetitions removed. 

The chemical evolution starts with production of light hy- 
drocarbons in the gas phase by radiative association of C + 
with H 2 , followed by hydrogen addition reactions: C + — > CH£ 
— >■ CH3 — S- CH5. An alternative route is the radiative associ- 
ation reaction between C and H 2 , leading to CH 2 , that can 
be further converted to C 2 H by addition of C. The protonated 
methane reacts with electrons, O, CO, C, OH, etc., forming 
CH4. Methane undergo reactive collisions with C + , produc- 
ing C 2 Hj and C 2 Hj. All these primal hydrocarbons disso- 
ciatively recombine on electrons or negatively charged grains 
(in the inner dark midplane), which leads to simple neutral 
species like CH, CH3, C 2 H, and C 2 H 2 . The neutral hydro- 
carbons readily react with ionized and neutral atomic carbon, 
forming other, heavier hydrocarbons, e.g. C 8 H 2 + C -> C9H + 
H and C 9 HJ + O -> C 8 H 2 + HCO + . The growth rate of chem- 
ical complexity of carbon chains is regulated by ion-molecule 



18 



Semenov & Wiebe 



and neutral-neutral exothermic reactions with atomic oxygen 
that lead to formation of CO, e.g. O + C3H2 — > C2H2 + 
CO (see also Fig 7 in iHenning et alj d2010h and Fig. 4 in 
Turner et al. (2000)). The associated chemical timescales are 
< 10 3 - 10 4 years, see Fig. [3]and Tables [3] and [4] but not ev- 
erywhere in the disk. In the disk atmosphere, at z ~ 2H r , 
the chemical evolution of CH2 and its daughter species, C2H 
and CH3, is subject to the evolution of H2 and O. As we dis- 
cussed in the previous subsection, at these high altitudes slow 
destruction of molecular hydrogen, water, and CO by He + lo- 
cally sets long evolutionary timescales (< 10 5 years) for many 
species, e.g. ethynyl, CH2, and CH3 radicals. 

To enhance further production of hydrocarbons in the gas, 
particular physical conditions have to be reached. The heavy 
carbon-bearing compounds can be produced at high densities 
and T > 800 K by p yrolysis of precursor hydrocarbons (e.g., 
Morgan et"aT] [T99l or in a Fischer- Trop sch-like process in- 
volvin g CO, H2, and catalytic surfaces (Tscharn uter & GalU 
120071) . These conditions are met at sub-AU radii in the very 
inner disk midplane that is not examined in the present study. 
Another possibility for the gas-phase production of hydro- 
carbons is to bring elemental carbon locked in CO back to 
the gas, while simultaneously maintaining low oxygen abun- 
dances. Due to self-shielding and mutual-shielding by H2, 
the UV photodissociation of CO is only important in di- 
lute di sk atmosphere (e.g., Ivan Dishoeckl [19881: IVisser et all 
2009a). However, in the lower part of the inner warm molec- 
ular layer (r < 200 AU, z ~ 1 H r ) the He + ions produced by the 
X-ray ionization slowly destroy CO, restoring gas-phase con- 
centrations of C + and O (tx ~ 5 10 3 years). Relatively high 
densities (n(H) ~ 10 8 — 10 9 cm" 3 ) and lukewarm temperatures 
(T ~ 30-75 K) in this region (Fig. [T]) allow rapid conver- 
sion of releazed oxygen into water. Thus, locally in the inner 
regions of the X-ray-irradiated disks the gas-phase C/O ratio 
may revert from the Solar value of ~ 0.43 to a > 1 value typ- 
ical of carbon-rich AGB shells (D. Hollenbach, priv. comm.), 
which facilitates gas-phase formation of heavy hydrocarbons. 

An alternative efficient route to accumulate complex hy- 
drocarbons in disks is surface chemistry coupled to the de- 
struction of hydrocarbon ices and other C-bearing species. 
In cold disk regions (T < 20 K) it is mainly hydrogenation 
of precursor species like C„H m (« = 2-9, m < 2). At low 
dust temperatur es only hydrogen is mobile to scan the grain 
surface (e.g.. Id'Hendecourt et al J IT985t lHasegawa & H erbst 
[T99HlKatzetalJI 19991) . Ultimately, this process leads to for- 
mation of saturated ices, e.g. CH4 and C2H6, which have 5% 
chance to escap e directly to the gas upon surface recombina- 
tion (Sect. 12.2b . Our model has a restricted set of hydrocar- 
bon chemistry involving species with < 10 carbon atoms and 
mostly having up to 4 hydrogen atoms only, and thus these 
processes do not appear in Table|9] In the inner disk midplane 
(r < 20 AU) warmed up by the accretion to the temperatures 
of ^40-50 K, desorption of CH4 dominates over its sticking 
to grains, forming an oasis of abundant gas-phase methane. 
At these elevated temperatures hea vier radicals like C, O. 
OH, etc. become partly mobile (e.g..lGarrod & Herbstll2006l 
iGarrod et all 120081; iHerbst & van Dishoeckll2009l) . However, 
the surface growth of hydrocarbons by addition of atomic car- 
bon is still not as efficient as the surface hydrogenation and 
the gas-phase chemistry, and accounts for only a few percent 
of the total surface production rate. The most favorable condi- 



tions for active surface hydrocarbon chemistry are met in the 
warm molecular layer, where atomic hydrogen is plentiful in 
the gas, and its accretion rate is fast even in comparison with 
rapid desorption. More importantly, complex ices are slowly 
destroyed by X-ray/CRP-driven photons as well as by partly 
absorbed stellar (inner disk) and interstellar UV (outer disk) 
photons, producing various reactive radicals. For example, 
CsH is produced by photodestmction of CsH2 ice, followed 
by evaporation (Table [9). In the outer disk partly transparent 
to the IS UV radiation these surface photoprocesses become 
important even in the midplane (see C2H2 in Fig. [8}. 

Note that our chemical network lacks the surface forma- 
tion of benzene apart from its accretion from the gas. The 
major production pathway for C6H6 is via ion-molecule re- 
actions of C3H4 and C3H4 (a = 7.5 10~ 10 cm 3 s" 1 ) as well as 
C 2 H 4 and C 6 H+ (a = 5.5 10" 11 cm 3 s" 1 ), followed by the dis- 
sociative recombination, and freeze-out at T < 150 K. In turn, 
C3HJ and CgHj are synthesized from neutral hydrocarbons by 
rapid charge transfer with ionized hydrogen and carbon or via 
protonation by H3, H30 + , HCO + , etc. The benzene building 
blocks, C3H3 and C3H4, are produced effectively via a num- 
ber gas-phase and surface reactions as described above. The 
characteristic chemical timescale for benzene is thus mainly 
determined by the evolution of C3H3 and C3H4, which in- 
volves slow surface routes. 

As we already discussed in Sect. 13.11 the surface chemistry 
timescales are the longest in the disk chemistry, typically ex- 
ceeding several Myr (see also Tables [3] and 0). Photoprocess- 
ing of ices on dust grains has similar timescale in the molec- 
ular layer and the disk midplane, r c h em ^ 10 6 years (Fig. [3}. 
The overall pace of the surface chemistry is also dependent 
on the amount of involved surface processes. Thus, the char- 
acteristic chemical timescales for hydrocarbons are beyond 
1 Myr, particularly for the heaviest and most saturated ones 
in the network (e.g., C8H4, etc.). Not surprisingly, the 2D- 
turbulent mixing strongly affects the hydrocarbon chemistry 
in disks. The evolution of heaviest carbon chains with more 
than 6 carbon atoms is so far from a steady-state that it be- 
comes sensitive to transport processes even in the slow mixing 
model (see, e.g., the C$H2 ice in Fig.[8j. The vertical transport 
brings up dust grains from the midplane to the upper, more 
irradiated and warmer disk regions, allowing more efficient 
photoprocessing of ices, photodesorption of frozen hydrocar- 
bons, and their photodissociation. These processes increase 
relative abundances of gas-phase hydrocarbons by up to sev- 
eral orders of magnitude. In turn, vertical transport down- 
ward allows to retain newly formed hydrocarbon radicals in 
the icy mantles and facilitates their slow hydrogenation. The 
radial mixing does not play a major role in the hydrocarbon 
chemistry. Consequently, abundances of the most saturated 
hydrocarbon ices in the network are increased by vertical mix- 
ing, in particular in the inner disk with the shortest dynamical 
timescales (e.g., benzene ice), whereas the intermediate prod- 
ucts become less abundant (e.g., C$H2 ice). Note also that the 
mixing allows CH2 and C2H to be more efficiently formed in 
the outer disk atmosphere, where their second layers of high 
relative abundances are developed. This is caused by the tur- 
bulent transport upward of the molecular hydrogen from the 
molecular layer, which is otherwise slowly destroyed in the 
atmosphere by the X-rays and the cosmic ray particles. 

4.3. Oxygen-containing molecules 
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R, AU R, AU 

Figure 9. The same as in Fig.|7]but for the O-containing species. Results are shown for CO, CO ice, CO2, CO2 ice, H2O, O, O2, and H2CO. 

Table 10 

Key chemical processes: O-bearing species 



Reaction 


a 





7 


vmin 


tmax 




[(cm 3 )s-'] 




[K] 


[yr] 


[yr] 


CO2 ice + hi^cRP — > CO ice + O ice 


1.71(3) 








1.00 


5.00(6) 


CO + UV -> + C 


0.20 (-9) 





3.53 


1.00 


5.00(6) 


C0 2 + UV -> CO + 


0.89 (-9) 





3.00 


1.00 


5.00(6) 


H 2 + UV OH + H 


0.80 (-9) 





2.20 


1.00 


5.00(6) 


O + grain — > O ice 


1.00 








1.00 


5.00(6) 


O2 + grain — ► O2 ice 


1.00 








1.00 


5.00(6) 


CO + grain -* CO ice 


1.00 








1.00 


5.00(6) 


CO2 + grain — » CO2 ice 


1.00 








1.00 


5.00(6) 


H 2 CO + grain -► H 2 CO ice 


1.00 








1.00 


5.00(6) 


H2O + grain — >■ H2O ice 


1.00 








1.00 


5.00(6) 


O ice -> O 


1.00 





8.00(2) 


1.00 


5.00(6) 


2 ice — > O7 


1.00 





1.00(3) 


1.00 


5.00(6) 


CO ice -> CO 


1.00 





1.15(3) 


1.00 


5.00(6) 


C0 2 ice -> CO2 


1.00 





2.58(3) 


1.00 


5.00(6) 


H 2 CO ice -> H 2 CO 


1.00 





2.05(3) 


1.00 


5.00(6) 


H ice + OH ice -> H^O 


1.00 








1.00 


5.00(6) 


H ice + HCO ice -> H 2 CO 


1.00 








1.00 


5.00(6) 


OH ice + CO ice — > CO2 ice + H ice 


1.00 





80.00 


1.00 


5.00(6) 


OH ice + CO ice -> C0 2 + H 


1.00 





80.00 


1.00 


5.00(6) 


CH 2 ice + O2 ice -> H,CO + 
CO + + H -> CO + H+ 


1.00 








1.96 


5.00(6) 


0.40 (-9) 








1.00 


5.00(6) 


C0 2 + H+ -> HCO+ + O 


0.30 (-8) 








1.00 


5.00(6) 


H 2 + H+ -> H 2 + + H 
O + H+ -> + + H 


0.73 (-8) 


-0.50 





1.00 


5.00(6) 


0.70 (-9) 





2.32(2) 


1.00 


5.00(6) 


O2 + H+ -> 2 + + H 


0.12 (-8) 








3.82 


5.00(6) 


HC0 2 + + CO -> C0 2 + HCO+ 


0.25 (-9) 


-0.50 





1.00 


5.00(6) 


0+ + H -> O + H+ 


0.70 (-9) 








1.00 


5.00(6) 


O- + CO -> C0 2 + e" 


0.65 (-9) 








7.48 


5.00(6) 


CO + He+ ->■ O + C+ + He 


0.1 6 (-8) 








1.00 


5.00(6) 


H 2 CO + H 3 + -> H^CO+ + H 2 


0.55 (-8) 


-0.50 





1.96 


5.00(6) 


H + OH -> H^O 


0.40(-17) 


-2.00 





7.48 


5.00(6) 


CO + OH -» C0 2 + H 


0.28(-12) 





1.76(2) 


1.00 


5.00(6) 


+ C 2 -> CO + C 


1.00 (-10) 








1.00 


6.11(3) 


+ CH 3 -» H 2 CO + H 


0.14 (-9) 








1.00 


5.00(6) 


+ HCO -> C0 2 + H 


0.50(-10) 








1.00 


5.00(6) 


O + OH -> 2 + H 


0.75 (-10) 


-0.25 





1.00 


5.00(6) 


O + H,CO -> CO + OH + H 


1.00 (-10) 








1.00 


5.00(6) 


o 2 + C 3 -> C0 2 + c 2 


1.00(-12) 








1.00 


5.00(6) 


o 2 + c -> O + CO 


0.47 (-10) 


-0.34 





1.00 


5.00(6) 
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Table 10 

Key chemical processes: O-bearing species 



H 3 + + e" -> H 2 + H 
HC0 2 + + e" -> C0 2 + H 
H 2 CO + + e" ->■ H^CO 
H 3 + + e" -> O + H + H, 
2 + + e" -y O + O 

H 1 + + grain(-) -+ H 2 + H + grain(O) 
CH 3 2 + + grain(-) -» C0 2 + H 2 + H + grain(O) 
2 + + grain(-) — s- O + O + grain(O) 



0.11 (-6) 


-0.50 





1.00 


5.00(6) 


0.60 (-7) 


-0.64 





1.00 


5.00(6) 


0.11 (-9) 


-0.70 





14.60 


5.00(6) 


0.56 (-8) 


-0.50 





3.82 


5.00(6) 


0.1 9 (-6) 


-0.70 





1.96 


1.75(5) 


0.25 








1.00 


5.00(6) 


0.50 








28.60 


3.42(5) 


1.00 








1.96 


1.75(5) 



The chemical evolution of O, O2, CO, CO2 and water is 
of particular attention in astrophysics since the recent puz- 
zling results obtain ed with Spitzer, Herschel, SWAS, and 
Odin satellite s (e.gjLarsson et al.ll2007t ICarr & Naiitall2008l 
Ber gin et alJ [2010; J0r gensen & van Dishoeckl 120101). and 
their potential relevance for astrobiology (e.g.. ISelsis et al] 
l200l iKaltenegger et all 12001 ISegura etafl 120071) . In this 
subsection we analyze in detail chemical and mixing pro- 
cesses responsible for the evolution of oxygen-bearing species 
in protoplanetary disks. We focus on simple neutral molecules 
that are composed of O, C, H only, and consider complex 01- 
ganics in separate subsection. 

There are 4 steadfast neutral O-bearing molecules (CO, 
OH, H2CO, and the water ice), and 3 steadfast ions (H2CO + , 
H3CO + , and H30 + ) (Table |5j. Their column densities in the 
laminar model and the fast mixing model differ by a factor 
of 3. Among the species sensitive to the turbulent transport 
(Table 13 there are 5 molecules (C3O, H 2 0, HCO, O, 2 ), 3 
ions (HCO + , + , OJ), and 4 solid species (CO, H 2 CO, O, OH 
ices). Their column densities are changed by up to 2 orders of 
magnitude by the turbulent transport. The hypersensitive O- 
bearing species (Table |7]i include 3 molecules (C02, H2O2, 
O3), 2 ions (OH + , H 2 + ) and 5 ices (C0 2 , H 2 2 , 2 , 2 H, 
O3). The turbulent diffusion alters their column densities by 
up to a factor of 4000 (CO2 ice). Apparently, neutral species 
with a large number of oxygen atoms (e.g., ozone) affected 
more strongly by the turbulent mixing than the molecules 
containing a single O (e.g., CO). Similar to the hydrocarbon 
chemistry, abundances and column densities of O-bearing ices 
are stronger altered by the turbulent diffusion compared to the 
gas-phase molecules. Once again this is an indicator that the 
chemical evolution of the multi-oxygen molecules in proto- 
planetary disks is at least partly governed by slow chemical 
processes (surface recombination, photodissociation of ices, 
etc.). 

In Fig. [9] the distributions of the relative molecular abun- 
dances and column densities at 5 Myr of CO, CO ice, CO2, 
CO2 ice, H2O, O, O2, and H2CO calculated with the laminar 
and the 2D-mixing models are presented. The abundance dis- 
tribution of atomic oxygen shows a 2-layered structure, with 
a maximum of relative concentrations (X(0) w 1.5 10~ 4 ) in 
the disk atmosphere and an upper part of the molecular layer, 
and strong depletion in the midplane (X(0) < 10 -11 ). On the 
other hand, CO and CO2 ices are concentrated in the mid- 
plane and a lower part of the molecular layer (at z w 1 H r ), 
with typical relative abundances of < 10~ 5 . The CO2 ice 
abundance is particularly high in the inner warm midplane 
(r < 20 AU). Its distribution is also 2-layered. The rest of 
considered O-molecules show a 3-layered abundance struc- 
ture similar to that of the polyatomic ions and hydrocarbons. 



The carbon monoxide spreads over a wide vertical h eights 
in the disk due to se lf-shielding (e.g.. I van Dishoeck 1988; 
iLyons &"Yo ung 2005l lVisser et al.l2009al) a nd low desorption 
energy of - 1 000 K (e.g.. lBisschop et alJl200d) . z « 1 - 2.5 
at r = 10 AU and z « 0.8-2 at r = 800 AU, with a typical 
abundance of 7.6 10~ 5 . The tip of low CO abundances in 
the inner low atmosphere at 2H r is due to enha nced X-ray- 
driven dissociation in this region (see also Sect. 13. It . Note 
that gas-phase CO in the laminar model is partially produced 
in the very inner disk midplane (r < 20 AU). The H2O is also 
wide-spread throughout the disk, though due to photodissoci- 
ation and rapid freeze-out at T < 120 K its peak abundances 
of 10~ 8 - 10~ 7 are confined to the bottom of the molecular 
layer, with a maximum in the outer disk region at r > 200 AU 
(Pig. O. Gas-phase water comprises less than 0.01% of the 
total water abundance that are locked in dust mantles in the 
midplane. The O2 abundance distribution has almost the same 
pattern as that of gas-phase water, with maximum values of 
X(0 2 ) f» 10" 7 -310" 5 . The peak C0 2 relative abundances 
are narrow and rather low, ~ 10~ 10 - 10~ 8 , and restricted to 
the very outer disk (r > 400 AU). In contrast, the H2CO has 
a radially-uniform, narrow (w 0.3H r ) molecular layer, with 

maximum abundances of ~ io~ 10 ~ 10 at z « 0.8 - 1 H r . 

Turbulent mixing does not considerably affect column den- 
sities of H2CO and CO in gas and solid phases, although it 
populates the viscously-heated inner midplane with carbon 
monoxide, and slightly widens their molecular layers upward 
(Fig. [9]). Unlike atomic ions, the mixing softens the vertical 
gradient of the atomic oxygen abundances, and lowers its col- 
umn density. Turbulent diffusion transports solid CO 2 from 
the inner midplane radially and vertically outward, enhanc- 
ing its column densities and abundances by more than 3 or- 
ders of magnitude at r > 30 AU. Finally, column densities of 
gas-phase CO2, H2O, and O2 are increased by the turbulence 
by factors of 10-1 000, with their molecular layers enriched 
and vertically expanded up to ~ 1 . 2 - 2 H r and more homoge- 
neously distributed in the radial direction. 

To better understand these results, we performed detailed 
analysis of the chemical evolution of CO, CO ice, CO2, CO2 
ice, H2O, O, O2, and H2CO in two disk vertical slices at r = 10 
and 250 AU (the laminar chemical model). The most im- 
portant reactions responsible for the evolution of their abun- 
dances in the midplane, the molecular layer, and the atmo- 
sphere are presented in Table [10] both for the inner and outer 
disk regions. The final list contains only top 25 reactions per 
region (midplane, molecular layer, atmosphere) for the entire 
5 Myr time span, with all repetitions removed. 

The chemical evolution of the selected O-bearing species 
is governed by a limited set of reactions. Atomic oxygen is 
present in the disk atmosphere and converted to CO, CO2 
and H2O in the midplane and the molecular layer. The rate 
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of its convertion in the molecular layer is partly regulated by 
the evolution of H3 and hydrocarbons, which is i n tur n deter- 
mined by the slow X-ray irradiation of H2 (Sect. 14.11 ) and the 
slow release of O from CO b y the X-ray-ionized helium atoms 
in the inner disk (Sect. 14. 2j , and the surface chemistry of O- 
bearing species. The resulting chemical timescale of > 1 Myr 
exceeds the transport timescale in the molecular layer (Ta- 
bles H and 0. 

The water ice forms on the dust surfaces via accretion of 
the gas-phase water in disk regions with T < 120 K, and via 
surface hydrogenation of frozen atomic oxygen and hydrox- 
yle (minor route). The destruction pathways for water ice 
are thermal or photoevaporation (major route), and further 
hydrogenation to hydrogen peroxide ice (very minor route). 
In the gas, water formation begins by production of OH + 
from O and H3, followed by subsequent hydrogen abstrac- 
tion reactions with H2 till H30 + is created. The protonated 
water dissociatively recombines with electrons or negatively 
charged grains into H 2 (25%) or OH (74%) or O (1%), or 
de-protonates by ion-molecule reactions with other abundant 
neutral molecules (CH4, CO, etc.). The gas-phase formation 
of water is assisted by slow radiative recombination reaction 
between H and OH, and photoevaporation of water ice in the 
warm molecular layer. The gas-phase removal channels for 
H2O include photodissociation, charge transfer reactions with 
H + followed by dissociative recombination, freeze-out, and 
reactions with He + (minor route). The solid water is a ter- 
minal species that serves as one of the sinks of the elemental 
oxygen in the disks. Thus, the key chemical processes lead- 
ing to the evolution of water are fast, r C h e m ^ 1 - 10 2 years, 
in the atmosphere and the midpane, whereas in the molecular 
layer it is regulated by the late-time evolution of H3 due to the 
X-ray ionization (r c hem ^ 10 4 years; see Sect. 14. 11 1. 

The molecular oxygen is produced in the gas by neutral- 
neutral exothermic reactions of OH and O, and destroyed in 
ion-molecule and neutral-neutral combustion reactions with 
various radicals (mostly dehydrogenated hydrocarbons, e.g., 
C, C3, etc.), by photodissociation, and via ion-molecule reac- 
tions of O2 with H3 and ionized C and H, followed by disso- 
ciative recombination of 2 into atomic oxygen. On dust sur- 
faces molecular oxygen is produced either directly from the 
recombination of oxygen atoms at conditions when surface O 
becomes mobile (T > 30 K; the inner disk midplane and the 
warm molecular layer), or via surface oxidation reactions: O 
(ice) + OH (ice) -> 2 H (ice) and O (ice) + 2 H (ice) 2 
(ice) + OH (ice). The surface O2 can be further converted to 
O3 ice. The characteristic chemical timescale for O2 is set by 
the slow surface chemistry timescale (r c h em ^ 10 6 years). 

CO molecules serve as a sink of almost all elemental car- 
bon and about a half of elemental oxygen in disks. Carbon 
monoxide is formed essentially in the gas-phase via reactions 
of atomic oxygen with CH, CH2, and C2, and other abundant 
hydrocarbons (see previous subsection and Table [9). The re- 
moval pathway of CO in the gas include freeze-out onto the 
grain surfaces (major channel), and the slow ion-molecule 
reaction with He + (minor channel). In the upper, heavily 
UV-irradiated and dilute disk atmosphere CO is also UV- 
photodissociated (z > 2.5//,). At T < 30 K, in the outer mid- 
plane CO sticks to the grains and partly converted to H2CO 
and CH3OH ices via hydrogenation reactions, whereas in the 
warm midplane (r < 30 AU, T> 30 - 40 K) CO is trans- 
formed into CO2 in a slightly endothermic reaction between 
the CO and OH ices (with a barrier of 80 K). The chemistry 



of CO gas reaches a steady-state in the molecular layer and 
the cold outer midplane within less than 10 2 — 10 4 years (see 
Fig.|2]i. The quasi-equilibrium CO chemistry in the molecular 
layer is restricted to protonation of CO molecules by H3 into 
HCO + , balanced by dissociative recombination. In the disk 
atmosphere CO chemistry is controlled by the photodissoci- 
ation slowed down by self- and mutual-shielding by H2, and 
Tchem ^ 10 4 - 10 6 years. Finally, in the inner warm midplane 
the surface conversion of CO into CO2 leads to a very long 
timescale of < 1 Myr. 

The carbon dioxide is a daughter molecule of CO, and is 
mainly produced in the gas phase via oxidation of HCO (O 
+ HCO — > CO2 + H), slow combustion of C3 and C2H, slow 
endothermic reaction of CO and OH with a barrier of 176 K 
(Table [10), and desorption of C0 2 ice at T > 60 K or UV- 
photodesorption. The main removal gas-phase pathways in- 
clude photodissociation, ion-molecule reactions with C + and 
H + (forming CO + and HCO + , respectively), and accretion to 
dust grains. The solid CO2 is mainly produced via endother- 
mic reaction of surface CO and OH in the disk inner midplane 
and the low part of the entire molecular layer, and through ac- 
cretion of the gas-phase carbon dioxide. The CO2 ice is de- 
stroyed by the X-ray/CRP-induced UV photons in the inner 
midplane, and via thermal and UV-desorption. Thus, the CO2 
chemistry has particularly long timescale associated with sur- 
face reaction of CO and OH, and slow photoprocessing of the 
C0 2 ice (r chem > 1 Myr). 

Finally, the H2CO molecule is mostly produced in the gas- 
phase through reaction of CH3 with oxygen atoms, and des- 
orption of formaldehyde ice. The major gas-phase removal 
routes are sticking to dust in the disk regions with T < 
40-50 K, photodissociation, reactions with ionized atomic 
C and H, oxidation by O (into CO, OH, and H), and proto- 
nation by H3 followed by dissociative recombination back to 
H 2 CO (33%), or CO (33%), or HCO (33%). The surface evo- 
lution of H2CO is governed by accretion and desorption pro- 
cesses, and a sequence of surface hydrogenation of CO into 
CH3OH where formaldehyde ice is an intermediate product. 
The timescale of key evolutionary processes for formalde- 
hyde, namely, oxidation of CH3 and H2CO, as well as accre- 
tion and evaporation, is fast, < 10 2 - 10 3 years (Tables[3}j4]l. 

As a result, the column densities of formaldehyde are only 
slightly affected by the turbulent transport (see Fig. |9}- The 
same holds true for very abundant gas-phase CO and the 
H2O ice. Their global chemical evolution is only slightly 
controlled by the surface chemistry, and the CRP/X-ray- 
irradiation becomes important only in the upper regions that 
do not contribute to the resulting vertical column densi- 
ties. The 2D-mixing enables more efficient production of 
formaldehyde in the atmosphere (at w l-2H r ) thanks to en- 
hanced abundances of CH3 (see Sect. |4.2| i. The turbulent dif- 
fusion lowers the column densities of atomic oxygen, and in- 
creases the column densities of gas-phase water, carbon diox- 
ide in all phases, and molecular oxygen. The pace of the 
conversion of atomic oxygen into other O-bearing species is 
partly governed by HI, which is sensitive to transport (see 
discussion in Sect. 14. 11 1. The gas-phase water production in- 
volves this slow process and slow desorption, and thus gas- 
phase water becomes sensitive to the mixing. Finally, CO2 
and O2 (and other multi-0 species in the model) are the 
molecules which chemical evolution is at least partly deter- 
mined by the slow surface reactions. This makes them sen- 
sitive to the turbulent mixing. Note also that the chemical 
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evolution of the CO2 ice and, to a less degree, the CO ice, is 4.4. Nitrogen-containing molecules 

influenced strongly by the radial mixing. 

Table 11 

Key chemical processes: N-bearing species 



Reaction 




P 


7 


tm in 


tlTMX 




[(cm 3 )s-'] 




[K] 


[yr] 


[yr] 


HC5N ice + hi/ CR p — > C4H ice + CN ice 


1.75(3) 








1.00 


5.00(6 


HCN + UV -> CN + H 


0.1 6 (-8) 





2.69 


1.00 


5.00(6 


HNC + UV — >• CN + H 


0.55 (-9) 





2.00 


1.00 


5.00(6 


HNO + UV — s> NO + H 


0.1 7 (-9) 





0.53 


1.00 


5.00(6 


OCN + UV — > + CN 


l.OO(-ll) 





2.00 


1.00 


5.00(6 


CN + grain — ► CN ice 


1.00 








1.00 


5.00(6 


HCN + grain — ► HCN ice 


1.00 








1.00 


5.00(6 


HNC + grain -> HNC ice 


1.00 








1.00 


5.00(6 


HNO + grain — ¥ HNO ice 


1.00 








1.00 


5.00(6 


NO + grain — > NO ice 


1.00 








1.00 


5.00(6 


OCN + grain — > OCN ice 


1.00 








1.00 


5.00(6 


HC5N + grain — ¥ HC5N ice 


1.00 








1.00 


5.00(6 


HCN ice — > HCN 


1.00 


2 


05 (3) 


1.00 


5.00(6 


HJNC ice — > HJNC 


1.00 


2 


05 (3) 


1.00 


5.00(6 


T T/~" XT : , T T/"^ XT 

HC5N ice — y HC5JN 


1.00 


6 


18(3) 


1.00 


5.00(6 


H ice 4- CN ice — > HCN 


1 m 
1.00 








1 .00 


5.00 (6 


H ice + C5N ice — ■> HC5N ice 


1 m 
1.00 








1 .00 


5.00 (6 


TT - 1 A 1 XT I ! T T/~" XT 

H ice + C5JN ice — > HCgN 


1.00 








28.60 


5.00(6 


TT ■ . TT/" 1 XT i ! TT f ^ XT : 

H ice + HC5JN ice — > H2C5N ice 


1.00 


1 


21(3) 


1.00 


5.00(6 


H ice + NO ice — » HNU 


1.00 








1.00 


5.00(6 


xt 1 1 f tt : 1 t t/~* xt : 

JN ice + C5H ice — > HC5N ice 


1.00 








1.00 


5.00(6 


U ice + CN ice — > OCN 


1.00 








4.18(2) 


5.00(6 


/" "\ • , TTXTA I ! \|/"\ , ATT 

U ice + HNO ice — >► JNU + OH 


1.00 








1.00 


5.00(6 


/~1T T /— 1 XT . T T-l- 1 T 1 /" 1 XT 1 /~iT T 4- 

CH3C5N + H + — *■ HC5N + LHj, 


0.20 (-7) 


-0.50 





3.82 


5.00(6 


HC5N + H + — > CsHN + + H 


0.40 (-7) 


-0.50 





1.00 


5.00(6 


HCN + H + — > HCN + + H 


0.28 (-7) 


-0.50 





1.00 


5.00(6 


HNC + H + — > HCN + H + 


0.25 (-7) 


-0.50 





1.00 


5.00(6 


HCN + + H — > HCN + H + 


0.37 (-10) 








1.00 


5.00(6 


XTn 1 /"I i /"'XT . / "\ 

NO + C — > CN + O 


0.60 (-10) 


-0.16 





1.00 


5.00(6 


OCN + C — > CN + CO 


1.00 (-10) 








1.00 


5.00(6 


/*1XT ■ /"A , AHXT 1 / \ 

CJN + O2 — > UCN + O 


0.14 (—10) 


—0.60 





1.00 


5.00(6 


[_| 7\[ f ~» . ( \ , XTfl . / \ 1_J 

HJNU + U — t JNU + UH 


A TO/ 1 rt\ 








1 .UU 


5.00 (6 


XTT T 1 f~\ 1 TTXTjTA 1 TT 

JNH2 + U— > HNO + H 


r\ OA/ 1 Ai 

u.8U(— 1U) 








1.00 


5.00(6 


C 7 N + 0->-OCN + C6 


U.4U (— 1U) 


U 


U 


1 An 
1 .UU 


j.UU (0 


/"•M -i_ (Ml 1 ( \j \ . it 

+ Un — r UL.1N + ri 


U. / VJ( — L\J} 




u 


1 nfi 

1 .UU 


J.UU \o 


I 1 \ | ) , IT v \ ( ) , IT 


n /is/ 1 m 






1 on 

1 .UU 


J.UU ^0 


N + CH —7- CN + H 




— 0.09 





1.00 


2.56(6 


N + CH 2 -> HNC + H 


0.40 (-10) 


0.17 





1.00 


5.00(6 


N + HCO -> OCN + H 


1.00 (-10) 








1.00 


5.00(6 


N + NO N 2 + O 


0.30 (-10) 


-0.60 





1.00 


5.00(6 


N + OH -> NO + H 


0.75 (-10) 


-0.18 





1.00 


5.00(6 


C 3 H,N+ + e~ -> HNC + C 2 H 


0.75 (-7) 


-0.50 





1.00 


5.00(6 


dHN+ + e - -> HCN + C 2 


0.30 (-6) 


-0.50 





1.00 


5.00(6 


C 5 H 2 N+ + e- -> HC 5 N + H 


0.15 (-6) 


-0.50 





1.00 


5.00(6 


C 5 H 3 N+ + e" -r HC 5 N + Ht 


1.00 (-6) 


-0.30 





7.48 


5.00(6 


C 6 H 4 N+ + e" -> HC 5 N + CH 3 


1.00 (-6) 


-0.30 





1.00 


5.00(6 


CH 2 CN+ + e" ->■ HCN + CH 


0.30 (-6) 


-0.50 





1.00 


5.00(6 


H 2 CN + + e" -> HCN + H 
H 2 CN+ + e" -> HNC + H 


0.1 9 (-6) 


-0.65 





1.00 


5.00(6 


0.1 9 (-6) 


-0.65 





1.00 


5.00(6 


H 2 NO + + e" -> HNO + H 


0.15 (-6) 


-0.50 





1.00 


5.00(6 


H 2 CN+ + e" ->• CN + H + H 


0.92 (-7) 


-0.65 





1.00 


5.00(6 



Similar to the oxygen-bearing molecules, the observa- 
tions of the N-bearing species provide useful diagnostics 
on the physical conditions and chemical state of the cos- 
mic objects. The ammonia do ublet lines are used to 
measure kinetic temperatures ( e.g.. IWalmslev & Ungerechtsl 
[1981 iChurchwell etaH IT990t Uiiina etal.lll999ir while the 
CN and HCN relative line st rengths are sensitive to the UV 
flux (e.g.. Be rgin et al.ll2003l) . The linear HC3N and heavier 
cyanopolyy nes with large dipo le moments serve as densito- 
meters (e.g. jPratap et al.l ll997). whereas OCN" has been pro- 
posed as a carrier of the 4.62/im absorption feature in inter- 



stellar ices dRaunier et al.l 120031: iBennett et aill2010l) . In this 
Section we analyze in detail chemical and mixing processes 
responsible for the evolution of the neutral nitrogen-bearing 
species in protoplanetary disks. 

There are 4 steadfast neutral N-bearing molecules (CN, 
HCN, HNC, NH 2 ), 2 ions (H 2 CN + and NHJ) (Tabled). Their 
column densities in the laminar and fast mixing models dif- 
fer by a factor of 3. The sensitive nitrogen species include 
19 molecules (e.g., C 2 N, C 3 N,..., C 9 N, HC 3 N), 6 ions (e.g., 
NO + , N 2 H + , CNC + ), and 22 solid species (e.g, CH 2 NH, 
HC 3 N,.., HC 9 N, CN, NO; Table H). Their column densi- 
ties are changed by up to 2 orders of magnitude by the tur- 
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Figure 10. The same as in Fig.[7jt>ut for the N-containing species. Results are > 

bulent transport. The N-bearing hypersensitive species are 
4 molecules (C 7 N, HC 7 N, N 2 0, and N0 2 ), N + , and 8 ices 
(e.g., CH3C3N, HNO, N 2 , NO2, etc.), see Table [7] The tur- 
bulent diffusion alters their column densities by up to 7 or- 
ders of magnitude (H5C3N ice). Similarly to hydrocarbons 
and oxygen-containing molecules, complex chains with mul- 
tiple N, C, or H atoms are more strongly affected by the 
turbulent transport than chemically simpler species, though 
the trend is not t hat clear (e.g., NO?, H NCO, N + , N 2 are 
the outliers). As IVasyunin et al.l (120081) have found, nitro- 
gen chemistry in disks involves a larger number of key re- 
actions, including many exothermic neutral-neutral reactions, 
surface processes, compared to the chemistries of the O- and 
C-containing species. We have selected several most interest- 
ing nitrogen species for detailed chemical analysis. 

In Fig. [TO] the distributions of the relative molecular abun- 
dances and column densities at 5 Myr of CN, HCN, HNC, 
NO, HNO, OCN, HC 5 N, and HC 5 N ice calculated with the 
laminar and the 2D-mixing models are presented. The relative 
abundance distributions of the all considered species show a 
3-layered structure, with peak concentrations in the molecu- 
lar layers at ps 0.8- lH r , with typical values of 10~ 1() - 10~ 8 
(with respect to the total amount of hydrogen nuclei). The 
molecular layers of N-bearing molecules are narrow, « 0.2- 
0.5 H r , similar to those of hydrocarbons, HCO + and H 2 CO 
(see Figs. [8] and |9). The photostable CN radical has a sec- 
ond molecular layer in the outer disk atmosphere, at z ~ 1.8 
(r > 200 AU), though it does not contribute much to the total 
CN column density. Note that HC5N ice is also concentrated 
in the molecular layer, at ~ 1 pressure scale height. 

Turbulent mixing does not affect column densities of stead- 
fast CN, HCN, and HNC species (see 4th panels in Fig.flOb. 
Nonetheless, their molecular layers are broadened by diffu- 
sion, and the second, upper molecular layer of CN becomes 
more prominent (compare 1st and 3rd panels in the Figure). 



5 1 — 1 — f — 1 — f — 1 — ai — 1 — u 16.0 




R, AU 

vn for CN, HCN, HNC, NO, HNO, OCN, HC5N, and HC 5 N ice. 

NO, HNO, OCN, HC 5 N, and HC 5 N ices are sensitive to the 
mixing, with their column densities increased by 2D-turbulent 
diffusion by up to a factor of 40 (cyanodiacetylene; Table |6}. 
The corresponding abundance distributions are vertically ex- 
tended up to the thicknesses of ~ 0.5- 1.5//,, and enhanced 
by the transport by up to several orders of magnitude. 

To better understand these result, we investigate the evolu- 
tion of CN, HCN, HNC, NO, HNO, OCN, HC 5 N, and the sur- 
face HC5N in the two disk vertical slices at r = 10 and 250 AU 
(the laminar chemical model). The most important reactions 
responsible for the time-dependent evolution of their abun- 
dances in the midplane, the molecular layer, and the atmo- 
sphere are presented in Table Q~T] both for the inner and outer 
disk regions. The final list contains only top 20 reactions per 
region (midplane, molecular layer, atmosphere) for the entire 
5 Myr time span, with all repetitions removed. 

The evolution of cyanopolyynes is tightly conn ected with 
the evolution of carbon chains discussed in Section l4~2l Their 
major production pathways in the gas include evaporation 
from icy mantles, neutral-neutral reaction of N with Renner- 
Teller hydrocarbons (C„H), slow or slightly endothermic re- 
action of CN with C„_iH 2 (e.g., CN + C 4 H 2 ->■ HC 5 N), and 
at later times dissociative recombination of their protonated 
analogs formed by the ion-molecule reactions with H3 (Ta- 
bleCnji. The main removal pathways are the freeze-out at tem- 
peratures < 70- 170 K, UV-photodissociation, reactive colli- 
sions with the ionized C, H, and He atoms. On dust surfaces, 
cyanopolyyne ices form either via addition of N to the frozen 
Renner- Teller hydrocarbons or through hydrogenation of the 
C„N ices. The major destruction routes for the cyanopolyyne 
ices are photoevaporation, photoprocessing by the X-ray- or 
CRP-driven UV photons (e.g., HC 5 N ice + h^ CRP -> C 4 H 
ice + CN ice), and surface conversion to even more complex 
species (e.g., HC5N ice + H ice — > H 2 CsN ice). In turn, de- 
ny drogenated carbon chains with attached nitrogen atom are 
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produced in the gas and on the dust surfaces by rapid neutral- 
neutral reactions of N with the C„H species, and dissocia- 
tion of complex molecules with multiple C and N atoms (e.g., 
CH3C5N ice + hv CRP — !■ C 5 N ice + CH 3 ice). Obviously, as in 
the case of complex carbon chains, characteristic timescales 
of the cyanopolyyne chemical evolution are regulated by the 
slow surface processes with r c h em exceeding million years. 

In contrast, chemical histories of CN, HCN, and HNC are 
closely related and governed by a small set of reactions. The 
production and destruction of CN proceeds entirely in the 
gas phase. The primal formation pathways are photodisso- 
ciation of hydrogen cyanide and isocyanide, rapid barrierless 
neutral-neutral reactions of N with CH, C2, and C2N, and 
NO with C (rate coefficients are ~ 10 _11 -610" 10 cm 3 s _1 ), 
and dissociative recombination of HCNH + at later times (t > 
10 2 — 10 4 years), see Table [TT] Another, less important forma- 
tion pathway for CN is the neutral-neutral reaction between 
C and OCN leading to CN and CO. The destruction of CN is 
mostly caused by photodissociation, fast neutral-neutral reac- 
tions with N, O, OH, and O2, and freeze-out in the disk mid- 
plane at T < 35-40 K. In the atmosphere CN reacts with H + , 
forming CN + , which is converted back to CN by the charge 
transfer with atomic hydrogen. The fact that the CN chem- 
istry involves only gas-phase routes implies its relatively short 
characteristic timescale, ~ 10 3 years (Tables[3]-|4]l. Only in the 
outer atmosphere r c h em (CN) > 10 5 years, because the CN evo- 
lution there depends o n the slow evolution of C2H (see Fig. [8] 
and discussion in Sect. I4.21 >. 

The key production route for gas-phase HCN and HNC are 
the neutral-neutral reaction of nitrogen atoms with CH2 and 
dissociative recombination of HCNH + (Table ITTb. More en- 
ergetically favorable isomer, HCN, is produced upon reactive 
collisions of ionized hydrogen with hydrogen isocyanide. The 
minor formation channels in the disk midplane are direct sur- 
face recombination of atomic hydrogen and cyanogen radical 
that leads to gas-phase HCN and H (with 5% probability), 
and the charge transfer reaction of H and HCN + . Due to their 
relatively large binding energies, thermal desorption of HCN 
and HNC does not occur until temperatures of T > 40 K are 
reached. The major destruction routes for the gas-phase HCN 
and HNC include accretion onto the dust grains, the charge 
transfer reaction with H + , the ion-molecule reaction with C + 
leading to CNC + or C2N" 1 ", and protonation reaction with H3, 
HCO + , and H 3 + (e.g., HCN + HCO + -> HCNH + + CO). An- 
other important destruction channel for hydrogen (iso)cyanide 
in the upper disk layers at ~ 1.5-2 scale heights is photodis- 
sociation. The characteristic timescales of the HCN and HNC 
evolution are < 10 3 - 10 4 years in the molecular layer. Their 
chemical timescales exceed 10 5 — 10 6 years in the midplane, 
where their evolution is partly controlled by the slow surface 
formation, and in the inner upper molecular layer/low atmo- 
sphere subject to the slow X-ray-driven dissociation of H2 and 
release of oxygen from water and CO. 

The evolution of NO, HNO, and OCN, similarly to that 
of CN, HCN, and HNC, is also governed by a set of rapid 
exothermic neutral-neutral reactions. The major formation 
pathways for nitrogen monoxide comprise reactive collisions 
of atomic nitrogen with hydroxyle, reactive collisions of 
atomic oxygen with nitrogen monohydride, gas-phase de- 
struction of HNO by atomic oxygen, and, in the inner mid- 
plane, surface destruction of HNO by O, CH2 and CH3, with 
products (NO and OH, CH3, or CH4) directly injected into 
the gas. In the disk atmosphere NO is produced by the pho- 



todissociation of HNO. In the inner disk midplane a source of 
the NO gas is thermal desorption of the NO ice. The ma- 
jor removal processes for NO include photodissociation in 
the atmosphere, accretion onto the dust grains in the disk 
regions with T < 30 - 35 K, neutral-neutral reactions with 
atomic carbon and nitrogen (leading to CO and N2), and 
charge transfer reactions with H + and C + . The characteris- 
tic timescale for NO is < 10 3 - 10 4 years in the molecular 
layer, and > 10 5 — 10 6 years in the midplane, where its evolu- 
tion is partly controlled by the surface processes. In the upper 
molecular layer at r ~ 100 AU the NO chemical timescale is 
again long, ~ 1 Myr, since it is related with the evolution of 
atomic oxygen that is slow at those disk heights. 

The chemical evolution of nitric acid (HNO) is very sim- 
ilar to the evolution of NO. The major production terms for 
HNO are neutral-neutral reaction of atomic oxygen with ni- 
trogen dihydride, and, in the warm molecular layer and the 
inner midplane, via surface hydrogenation of NO directly to 
the gas-phase. The evaporation of the solid HNO plays only 
a minor role for production of gas-phase HNO (in the lam- 
inar model). The major destruction terms for NHO include 
photodissociation in the atmosphere, accretion onto the dust 
grains in the disk regions with T < 40-50 K, neutral-neutral 
reactions with atomic hydrogen and oxygen (leading to NO), 
and ion-molecule reaction with H + (forming ionized NO and 
H2). The distribution of the characteristic timescale for HNO 
over the disk is close to that of NO. 

Finally, the chemical evolution of cyanate (OCN) is more 
diverse than that of NO and HNO. The key formation routes 
comprise of neutral-neutral reactions of CN with OH and 
O2, N and HNO, and oxidation reactions with abundant 
cyanopolyynes, e.g. C7N + O — >• OCN + C$. It can also 
be produced by direct surface recombination in O and CN 
ices. Similar to HNO, desorption of the solid OCN is not a 
key formation process for gas-phase cyanate (in the laminar 
model). The key removal routes include photodissociation in 
the atmosphere, sticking to the dust grains in the disk regions 
with T < 50-60 K, neutral-neutral reactions with atomic oxy- 
gen and carbon (leading to NO or CO), slow combustion re- 
action with O2 (producing CO2 and NO), and ion-molecule 
reaction with C + (forming ionized CO and CN). The charac- 
teristic timescale of the OCN evolution is < 10 3 - 10 4 years in 
the molecular layer, and > 10 5 — 10 6 years in the inner and 
outer midplane, where its evolution is partly controlled by 
the surface processes. In the atmosphere the OCN chemical 
timescale increases from 10 4 till 10 6 years, similar to the CO2 
timescale shown in Fig. [2] 

Consequently, column densities of CN, HCN, and HNC are 
not much affected by diffusion as their key evolutionary pro- 
cesses in their molecular layers proceed rapidly in the gas- 
phase (Fig. [Tol l. As for many molecules in the model, the 
turbulent transport expands their molecular layers in vertical 
direction. Since the chemistry of CN involves C2H, which 
abundances are increased by the mixing in the atm osphere 
at z > 2H r (see Fig. [8] and discussion in Sect. 14.21 ). in the 
2D-mixing case (Sc =1) the second molecular layer of CN 
is formed. Similarly, the hydrogen cyanide and isocyanide 
evolution is related to the evolution of CH2, which is sensitive 
to the transport in the same upper disk region, though not as 
strong as ethynyl radical. The chemistry of the NO-containing 
species (e.g., NO, HNO, OCN) is sensitive to mixing as it 
depends on the evolution of O and O2 that is in turn influ- 
enced by the transport, and since their production involves 
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minor slow surface routes. Moreover, thermal evaporation of Finally, due to extreme importance of surface chemistry for 

these heavy species is inefficient, and thus photoevaporation the chemical evolution of cyanopolyynes, their abundances 

of their ices becomes an important production pathways for and column densities are increased by the turbulent mixing, 

the gas-phase counterparts. Vertical diffusive mixing allows albeit not as strongly as for some heavy hydrocarbons, 
more proficient evaporation of NO, HNO, and OCN ices in the 

warm molecular layer, enhancing their gas-phase abundances. 4.5. Sulfur-containing molecules 



Table 12 

Key chemical processes: S-bearing species 



Reaction 
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tmax 
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1 50(3) 








1 00 


5 00(6 
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1 00 


5 00 (6 


TtS + TIV — v Pi + S 
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l .\jyj \ 1 \j j 





2.00 


1.00 


5.00(6 


H2CS + UV — > CS + H2 


[ 00 (-9) 





1.70 


1.00 


5.00(6 


H 2 S + UV — y HS + H 


0.31 (-8) 





2.27 


1.00 


5.00(6 


OCS + UV — y CO + S 


0.37 (-8) 





2.07 


1.00 


5.00(6 


S0 2 + UV — y SO + 


O.i 9 (-g) 





2.38 


1.00 


5.00(6 


CS + grain —y CS ice 


1.00 








1.00 


5.00(6 


C2S + grain — » C2S ice 


1.00 








1.00 


5.00(6 


C 3 S + grain — y C3S ice 


1.00 








1.00 


5.00(6 


OCS + grain — y OCS ice 


1.00 








1.00 


5.00(6 


SO? + grain — y SO2 ice 


1.00 








1.00 


5.00(6 


H2S + grain — y H7S ice 


1.00 








1.00 


5.00(6 


H 2 CS + grain -y H 2 CS ice 


1.00 








1.00 


5.00(6 


CS ice -> CS 


1.00 





1.90(3) 


28.60 


5.00(6 


C 2 S ice C 2 S 


1.00 





2.70(3) 


1.00 


5.00(6 


C 3 S ice ->■ C 3 S 


1.00 





3.50(3) 


1.00 


5.00(6 


C2S ice + C ice — > C 3 S ice 


1.00 








1.00 


5.00(6 


Hice + HS ice->H,S 


1.00 








1.00 


5.00(6 


O ice + CS ice -+ OCS 


1.00 








1.00 


5.00(6 


O ice + SO ice -y S0 2 


1.00 








1.00 


5.00(6 


S ice + CO ice -> OCS 


1.00 








1.00 


5.00(6 


CS + H+ -> CS+ + H 


0.1 8 (-7) 


-0.50 





1.00 


5.00(6 


C 2 S + H+ -!- C 2 S+ + H 


0.11 (-7) 


-0.50 





1.00 


5.00(6 


H 2 S + H+ -4 H 2 S+ + H 


0.38 (-8) 


-0.50 





1.00 


5.00(6 


H 2 CS + H+ -y H 2 CS+ + H 


0.64 (-8) 


-0.50 





1.00 


5.00(6 


OCS + H + -)• HS + + CO 


0.65 (-8) 


-0.50 





1.00 


5.00(6 


S0 2 + C + -y SO+ + CO 


0.20 (-8) 


-0.50 





1.00 


5.00(6 


S0 2 + H 3 + -> HS0 2 + + H 2 


0.37 (-8) 


-0.50 





1.96 


5.00(6 


O + SO -> so 2 


0.32(-15) 


-1.60 





1.00 


5.00(6 


O + HCS -> OCS + H 


0.50 (-10) 








1.00 


5.00(6 


S + CO -)• OCS 


0.16(-16) 


-1.50 





1.00 


5.00(6 


S + CH 3 H,CS + H 


0.14 (-9) 








1.00 


5.00(6 


S0 2 + C -4 CO + so 


0.70 (-10) 








1.00 


5.00(6 


OH + SO -> S0 2 + H 


0.86 (-10) 








1.00 


5.00(6 


C 3 S+ + e" -> C 2 S + C 


1.00 (-7) 


-0.50 





1.00 


5.00(6 


C 3 S+ + e" -> CS + C 2 


1.00 (-7) 


-0.50 





1.00 


1.20(4 


C 4 S + + e--i-C,S + C2 


1.00 (-7) 


-0.50 





1.00 


5.00(6 


H 3 CS+ + e" H 2 CS + H 


0.30 (-6) 


-0.50 





1.00 


5.00(6 


H 3 S + + e~^HiS + H 


0.30 (-6) 


-0.50 





1.00 


5.00(6 


HC 2 S + + e--i-C,S + H 


0.15 (-6) 


-0.50 





1.00 


5.00(6 


HCS+ + e" -> CS + H 


0.1 8 (-6) 


-0.57 





1.00 


5.00(6 


HOCS+ + e ~ -> CS + OH 


0.20 (-6) 


-0.50 





1.00 


5.00(6 


OCS+ + e ~ -y CS + 


0.48 (-7) 


-0.62 





1.00 


1.20(4 


H 2 CS + + e" -> H 2 CS 


0.11 (-9) 


-0.70 





14.60 


5.00(6 


H 2 S + + e" -> H 2 S 


0.11 (-9) 


-0.70 





1.00 


5.00(6 


H 3 CS+ + e" -> CS + H + H 2 


0.30 (-6) 


-0.50 





1.00 


5.00(6 


C 3 S+ + grain(-) -> C 2 S + C + grain(O) 


0.33 








3.82 


5.00(6 


HOCS+ + grain(-) -y CS + OH + grain(O) 


0.50 








3.82 


5.00(6 


H 3 CS + + grain(-) -y H 2 CS + H + grain(O) 


0.50 








3.82 


5.00(6 


HC 2 S + + grain(-) -y C 2 S + H + grain(O) 


0.50 








1.96 


5.00(6 



In this Section we analyze in detail chemical and mix- 
ing processes responsible for the evolution of sulfur-bearing 
species in protoplanetary disks. The chemistry of sulfur- 
bearing molecules is the least understood in astrochemistry. 
As other heavy elements, sulfur is depleted from the gas, but 
the magnitude of the effect and in which form sulfur is bound 



is not yet clearly determined (e .g., lLoddersll2003t iFlvnn et all 
2006: ^Goicoechea et alj [2006). Moreover, many processes 
and reaction rate data for the S-bearing species are lacking 
accurate laboratory measurements or theoretical foundations, 
leadin g to considerable uncertainties in modeled abundances 
(e.g.. l Vasvunin et alj|2008t IWakelam et al.ll2010bl) . Thus, re- 
sults in this subsection should be taken with special care. So 
far there are only two sulfur-containing molecules detected in 
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Figure 11. The same as in Fig.[7Jbut for the S-containing species. Results are shown for CS, C2S, C2S ice, C3S ice, tbS, OCS, H2CS, and SCh. 



disks, namely, CS ( e.g jDutrev et al.|[l997l) and, recently, SO 
dFuente et alj|2010h . In the interstellar medium and the shells 
of AGB stars other s pecies have been disco vered, e.g. C 2 S, 
C 3 S and H,CS (e.g..lCernicharo et alJ1987t lYamamoto etail 
I1987L llrvine et al.lll989l) . So, here we restrict ourselves with 
discussion of neutral species only. 

There are 2 steadfast neutral S-bearing molecules (H2S and 
HCS), 5 ices (HS, H 2 S, H 2 CS, S, and S 2 ), and no ions (Ta- 
ble [3). Their column densities in the laminar and the fast 
2D-mixing model differ by a factor of < 3. The sensitive 
sulfur-bearing species include 7 gas-phase molecules (e.g., 
CS, H 2 CS, H 2 S 2 , HS, S), 2 ions (S + and SO + ), and 7 ices 
(e.g, CS, C 3 S, C 4 S, H 2 s 2 , SO; Table©. Their column densi- 
ties are changed by up to 2 orders of magnitude by the turbu- 
lent transport. Among hypersensitive species (Table [7} are 4 
gas-phase molecules (C 2 S, C3S, SO, and S0 2 ), no ions, and 
3 ices (C 2 S, NS, and SiS). The turbulent diffusion alters their 
column densities by up to a factor of 7 000 (S0 2 ice). Simi- 
larly to hydrocarbons and oxygen-containing molecules, com- 
plex S -chains with multiple carbon atoms or oxygen are more 
strongly affected by the turbulent transport than chemically 
simpler species. Unlike the case of C- and N-bearing species, 
there is no trend that the concentrations of the sulfur-bearing 
ices are stronger altered by the transport than the abundances 
of their gas-phase analogs. 

In Fig. [TT] the distributions of the relative molecular abun- 
dances and column densities at 5 Myr of CS, C 2 S, C 2 S ice, 
C 3 S ice, H 2 S, OCS, H 2 CS, and S0 2 calculated with the lam- 
inar and the 2D-mixing models are shown. In the laminar 
model the relative abundance distributions of the gas-phase 
species show a 3-layered structure, with peak concentrations 
in the molecular layers at w 0.8 - 1 H r . In contrast, the C 2 S and 
C3S ices are more concentrated toward the disk midplane and 
the bottom of the molecular layer, with typical abundances 
of 10~ 12 - 10~ 10 . The molecular layers of gas-phase sulfur- 



bearing molecules have various thicknesses, from ss 0.05//, 
for C 2 S and up to ~ 2H r for CS, and are not as pronounced 
as for carbon-, oxygen- and nitrogen-bearing species (see 
Figs. [HT-fTOb. Their typical peak relative abundances are only 
10~ 12 - 10~ 10 , though in the outer molecular layer the H 2 S 
abundances are as high as ~ 10~ 8 . 

Among the above molecules, only H 2 S is steadfast (see 4th 
panels in Fig. [TTJ. The mixing enhances the H 2 S concentra- 
tion in the upper molecular layer at z > l.3H r that does not 
contribute to the total column density (compare 1st and 3rd 
panels in the Figure). The CS, C 2 S, C 2 S ice, C 3 S ice, OCS, 
H 2 CS, and S0 2 are sensitive and hypersensitive to the mixing. 
The corresponding molecular layers are vertically broadened 
up to r; 2H r , and strongly molecularly enriched by several 
orders of magnitude. 

To better understand these results, we investigate the evolu- 
tion of CS, C 2 S, C 2 S ice, C 3 S ice, H 2 S, OCS, H 2 CS, and S0 2 
in the two disk vertical slices at r = 10 and 250 AU (the lam- 
inar chemical model). The most important reactions respon- 
sible for the time-dependent evolution of their abundances in 
the midplane, the molecular layer, and the atmosphere are pre- 
sented in Table 1 121 both for the inner and outer disk regions. 
The final list contains top 20 reactions per region (midplane, 
molecular layer, atmosphere) for the entire 5 Myr time span, 
with all repetitions removed. 

The chemical evolution of hydrogen sulfide in the atmo- 
sphere begins with radiative association of ionized sulfur and 
H 2 , leading to H 2 S + that slowly radiatively recombines into 
H 2 S. The major destruction reactions in the atmosphere is 
photodissociation and rapid ion-molecule reaction with C + 
forming protonated CS and H (75%) or H 2 S + and C (25%), 
and charge transfer reaction with H + . In the molecular layer 
and midplane hydrogen sulfide is initially produced as in the 
atmosphere, but after 1-100 years the top formation pathway 
for gas-phase H 2 S is a direct surface recombination of HS and 
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H ices. This is caused by a relatively large binding energy of 
H2S that prevents H2S ice to return to the gas phase unless 
temperatures exceed about 45-50 K. The HS radical is formed 
via neutral-neutral reaction of S with OH, and also via di- 
rect surface recombination of H and S. The major destruction 
pathways for H2S in the molecular layer and midplane are ac- 
cretion onto dust grain surfaces, ion-molecule reactions with 
HCO + and H| (leading to H3S" 1 "), and charge transfer reaction 
with atomic hydrogen. The characteristic timescale of the 
H2S evolution is short in the midplane and molecular layer, 
< 10 - 10 4 years since its initial production and freeze-out are 
relatively fast. Later, accretion is balanced out by the direct 
surface recombination, while desorption of the H2S is not sig- 
nificant. In the atmosphere the H2S reaches a quasi steady- 
state within 1 4 — 1 6 years, which is determined by the slow 
evolution of H2, Hj, and H3 (see discussion in Sect. 14. II ). 

The chemical evolution of CS proceeds entirely in the gas 
phase. The carbon monosulfide is synthesized by dissociative 
recombination of HCS\ C 2 S\ C 3 S\ HOCS + , and H 3 CS + on 
electrons or negatively charged grains in the inner midplane, 
and via desorption of the CS ice at T > 35 -40 K. The major 
destruction pathways are photodissociation in the atmosphere, 
depletion in the molecular layer and the midplane, slow en- 
dothermic oxidation reaction (CS + O — > CO + S), charge 
transfer with ionized hydrogen atoms, and ion-molecule reac- 
tions with primal ions, e.g. HCO + , H30 + , and H3, leading to 
HCS + (Table [T2"ii. In turn, the protonated carbon monosulfide 
is initially synthesized through an ion-molecule destruction 
of H2 by C + . The C2S" 1 " is produced either by an ion-molecule 
reaction of ethynyl radical with ionized atomic sulfur or via 
charge transfer reactions of C2S with H + and C + . The C3S + 
chemistry is similar to that of C2S + , with the production path- 
way that involves C3H instead of ethynyl. The HOCS + is 
protonated OCS, and is formed via ion-molecule reactions of 
OCS with dominant polyatomic ions (H3, HCO + , etc.). Fi- 
nally, HsCS + is synthesized via protonation of H2CS and via 
ion-molecule reactions of ionized sulfur or HS + with methane. 
The steady-state for CS is attained at about 10 2 - 10 4 years 
in the lower molecular layer/upper midplane, and later in 
other disk parts. While the CS chemistry does not include 
surface processes directly, it is related with the evolution of 
C2S + , C3S + , and HOCS + , and, thus, their chemical "parents", 
namely C2S, C3S, and OCS. The latter molecules are partly 
synthesized via surface reactions, making their evolution slow 
in the cold midplane. 

The chemical evolution of all the C„S species in the model 
(n < 4) is well coupled. For example, like in the case of CS, 
major formation channels for gas-phase C2S include disso- 
ciative recombination of C 3 S + , C 4 S + , HC 2 S + , and HC 3 S + on 
negatively charged grains in the inner midplane and electrons. 
The HC2S" 1 " and HC3S" 1 " ions are formed later in the disk by 
protonation of C2S and C3S. The binding energy adopted for 
C2S is 2 700 K, and higher for more massive carbon chain 
sulfides, which precludes effective thermal desorption of their 
ices in the disk midplane and the molecular layer beyond 
10 AU. The major destruction channels for C2S are photodis- 
sociation, accretion onto dust grains, charge transfer reactions 
with H + and C + , and protonation by abundant polyatomic ions 
(leading to HC2S" 1 "). The evolution of the C2S ice is signifi- 
cantly simpler. It involves only 3 major processes: accretion 
of the gas-phase C2S, surface reaction of frozen atomic carbon 
with the C2S ice producing the C3S ice, and photodissociation 
of the C2S ice by the CRP/X-ray-induced UV photons (lead- 



ing to C and CS ices). The evolution of the C3S ice includes its 
surface synthesis from C2 and C, accretion of gas-phase C3S, 
and photodestruction to the C2 and CS ices. Consequently, the 
characteristic chemical timescales for gas-phase and, particu- 
larly, for surface C„S molecules are long, T C h e m ^ 10 6 years. 

Thioformaldehyde, H2CS, like CS and H2CO, is synthe- 
sized in the gas. Its major production routes consist of neutral- 
neutral reaction between S and CH3, and dissociative recom- 
bination of protonated H2CS at later times. The key destruc- 
tion routes are photodissociation, accretion onto dust grains 
at T < 45 - 50 K, charge transfer reaction with ionized atomic 
hydrogen, ion-molecule reaction with ionized atomic carbon, 
and proton addition in reactive collisions with HCO + , H3, etc. 
The distribution of the H2CS characteristic timescale over the 
disk closely resembles T c h em (CS). 

The chemical evolution of carbonyl sulfide (OCS) begins 
with its formation upon oxidation of HCS and via radiative 
association of S and CO. At later times OCS is also produced 
by the direct surface recombination of the S and CO ices, and 
the CS and O ices. The destruction of OCS involves pho- 
todissociation in the disk atmosphere, accretion onto dust in 
the molecular layer and the midplane, and ion-molecule re- 
actions with ionized C and H (leading to either CS + and CO 
or HS + and CO). The protonation of OCS by the polyatomic 
ions leads to HOCS + , which dissociatively recombines into 
OCS and H (50%) or CS and H (50%). The characteristic 
timescale for carbonyl sulfide is < 10 2 — 10 6 years in the disk 
midplane, ~ 10 3 — 10 4 years in the lower molecular layer, and 
> 1 4 — 1 6 elsewhere in the disk. 

Finally, sulfur dioxide (SO2) is synthesized by exothermic 
neutral-neutral reaction between SO and OH, slow radiative 
association of SO and O, and the same reaction catalyzed 
by dust grains, with SO2 directly desorbed back to the gas 
(with 5% probability). Its major removal channels are pho- 
todissociation, freeze-out in disk regions with T < 60- 80 K, 
neutral-neutral reaction with atomic carbon (leading to SO 
and CO), and ion-molecule reaction with C + (producing SO + 
and CO). Later a quasi-equilibrium cycling of SO2 prevails, 
which includes protonation of sulfur dioxide by abundant 
polyatomic ions (H3, HCO + ) followed by dissociation recom- 
bination back to SO2 with the 66% probability. In turn, sulfur 
monoxide (SO) is produced by slow combustion of S, neutral- 
neutral reactions between S and OH and O with HS, and direct 
surface reaction of the HS and O ices. Consequently, the char- 
acteristic chemical timescale of SO2 is governed by the slow 
evolutio n of molecular oxygen due to surface processes (see 
Section 14.31 ), and exceeds 1 Myr in the disk midplane, and 
10 4 - 10 5 years in the molecular layer. 

Since chemical network for S-bearing species, especially 
complex chains, in our model is more limited than that for the 
N- and C-containing molecules, the turbulent transport en- 
hances abundances and column densities even for transient 
S-species. All species considered above either formed via 
surface reactions or their synthesis involves surface-produced 
radicals, and their chemical evolution proceed slower than the 
mixing (see Fig. ITTb . As for many molecules in the model, 
the turbulent transport expands their molecular layers in ver- 
tical direction. Complex S-ices are more readily produced in 
the mixing case as their formation is based on heavy radicals, 
which become mobile at T > 20-40 K. Also, in the mixing 
model photodissociation of these ices is more effective (see 
a peak in the relative abundances of C2S ice at r ~ 100 AU, 
z ~ \ .5H r ). The diffusive transport most strongly increases 
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concentration of SO2 molecules by producing more molecu- as its characteristic evolutionary timescale is shorter than the 
lar oxygen needed for its synthesis. For SO2, its radial distri- mixing timescale in the disk region where it is most abundant 
bution becomes sensitive to the transport. The least affected (see Tables [3]-|4]l. 
species is H2S (as well as HCS, and ices of S, HS, and H2S) 

4.6. Complex organic molecules 



Table 13 

Key chemical processes: organic species 



Reaction 


a 





7 


tin in 


traax 




[(crn^s" 1 ] 




[K] 


[yr] 


[yr] 


CH 2 CO ice + hf C RP — > CH 2 ice + CO ice 


9.15(2) 








1.00 


5.00(6) 


CH2CO ice + hi/cRP — ► C 2 ice + H2O ice 


4.07 (2) 








1.96 


5.00(6) 


CH3CHO ice + h^cRp — > CH 3 ice + HCO ice 


5.25(2) 








1.96 


5.00(6) 


CH3CHO ice + hi/ C Rp -> CH 4 ice + CO ice 


5.25(2) 








1.96 


5.00(6) 


CH 3 OH ice + hi/ C Rp — ► CH 3 ice + OH ice 


1.50(3) 








1.00 


5.00(6) 


CH3OH ice + hi^cRP -> H 2 CO ice + H 2 ice 


3.17(3) 








1.00 


5.00(6) 


C2H5OH ice + hi/ C Rp — > CH3CHO ice + H 2 ice 


6.85(2) 








2.34(4) 


5.00(6) 


HCOOH ice + hz/cRp — > CO2 ice + H ice + H ice 


6.50(2) 








1.00 


5.00(6) 


HCOOH ice + Im/crp — ► HCO ice + OH ice 


2.49 (2) 








3.82 


5.00(6) 


HNCO ice + hfcRp -> NH ice + CO ice 


6.00(3) 








1.00 


5.00(6) 


CH3OH ice + UV ->• H 2 CO ice + H 2 ice 


0.72 (-9) 





1.72 


1.00 


5.00(6) 


C 2 H 5 OH ice + UV — > CH3CHO ice + H 2 ice 


0.13 (-9) 





2.35 


4.18(2) 


5.00(6) 


HCOOH ice + UV — >■ HCO ice + OH ice 


0.28 (-9) 





1.80 


1.00 


5.00(6) 


HNCO + Ii^crp — > NH + CO 


6.00(3) 








1.00 


5.00(6) 


CH3CHO + UV -> CH 3 CHO+ + e - 


0.26 (-9) 





2.28 


1.00 


5.00(6) 


CH3CHU + UV — > CH4 4- CO 


0.34 (-9) 





1.52 


1.00 


5.00(6) 


CH3CHO + UV — >■ CH3 + HCU 


0.34 (-9) 





1.52 


1.00 


5.00(6) 


T T/~*/~\/~\T T 1 T T\ T ! T 1 / M T+ i — 

HCOOH + UV — ► HCOOH + e 


0.17 (-9) 





2.59 


1.00 


5.00(6) 


HCUUH + UV — > HCU + OH 


0.28 (-9) 





1.80 


1.00 


5.00(6) 


CH2CO + grain — > ChbCO ice 


1.00 








1.00 


5.00(6) 


HCUUH + grain — > HCOOH ice 


1.00 








1.00 


5.00(6) 


T |MC Y1 . _, ■ , ITMfA - 

HJNCO + grain — > HNCO ice 


1.00 








1.00 


5.00(6) 


Cri2CU ice — > CrnCO 


1 .UU 


u 


Z.ZU {3} 


1 .uu 




CH3CHO ice -> CH3CHO 


1 nn 
1 .uu 


u 


Z.5/ (3) 


1 nn 
1 .uu 


J . UU ) 


rirN^U ice — ? rllNx^U 


1 nn 
l.UU 


n 
u 


z.oj (_T ; 


1 nn 
1 .uu 


j.uu^o; 


11 . IT f ill \ri* ^ fT-T^fTR \rp 

rl ice -r L^rlTVJrl ice — r V n 1 \.> 11 ice 


1 .uu 


n 
u 


n 
u 


1 nn 
1 .uu 


J . UU (p) 


H ice + HCtO ice — > CH 2 CO ice 


1.00 








1.00 


5.00(6) 


H ice + OCN ice -> HNCO ice 


1.00 








1.00 


5.00(6) 


H ice + OCN ice -> HNCO 


1.00 








1.00 


5.00(6) 


OH ice + CH 3 ice -> CH3OH ice 


1.00 








1.00 


5.00(6) 


OH ice + HCO ice -> HCOOH ice 


1.00 








1.09(2) 


5.00(6) 


HCOOH + H+ -> HCOOH+ + H 


0.28 (-8) 


-0.50 





1.00 


5.00(6) 


HNCO + H+ -> NH 2 + + CO 
CH3CHO + H 3 + -> CH 3 CH 2 0+ + H 2 


0.15 (-7) 


-0.50 





1.00 


5.00(6) 


0.62 (-8) 


-0.50 





1.00 


5.00(6) 


CH3CHO + HCO+ -> CH 3 CH 2 + + CO 


0.25 (-8) 


-0.50 





1.00 


5.00(6) 


HCOOH + H3 + -> HCO+ + H 2 + H 2 


0.23 (-8) 


-0.50 





1.00 


5.00(6) 


HCOOH + HCO+ -> CH 3 2 + + CO 


0.13 (-8) 


-0.50 





1.00 


5.00(6) 


O + C1H5 H. CH3CHO + H 


0.13 (-9) 








1.00 


5.00(6) 


OH + H 2 CO -> HCOOH + H 


0.20(-12) 








1.00 


5.00(6) 


CH 3 CH 2 0+ + e" -> CH3CHO + H 


0.15 (-6) 


-0.50 





1.00 


5.00(6) 


CH 3 2 + + e - -> HCOOH + H 


0.15 (-6) 


-0.50 





1.00 


5.00(6) 


CH 3 CH 2 + + grain(-) -> CH3CHO + H + grain(O) 
CH30 2 + + grain(-) -» HCOOH + H + grain(O) 


0.20 








55.90 


5.00(6) 


0.50 








1.00 


5.00(6) 


CH3CHO + grain -> CH3CHO ice 


1.00 








1.00 


5.00(6) 


CH3OH + grain -> CH3OH ice 


1.00 








1.00 


5.00(6) 



The presence of amino acids and other complex organics 
in the early Solar system is a well-established through de- 
tailed mass-spec trometry of carbonaceous meteorites (e.g., 
iGlavin et al.ll20ldl) . via analysis of the cometary dust sam- 
pled by the Giotto spacecraft in the Halley comet (e.g., 
Jessberger et al. 1988), and the recent i dentification of g lycine 
in the Stardust cometary dust samples (Elsila et al. 2009J). The 
ground-based search for simple gas-phase organic species, 
such as methanol and formic acid, in nearby protoplanetary 
disks with radiotelescopes have so far been frui tless, though 
formal dehyde has be en detected in DM T au (e.g. JDutrev et al.l 
I2007al) and LkCa 15 dAikawa et al.l2003l) . The simple organic 



ices, e.g. HCOOH ice, has been identified in t he Spitzer spec- 
tra of several low-mass Class I/II objects (e.g., Zasow ski et alj 
2009). In contrast, rich organic gas-phase molecules have 
been detected in hot massive cores and corinos, where condi- 
tions are appropriate for st eady sublimation and production of 
heavy complex ices (e.g., Belloche et alj|2008l : iGarrod et al.l 

2008) . Nowadays synthesis of complex organics in cosmic 
objects is thought to proceed solely via surface processes, 
such as photoproduction of reactive radicals by high-energy 
irradiation of complex precursors an d their surface recombi- 
nation at elevated temperatures, (e.g. jHerbst & van Dishoeckl 

2009) . In this Section we analyze in detail chemical processes 
responsible for the formation and destruction of the complex 
organic species in protoplanetary disks. These are meant to 
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be polyatomic neutral molecules consisting of at least several 
H, C, and O atoms, and heavier elements. 

There is a single organic species insensitive to the tur- 
bulent diffusion, namely, formaldehyde (see Sect. I4.31 >. 
The transport-sensitive organic species include 2 gas-phase 
molecules (CH 2 CO and HNCO) and 5 ices (CH 2 CO, HNCO, 
CH 3 OH, H 2 CO, NH 2 CHO; Table 0. Their column densi- 
ties are changed by up to 2 orders of magnitude. One gas- 
phase molecule (HCOOH) and 3 ices (HCOOH, H 2 C 3 0, and 
CH3CHO) are among hypersensitive species in Table [7] The 
turbulent diffusion alters their column densities by up to a fac- 
tor of 3 000 (solid formic acid). The chemical evolution of all 
the considered organic species is influenced by the transport 
as their major production and removal routes require surface 
processes. 

In Fig. [12] the distributions of the relative molecular abun- 
dances and column densities at 5 Myr of HCOOH, HCOOH 
ice, CH3OH, HNCO, HNCO ice, CH 3 CHO, CH3CHO ice, 
and CH2CO calculated with the laminar and the 2D-mixing 
models are presented. In the laminar model the relative 
abundance distributions of the gas-phase species show a 3- 
layered structure, with very narrow molecular layers located 
at « 0.8 - 1 H r . Note that their abundances and column den- 
sities are substantially lower than those of simpler C- or O- 
bearing species considered above, so that some heavy organic 
species are not included in Tables EH7] The HCOOH, HNCO, 
and CH2CO molecular layers are wider beyond ~ 200 AU, 
and have higher abundances. Remarkably, abundance distri- 
butions of complex ices are maximum either at the bottom of 
the molecular layer (HNCO ice, HCOOH ice) or in the in- 
ner warm midplane (CH3CHO ice, HCOOH ice). The overall 
tendency can be easily explained as heavy ices are hard to 
evaporate thermally, so photoevaporation is necessary, while 
gas-phase organic molecules are rather photofragile. For ef- 
ficient production of heavy ices surface mobility of radicals 



(O, C, CH, etc.) necessitates warm temperatures (T > 30 K) 
and/or photoprocessing of their precursor molecules by the 
X-ray /CRP-induced UV photons. 

To better understand these results, we investigate the evolu- 
tion of HCOOH, HCOOH ice, CH 3 OH, HNCO, HNCO ice, 
CH3CHO, CH3CHO ice, and CH 2 CO in two disk vertical 
slices at r = 10 and 250 AU (the laminar chemical model). The 
most important reactions responsible for the time-dependent 
evolution of their abundances in the midplane, the molecular 
layer, and the atmosphere are presented in Table Q~3] both for 
the inner and outer disk regions. The final list contains only 
top 25 reactions per region (midplane, molecular layer, atmo- 
sphere) for the entire 5 Myr time span, with all repetitions 
removed. 

The chemical evolution of formic acid (HCOOH) starts 
with a single production channel via dissociative recombi- 
nation of CH3O2 either on electrons or negatively charged 
grains (Table [T3l , In turn, protonated formic acid is pro- 
duced by radiative association of HCO + and H 2 0, and by ion- 
molecule reaction of methane with ionized molecular oxy- 
gen. The primal removal channels for HCOOH are photodis- 
sociation and photoionization, depletion onto dust grains at 
T < 100- 120 K, charge transfer with ionized atomic hydro- 
gen, and protonation by abundant polyatomic ions (HCO + , 
H3, H3O" 1 "). The HCOOH ice is essentially synthesized by ac- 
cretion of gas-phase formic acid, and destroyed by secondary 
UV photons in the disk midplane and molecular layer, and 
photoevaporation. The surface as well as gas-phase formation 
of HCOOH via neutral-neutral reaction of OH and H 2 CO is 
only a minor channel in the laminar model. Consequently, the 
distribution of evolutionary timescales for gas-phase and sur- 
face HCOOH is similar and dominated by r c h e m ^ 10 5 years. 

Gas-phase methanol (CH3OH) is produced by direct sur- 
face recombination of frozen H and CH 2 OH as well as frozen 
OH and CH3, evaporation of methanol ice, dissociative re- 
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combination of protonated methanol, and dissociative recom- 
bination of H5C2O2 (minor route). Its key removal pathways 
are surface accretion in the disk regions with T < 100- 120 K, 
photodissociation and photoionization, and protonation by the 
dominant polyatomic ions. On dust surfaces methanol is 
produced via a sequence of hydrogenation reactions, starting 
from the CO ice, and is photodissociatied by the UV radiation. 
Similarly to HCOOH, the characteristic chemical timescale 
for methanol is fully controlled by the slow surface chemistry 
(Tchem > 10 5 years). 

The chemical evolution of the isocyanic acid (HNCO) is 
also dominated by surface processes. In the gas-phase it is 
produced either by evaporation of HNCO ice at T < 50-60 K 
or by direct surface recombination of surface H and OCN. The 
major gas-phase destruction pathways are accretion to grains 
in the midplane and the molecular layer, and ion-molecule 
reaction with H + (leading to NH£ and CO; Table [B]). The 
HNCO ice is produced by the surface reaction involving H 
and OCN, re-accretion of the HNCO gas, and destroyed by 
the UV. Their characteristic timescales alike those of HCOOH 
and methanol. 

The chemistry of acetaldehyde (CH3CHO) involves gas- 
phase production by reactive collisions between O and C2H5, 
accretion to dust grains at T < 50 - 60 K, and desorption at 
T < 60 K, and removal via charge transfer reactions with C + 
and ion-molecule reactions with C + , H3, and H3O" 1 ". At later 
times acetaldehyde is re-produced from its protonated analog 
by dissociative recombination (albeit with low probability). 
The frozen acetaldehyde is synthesized by accretion of gas- 
phase CH3CHO, via surface recombination of CH3 and HCO 
ices, and destroyed by the UV-dissociation. As for other com- 
plex organics, the chemical timescale for acetaldehyde is typ- 
ically longer than 1 Myr in the disk. 

Finally, the gas-phase ethenone (CH2CO) is produced via 
oxidation of C2H3, direct surface recombination of the H 
and HC2O ices, and dissociative recombination of protonated 
ethenone. Thermal evaporation of ethenone ice is effective 
when dust temperatures exceed ~ 100 K. Key removal chan- 
nels include photodissociation and photoionization, freeze- 
out in the midplane and the molecular layer, charge transfer 
reaction with C + and H + , and, at later times, protonation by 
polyatomic ions. The chemical timescale for CH2CO exceeds 
10 5 - 10 6 years. 

Not surprisingly that turbulent mixing enhances abun- 
dances and column density of these organic species, given 
their long evolutionary timescales governed by surface reac- 
tions and photodissociation of ices by the CRP- and X-ray- 
induced secondary UV radiation field. The HCOOH and the 
HCOOH ice are more abundant since water ice abundances 
and HCO + abundances are increased by the mixing, leading 
to more efficient production of the parental ion, CH3O" 1 ". Both 
for gas-phase and solid formic acid radial transport is impor- 
tant. The effect is less pronounced for HNCO as its precursor 
species, OCN, is enhanced by transport only by less than 1 
order of magnitude. The HNCO ice, produced in the lower 
part of the molecular layer, is transported by diffusion to the 
cold midplane where it cannot be synthesized otherwise. Rel- 
ative abundances of gas-phase acetaldehyde are greatly en- 
riched by the disk mixing in the middle of the molecular layer 
(z « 1.5 H r ), within the inner 100 AU. As its synthesis pro- 
ceeds via surface recombination of CH3 and HCO, in the lam- 
inar model appropriate conditions are only met in the very in- 
ner midplane/molecular layer. In the fast mixing model larger 



quantities of solid acetaldehyde can be accumulated as more 
icy grains reach the warm disk regions. The vertical mixing 
brings the CH3CHO ice in the inner midplane to the molec- 
ular layer, where it photoevaporates, and then radially trans- 
ported to larger distances. The ethenone abundances are ver- 
tically broaden by the mixing, and moderately enriched as the 
ethenone ice production is enhanced by the turbulent mixing. 

5. DISCUSSION 

5.1. Comparison with previous studies and future 
developments 

In this Section we discuss the results and the drawbacks 
of our model in the context of other studies of the chemo- 
dynamical evolution of protoplanetary systems. 

(Models of the inner Solar nebula: the role of advective 
transport) Historically, the interest to this topic has been ini- 
tiated by the cosmochemical community studying the initial 
stages of the formation of the Solar system. The presence of 
high-temperature condensates, like CAIs and crystalline sili- 
cates, in pristine chondritic meteorites (formed within a few 
AU from the Sun) and comets (formed within 10-20 AU) 
as well as almost perfect isotopic homogeneity of the in- 
ner Solar nebula (r < 10-20 AU) at a bulk level both re- 
quire an efficient mixing mechanism with a transport speed of 
< 1 AU per 5 10 4 yr or 10 cms" 1 (e.g., iWooden et alJl2007t 
iBrownlee et~aLll2008b ICieslall2009h . According to the current 
knowledge, the Solar nebula and other protoplanetary disks 
represent an outcome of the viscous evolutionary stage that 
follows the initial collapse of a molecular cloud core (e.g., 
ICassen & MoosmanlfT98lt ICameronlfT99l lLarsoril2003l 

The ability of turbulent or advective transport to cause ra- 
dial mixing of materials in the inner early Sol ar system has 
been proposed and in vestigated by Morfill (e.g., Morfil ll983t 
Morfil l & Volkl 19841) . with a simple ID analytical disk model 
and passive tracers. In a similar manner, 2D radial mix- 
ing of gas eous and solid wa ter in the inner nebula has been 
studied by ICvr et al] dl998l) . In a long series of papers the 
group of Prof. H.-P. Gail (Heidelberg University) has investi- 
gated various aspects of the chemo-dynamical evolution of 
the inner protosolar nebula (1-10 AU), progressing from 
stationary accretion disk models and crude gas-phase chem- 
istry toward a self-consistent 2D radiative-hydrodynamical 
mod el with the coupled C-, H-, O- gas - phase chemistry 
(e.g IBauer et all 19971; IFinocchi et all 1 997t iGaill 1 998l|200 1L 
| 2002tlWehrstedt & GaiJ2002tlGaiJ2004tlKeller & Gaill2004t 
iTscharnuter & Gaill 120071) . In IBauer et al.l (119971) the gas- 
phase C-, H-, N-, O-chemistry driven by dust destruction 
and evaporation of ices in the presence of slow radial trans- 
port has been mod eled, utilizing the semi-analytical one- 
zone disk model of iDuschl et al. (1996). The authors have 
found that silicates tend to evaporate under equilibrium con- 
ditions, whereas carbon dust is slowly destroyed by combus- 
tion by OH, and that at r < 1 AU the CO concen trati on is 
higher than that of water. In IFinocchi et al.l (119971) and IGaill 
(2001) the one- zone station ary disk model has been con sid- 
ered, whereas in IGaill (120021) and lWehrstedt & Gaill (12002b the 
one-zone time-dependent a-disk model with a larger chemi- 
cal network of 106 species, annealing of amorphous silicates, 
and more detailed description of carbon dust oxidation by O 
and OH has been used. The major finding is that the ra- 
dial transport enriches the outer, r > 10 AU nebular regions 
with methane and acetylene produced by oxidation of carbon 
dust at r < 1 AU, as observed in comets of Hyakutake and 
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Hale-Bopp. The advection also transports cryst alline silicates 
to the 10-20 AU region In th e later papers (Keller & Gaill 
120041; ITscharnuter & Gail 2007), a more accurate model has 
been developed, with a self-consistent 2D prescription of ad- 
vective transport and turbulent diffusion, grey radiative trans- 
fer, proper dust and gas opacities, and disk hydrodynamical 
evolution with parametric viscosity. The model was coupled 
to a limited gas-phase C-, H-, O-chemistry involving nearly 
90 neutral-neutral reactions. They have found that quasi- 
stationary accretion flows develop a 2D structure such that in 
the disk midplane gas moves outward, carrying out the angu- 
lar momentum, whereas inward mass accretion goes through 
the surface layers. The radial advection dominates the diffu- 
sive mixing in outer disk regions at r > 5 AU. Thus, chemical 
species produced by the "warm" chemistry in the inner nebula 
can reach its outer region, where they are intermixed by tur- 
bulence with surrounding matter and freeze-out. They have 
considered neither X-ray- nor UV-driven processes, no ion- 
molecule and surface chemistry, and assumed sublimation- 
evaporation equilibrium for ices. Their key result is that the 
slow water-shift reaction transforming CO and water into CO 2 
and H2 at high pressure and temperature is important for the 
oxygen chemistry. The H2 and CO evolution is insensitive to 
transport, while abundance distributions of e.g. O, O2, CO2 
are altered by the disk dynamics within 300 years of the evo- 
lution. 

Contrary to these studies, our modeling is related to outer 
nebular regions beyond 10 AU and does not include advective 
transport. However, in the presence of steep chemical gra- 
dient for a molecule, radial diffusive mixing is able to trans- 
port it from 100 AU to 1 AU within about 1 — 2 1 6 years, and 
from 1 AU to 10 AU within about < 3 10 5 years, which is still 
less than a typical disk lifetime of several Myrs (see Eq. (|7} 
and Fig. [TJ. Our chemical network incorporates much more 
species and reactions, including surface chemistry, and the 
high-energy processing of the gas and icy mantles. Nonethe- 
less, the calculated H2 and CO column densities are not sen- 
sitive to the diffusive mixing in our mo del, similar to what 
was found bv lfscharnuter & Gaill (120071) (Fig.|9]and TableE). 
Also, in agreement with Gail et al. results, we show that the 
synthesis of hydrocarbons, including CH4 and C2H2, is inten- 
sified by the mixing, particularly in the zone of comet forma- 
tion around 10-20 AU (Fig. [8). Furthermore, CO2 ice, which 
is produced via surface reactions in the warm inner midplane 
and can be considered as an example of a "warm" species, 
is transported outward by the 2D-turbulent mixing. In addi- 
tion, we also found that oxidation chemistry is important for 
the evolution of O-bearing species, thoug h we do not account 
for th e carbon dust destruction, unlike iTscharnuter & Gaill 
d2007h . 

(Models of protoplanetary disks: the role of advective 
transport) Next, we focus on disk models that has been used 
to simulate chemical e volution of protoplan etary disks in the 
astrophysical context. Aikaw aet al.l (11999b have util i zed an 
isothermal a-disk model of iLvnden- Bell & Pringlel (fl974h 
with M= 10~ 8 M Q yr _1 , and calculated its density structure 
assuming hydrostatic equilibrium in vertical direction. They 
have used the early version of the UMIST 95 chemical ratefile 
(Milla ret al.lll997l) . considered sticking and thermal desorp- 
tion of molecules and surface formation of molecular hydro- 
gen, and grain re-charging (about 250 species and 2 300 reac- 
tions). No X-ray-driven and surface chemistry, other than H2 
formation, has been included. 



These authors have modeled the chemical evolution in the 
presence of steady inward accretion, using both atomic and 
molecular initial abundances, and found similar results. We 
also find that the difference in the calculated abundances be- 
tween the adopted "low metals" atomic initial abundances 
and those from a mo lecular cloud is negligible in our model. 
lAikawa etaT] {[999) have found that it takes about 3 Myr to 
transport disk material from 400 AU to 10 AU (similar to the 
timescale of radial mixing in our model). The inward advec- 
tion leads to higher concentrations of heavy hydrocarbons at 
< 20 AU, whereas methane is a dominant hydrocarbon in the 
outer disk region. In our laminar model methane prevails over 
heavier hydrocarbons in the entire disk, while in the presence 
of radial and vertical mixing a substantial fraction of CH4 is 
converted i nto heavy carbon ch ains (especially at 10-50 AU; 
see Fig. [8). lAikawa et al.l (1 19991) have concluded that the ra- 
dial transport leads to simultaneous existence of the reduced 
(e.g. CH4) and oxidized (e.g. CO2) ices, as observed in 
comets, and as found by our chemo-dynamical modeling. Fi- 
nally, they have studied in detail chemistries of O-, N-, and C- 
bearing species and found that H3 and He + produce e.g. CO2, 
CH4, and NH3 from CO and N2, and that grain properties 
and ionization rate are crucial factors for molecular evolution. 
This is consistent with our findings, but with an additional 
notion that X-rays further strengthen He + and H3 influence on 
the che mical evolution of turbulen t protoplanetary disks. 

Later, IWoods & Willacvl (120071) have investigated the gen- 
esis of benzene in protoplanetary disks at r < 35 AU. They 
have used a D' Alessio-like flaring a-disk model with accurate 
UV radiative transfer and detailed calculations of the heating 
and cooling balance of gas. The X-ray ionization and disso- 
ciation have not been taken into consideration. The chemical 
network has been ad opted from a sub-set of the UMIST 99 
dLe Teuff et al.l l2000) ratefile augmented with gas-grain and 
surface reactions (in total 200 species and 2400 reactions). 
The passage of a gas par cel from > 35 AU onto the central 
star has been calculated. Woods & Willacv (2007) have found 
that the radial transport results in efficient production of ben- 
zene at < 3 AU, mostly due to ion-molecule reactions between 
C3H3 and C3H4, followed by grain dissociative recombina- 
tion. These results are very sensitive to the adopted value 
of C6H6 binding energy. Similarly, we find that benzene is 
produced via the same ion-molecule reaction because the sur- 
face network is limited for heavy hydrocarbons and does not 
include surface routes to C(,H^, whereas the gas-phase pro- 
duction is efficient. Nevertheless, as in the paper of Woods 
& Willacy, the abundances and column densities in the 
inner disk at 10-30 AU are increased by vertical turbulent dif- 
f usion in our mo d el by up to 2 orders of magnitude, see Fig. [8] 

iNomura et al] d2009), using a disk chemical model with 
radial advection, studied the evolution of the inner region 
(< 12 AU). The a-model of a disk around a 6000 K star 
with the mass accretion rate of 1O~ 8 M0 yr -1 has been consid- 
ered, with gas and dust temperatures computed independently. 
They have modeled penetration of the U V and X-ray radiation 
into the disk, and utilized the RATE 06 dWoodall et al.ll2007l) 
with gas-grain interactions. The chemical evolution of gas 
parcels following various inward trajectories has been calcu- 
lated. The authors have found that the abundances are sensi- 
tive to the transport speed, and that the fast transport allows 
gaseous molecules to reach disk regions where they would 
otherwise be depleted. They have concluded that dynamical 
disk model facilities synthesis of methanol, ammonia, hydro- 
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gen sulfide, acetylene, etc., and predicted the observability of 
the J = 3o-2 methanol line at 145 GHz with ALMA. As we 
have shown above, our 2D-mixing model follows the same 
tendency, and enriches molecular content of the protoplane- 
ta ry disk with complex species. 

Visseretal. (2009b) has developed a 2D semi-analytical 
evolutionary model that follows collapse of a dense core into 
a disk (a transition from the Class to Class I phase) and 
tracks by trajectories infall motions of gas. They have used 
time-dependent RT to compute dust temperature consistently, 
and considered only freeze-out and evaporation of CO and 
water ices. Their major results are that the CO ice evaporates 
during collapse and later re-accretes to dust, while the wa- 
ter ice remains almost intact (unless heated by a protostar at 
r < 10 AU). Material inside 5-30 AU in the early disk spends 
< 10 4 years at temperatures of ~ 20-40 K, which is sufficient 
to produce complex organics. With our disk physical model 
that is at a steady-state and incorporates a detailed gas-surface 
chemistry we indeed see a rise in abundances of many com- 
plex molecules, as their surface production at T < 30 K and 
evaporation are enhanced by the diffusive mixing. 

(Models of protoplanetary disks: the role of turbulent trans- 
port) Now we compare our results to more closely related 
studies of disk chemical evol ution with turbulen t mixing. In 
the first paper on this topic, Ilgne r et ail (120041) have inves- 
tigated the influence of the ID vertical mixing and the ra- 
dial advection on chemistry of a steady 1+lD-a-disk model 
(a = 0.01, M = 10~ 8 M©yr _1 ). The modeling of a reac- 
tive flo w system has be en discussed in detail. They have 
adopted Xie et al. (1995) description for the diffusion and La- 
grangian desc ription for advective transport (trajectories from 
10 to 1 AU). Illgner et all (120041) have used a subset of the 
UMIST95 database with gas-grain interactions and no sur- 
face chemistry (w 240 species, 2 500 reactions). Ilgner et al. 
have considered 3 steady-state solutions for the disk physical 
structure, and investigated the disk evolution by adding ad- 
vective and turbulent transport processes. Their key findings 
are that the chemistry is sensitive to the disk thermal profile, 
that the vertical mixing removes vertical abundance gradients, 
and that local changes in species concentr ations due to mixing 
can be radially transported by advection. Illgner et ail (2004) 
have concluded that diffusion has a limited effect on the disk 
regions dominated by gas-grain kinetics, though it enhances 
abundances of atomic oxygen and thus alters the evolution 
of related species (SO, SO2, CS, etc.). The advective trans- 
port without mixing results in destruction of oxygen and OH 
at r < 5 AU. Our model, despite being limited to outer disk 
regions (beyond 10 AU), confirms their results that oxygen 
abundances are increased by diffusion (Pig |9j, further propa- 
gating into the chemistry of many related species, particularly 
sulfur-bearing species like SO and SO2. We find that turbu- 
lent mixing may steepen or soften chemical gradients (e.g., 
for atomic ions and most of molecules, respectively), and it 
definitely affects more strongly the evolution of species pro- 
duced via gas-surface kinetics. This fundamental feature of 
the presented chemical model is due adopted extended set of 
surface photo- and recombination processes. This further em- 
phasizes the need in the extended surface reaction network. 
Some conclusions, especially those related to complex organ- 
ics do depend on the complexity of the adopted surface chem- 
istry mo del. 

Next, IWillacy et al.1 (120061) have utilized a steady-state a- 
disk model similar to that of Ilgner et al. and used the Xie 



approach to account for the disk viscosity. A plane-parallel 
model of the UV irradiation, no stellar X-ray radiation, and a 
subset of the UMIST 95 d atabase with gas - grain interactions 
and surface kinetics from lHasegawa et al.l (11992) have been 
adopted. Our chemical network is about 3 times larger, and 
uses recently updated reaction rates. They have studied the 
impact of the ID-vertical mixing on the chemical evolution 
of the outer disk (r > 100 AU). Overall, they have found that 
the vertical transport can increase column densities by up to 2 
orders of magnitude. Another result is that the 3-layered disk 
structure is preserved in the mixing model, albeit depths of 
many molecular layers are increased. They have observed that 
the higher the diffusion coefficient, the higher the impact of 
the turbulent mixing. Our results confirm these general find- 
ings, with the chemical network, that is about 3 times larger 
and based on recently updated reaction rates. In contrast to 
the results of Willacy et al., we obtain that the evolution of 
the ionization degree, ammonia and N2H + are sensitive to the 
mixing, whereas column densities of CO, H2CO, CN, C2H are 
steadfast (Tables |5]-[7]i. This is because in our model the X-ray 
irradiation of the disk, production of reactive radicals by pho- 
todissociation of ices, and UV-photodesorption are taken into 
account. In the absence of photoevaporation, e.g., the CO pro- 
duction in the model of Willacy et al. becomes responsive to 
transport of ices from midplane and atomic C and O from at- 
mosphere. In our model this is true only for heavy ices with 
large binding energies (< 3000-5000 K), like hydrocarbons. 
Also, our disk study focuses on the large distances from 800 
till 10 AU that has not been considered by Willacy et al. If 
we restrict our model to the distances beyond 100 AU, the 
ionization degree and column densities of many other species 
become steadfast. Furthermore, our study supports the result 
that the column densities of complex molecules like methanol 
are greatly enhanced by the disk dynamics. 

In 2006 Ilgner & Nelson have investigated in detail ion- 
ization chemistry in disks and its sensitivity to various phys- 
ical and chemical effects at r < 10 AU, like the X-ray 
flares from the y oung T Tauri star (Illgner & Nels on 2006b), 
vertical mixing (Illgner & Nelson! I2006al) . amount of gas- 
phase metals, and vali dity of various chemical networks 
(Illgner & Nelson] I2006d) . Using an a-disk model with stel- 
lar X-ray-irradiation, they have d emonstrated that the simple 
Qppenhei mer & Dalgarnol (119741) network made of 5 species 
tend to overestimate the ionization degree since metals ex- 
change charges with polyatomic species, and that magneti- 
cally decoupled "dead" region may exist in disks unless small 
grains and metals are removed from gas. Next, the influence 
of vertical diffusion on the ionization fraction has been stud- 
ied. They have found that the mixing has no effect on X(e~) 
if metals are absent in the gas since recombination timescales 
are fast, whereas at X(Me) < 10~ 10 - 10~ 8 r c h e m > Tmix and dif- 
fusion drastically reduces the size of the "dead" zone. Finally, 
it has been shown that in the disk model with sporadic X-ray 
flares (by up to a factor of 100 in the luminosity) the outer 
"dead" zone disappear, whereas the inner "dead" zone evolves 
along with variations of the X-ray flux. Indeed, our extended 
disk chemical model shows that ionization degree and abun- 
dances of charged atoms and molecules are sensitive to the 
transport in the disk regions with r < 100-200 AU, and that 
poly atomic ions are impor tant charge carriers in the disks (see 
also ISemenov et aLl l2004). Also, the role of the X-ray pho- 
tons is crucial for the dis k ionization fraction and its molecu- 
lar composition (see also Schreveret al. 20081 [Pfenning et al] 
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l2010h . 

The first attempts to model self-consistently disk chemical, 
physical, a nd turbulent structu res i n full 3D have been per - 
formed by iTurner et all {2006) and lllgner & Nelsonl (12008b . 
Both studies have employed a shearing-box approximation to 
calculate a patch of a 3D MHD disk at radii of ~ 1 -5 AU, 
treated the development of the MRI-driven turbulence, and 
focused on the multi-fluid evolution of the disk ionization 
state. In the study of Turner et al. ohmic resistivity, verti- 
cal stratification, the CRP-ionization, and the Oppenheimer- 
Dalgarno time-dependent ionization network have been con- 
sidered. They have found that turbulent mixing transports 
free charges into the dark disk midplane faster than these re- 
combine, coupling the midplane matter to magnetic fields. 
As a result, accretion stresses in the disk "dead" zone are 
only several times lower than i n surface layers, with a typi- 
cal timescale of about 1-5 Myr. lllgner & Nelsonl (120081) have 
adopted a similar approach but considered the X-ray ioniza- 
tion with Lx = 10 31 ergs" 1 . They hav e re-confirmed their ear- 
lier findings (lllgner & Ne lson 2006a c) that turbulent mixing 
has no effect on the disk ionization structure in the absence 
of the gas-phase metals. The presence of metals, however, 
prolongs the recombination timescale, and the mixing is thus 
able to enliven the "dead" zone at r > 5 AU (with the result- 
ing a = 1 -5 10~ 3 ). Both these studies are based on a pure 
gas-phase chemical network, whereas in our extended gas- 
grain model metals are fully depleted in the disk regions with 
T < 150-200 K. Cons equently, our modeling c orresponds to 
the "no metals" case of Ilgner & Nelson| (120081) . so turbulent 
transport does not much affect the ionization fraction in the 
disk midplane (see Fig. [7] 1st panel). 

Finally, in the recent paper by iHeinzeller et al.l (1201 ll) the 
chemical evolution of the protoplanetary disk along with ra- 
dial viscous accretion, vertical mixing, and vertical wind 
transport has been investigated. The steady-state disk model 
with a = 0.01, M = 1O" 8 M yr" 1 , and L x = 10 30 ergs" 1 has 
been adopted. The heating and cooling processes have been 
included to calculate consistently the gas temperature in the 
disk surface region. In the inner disk region the gas and dust 
temperatures start to diverge when the zjr ratio reaches about 
0.2, which for our disk model translates to about 1 .5H r . Thus 
our approximation of equal dust and gas temperatures is accu- 
rate in the disk midplane and the lower part of the molecular 
layer, and becomes unrealistic in the disk atmosphere. Yet it 
barely affects column densities of many molecules as these 
are dominated by the molecular concentrations in the dense 
re gions where Tjw r yas 

IHeinzeller et al.l (120111) hav e adopted the gas-grain 
RATE 06 dWoodall et al!Tl2007h network, no surface re- 
actions apart from the H2 formation, and the X-ray and 
UV-photochemistry (375 species and about 4350 reactions). 
Heinzeller et al. have concluded that water and hydroxyle 
abundances in the disk surface regions may increase substan- 
tially via neutral-neutral reactions with molecular hydrogen 
produced via chemisorbed-assisted surface recombination. 
It has been found that the disk wind has a negligible effect 
on disk chemistry as its upward transport speed is too low 
compared to the longest chemical timescale attributed to the 
adsorption. Note that in our model the adsorption timescale 
is also long, but the surface chemistry timescales are even 
longer (> 1 Myr). They have pointed out that the radial 
accretion flow alters the molecular abundances in the cold 
midplane, whereas diffusive turbulent mixing affects the 



disk chemistry in the warm molecular layer (r = 1.3 AU, 
T ~ 200 K). We find that the radial turbulent transport can 
affect abundances in the midplane at r > 10 AU (see, e.g., 
CO2 ice in Fig. |9), while effective vertical mixing operates 
from midplane well through the molecular layer and lower 
disk atmosphere. They have concluded that diffusive mixing 
smoothens the chemical gradients and that the abundances 
of NH3, CH3OH, C2H2 and sulfur-containing species are the 
most enhanced. This has been related to increased ammonia 
abundances in the transport model, in which NH3 is effec- 
tively produced via gas-phase reactions with more abundant 
oxygen. In our simulations temperatures are usually well 
below 200 K, but we also see that sulfur-bearing molecules 
are among the most sensitive species to the disk dynamics, 
along with complex organics (like methanol), hydrocarbons 
(like acetylene), and other species (like ammonia). This 
is due to their dependence on slow surface processes and 
evaporation (for heavy gaseous species). 

(Future directions of development for disk chemical mod- 
els) Overall, there is a considerable progress over the last 
decade in constructing more feasible chemo-dynamical mod- 
els of the early Solar nebula and other protoplanetary disks . 
It will gain further momentum in light of the forthcoming 
ground-breaking Atacama Large Millimeter Array (ALMA) 
currently under construction in Chile. In what follows, we 
surmise possible future directions of the development of the 
nebular and protoplanetary disks chemical models. The full 
3D chemo-MHD models of the entire disk evolution within 1- 
5 Myr are beyond possible even with modem computational 
resources available for astronomers. Only local 3D mod- 
els with radically reduced chemistry and a short evolution- 
ary time span of 1 000- 10000 years or semi-analytical mod- 
els of the entire disk with extended chemistry and t ~ 1 Myr 
will be feasible in the foreseeable future. The inner, planet- 
forming disk regions subject to disk-planet(s) interactions 
(gaps, spiral waves, shocks) may well be far from axial sym- 
metry and require 3D modeling approach (e.g ., IWolf et all 
120021 iFukagawaetail 120041 iPiitu et all 120051 The accu- 
rate opacities of dust and gas covering a wide range of tem- 
peratures, densities, gas composition, and dust topological 
and mineralogical properties will have to be developed (e.g., 
iSemenov et al]l2003tlHelling & Lucasll2009l) . 

Likely, the steady-state disk physical model adopted in 
the present study needs to be re-adjusted into an evolution- 
ary model with consistent calculations of grain evolution 
(e.g. iBirnstiel et aTll2010t iFogel et al.ll201 lUVasvunin et alJ 
1201 ll) . photoevaporat ion ("e.g..lGorti et al.ll2009l) , and accurate 
gas te mperature (e.g.. lGorti & H ollenbach 2004; Woit ke et alJ 
2009). The grain coagulation, fragmentation, and sedimen- 
tation lead to re-distribution of the total dust surface area, 
higher UV ionization rates, and thus to shift of molecular lay- 
ers toward the midplane, with a significant increase of molec- 
ular column densities in t he intermediate layer and decrease 
of th eir depletion zones (Fo gel et aD 1201 It IVasyunin et al.1 
1201 ll) . In the evolutionary disk models the total disk mass 
and size, and thus its density, thermal, and ionization struc- 
tures change with time, altering conditions at which chem- 
ical processe s proceed both in the gas phase and onto the 
dust surfaces dGorti et al.ll2009l) . In evolutionary models cov- 
ering formation and build-up disk phases, a large amount of 
gas and dust materials experience events of heating, cooling, 
re-co ndensation, annealing, and varying irr adiation intensities 
(e.g.. lG^l200ll2004t Ivisser et al.ll2009bl) . 

A more realistic 2D or full 3D prescription of the X-ray 
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and UV radiation transfer modeling with scattering has to 
be included. The accurately calculated UV spectrum with 
the L a line is vital to compute realistic photodissociation 
and photoioniza tion rates, and shielding factors for CO and 
H? (see, e.g.. Ivan Zadelhoff et all 12003? iBethell et afl 120071: 
IVisser et aLll2009ah . The realistic (may be, variable) stel- 
lar X-ray spectrum is required to calculate gas temperature 
in the disk atmosphere, and for modeling disk ionization 
structure and ion-molecule chemistry (e.g.. iGlassgold et alj 
12005b iMeijerink et al.ll2008b lAresu et alJl201 ll) . At densities 
of < 10 4 — 10 6 cm" 3 gas temperature decouples from that of 
dust, which will affect the chemi cal and dynamical evolu- 
tion in the upper disk layers (e.g., Kamp & Dullemond1 l2004t 
lOwen et al. 2011). The dissipation of Alfvenic waves gener- 
ated by the MHD processes i nside the disk can also heat gas 
in the disk atmospheres (e.g., Hirose & Turner 2011b. 

In comparison to the status of the disk physical structure 
where many key issues have been realized and partly solved, 
the advances in modeling disk chemistry are less straight- 
forward. First of all, chemical studies employ various net- 
works fully or partly based on astrochemical ratefiles that 
have been developed to model distinct astrophysical environ- 
ments. Hopefully, this will be relieved with the advent of the 
public Kinetic Database for Ast rochemistry (KIDAfl that in- 
cludes state-of-the-art rate data (Wakel am et al.ll2010bl) . Un- 
fortunately, the intrinsic uncertainties in gas-phase reaction 
rates hamper the accuracy of chemical pred ictions both on 
abundances and column densities (see, e.g.. IWakelam et alj 
l2006tlVasyunin et al.ll200l IWakelam et alJl2010al) . The pace 
at which these quantities are measured in laboratories or 
calculated by quantum chemical models precludes the rapid 
progress for thou sands of astrochemically-relevant processes 
(iSavin et al.ll20TTh . 

In the inner (< 5 AU) disk regions with T > 1 00- 300 K and 
densities exceeding 10 12 cm" 3 three-body processes, many 
reverse reactions, and neutral-neutral reactions with large 
barriers are activated, and have to be fully taken into ac- 
count (e.g. lAikawaetaLl 12002b iTscharnuter & Gaill 120071; 
lHarada et al.ll2010l) . The X-ray-driven ionization and disso- 
ciation of molecules other than H2 are poorly understood , 
but important for disk chemistry (e.g., Glassgold et al. 2009). 
In the upper disk regions, upon gas-phase or surface recom- 
bination or due to ionization/dissociation an excess of en- 
ergy may translate into (ro-)vibrational excitation of a product 
molecule, which may rea ct differently with other species (e.g., 
iPierce & A'Hear n 2010). This aspect has so far been almost 
completely neglected in astrochemical models. 

In cool outer disk regions gas-grain interactions and sur- 
face reactions are essential, yet the latter are often disre- 
garded in disk models. The same is true for nuclear-spin- 
dependent chemical reactions involving ortho- a nd para-states 
of key species, such as Ht, H2, H 2O, etc. (e.g., Paga ni et alj 
1 19921120091: iCrabtree etatll2011l) . Usually weak physisorp- 
tion of molecules is considered in chemical models, while 
molecules can also form chemical bonds with dust surfaces, 
leading to heterogeneous surface chemistry active both in 
cold and hot regio ns, particula r ly on PAHs and carbona- 
ceous grains (e.g., iFraser et al.l 12005b iCazaux et al.l 120051 
Cuppen & Hornekaer 2008}. Despite the recent laboratory ef- 
forts to measure binding energies of key molecules like CO, 
N2, water, etc. to various astrophysical ices and their pho- 

6 http://kida.obs.u-bordeauxl.fr 



todesorption yields, many of these values are still lacking ac- 
curate estimates dBisschop et al.ll2006b lOberg et al.ll2009bllah . 
The dynamics, reactivity, photodissociation, and desorption 
of ices embedded into dust mantles are hard to measure or 
model, a nd even harder to interpret (see, e.g., results for the 
water icejAndersson et al. 2006t lAnders son & v an Dishoeckl 
12008b [Bouwrna net alj|201 ll) . For example, sub-surface diffu- 
sion may increase desorption efficiency in the case of well- 
mixed ices within the H2O ice matrix. Along with ther- 
mal, UV-, CRP- and X-r a y-triggered desorp tion (Leg er etail 
[19851 iNaiita et alj l200lt iWalshetai] 12010). typically con- 
sidered in disk chemical studies, other non-thermal mecha- 
nisms, like grain-grain destruc tion and explosive desorption 
(Shalabiea & Greenberg 1994) can be operative. The sur- 
face recombination on porous grain surfaces or within het- 
erogeneous ices is also non-trivial to model accurately (e.g . , 
iCuppen et ail 12009b iFuchs et all 120091 lloppolo et all I261 lb . 
The surface chemistry can be restricted to several upper- 
most monolayers of the grain, while often in chemical mod- 
els the corresponding rates are calculate d assuming it is ac- 
tive e verywhere in the ice mantle (e.g., Hasegawa & Herbst 
1 19931) . In addition, the presence of large, photostable PAHs 
in the gas assists synthesis of polyatomic molecules by pro- 
viding surface area for surface recombination, ability to re- 
radiate energy releazed upon non-destructive ion-PAH recom- 
bination, and du e to the combustion chemistry in the inner 
disk region (e.g.. IWakelam & Herbstl l2008l iKress e taT1l2010b 
Perez-Bec ker & Chiangll201 lb . Finally, calculated molecular 
concentrations depend on the adopted set of elemental abun- 
dances that may vary from region to region, and underlying 
assumption s about depletion of he avy elements from the gas 
phase (e.g., Wakelam et al. 2010a). 

5.2. Cold molecules in DM Tau 

In the following two subsections we elaborate on feasibil- 
ity of the present study to explain observations of diagnostic 
molecules and ices in protoplanetary disks and the Solar sys- 
tem. 

The kinetic temperature distributions have been measured 
in near by disks of DM Tau LkCa 15, MWC 4 8 0, and AB 
Aur bv IDartois et al l (120031) : IPietu et al.l (120051 120071) and 
iHenning et alj d2010l) . using the Plateau de Bure Interferom- 
eter and the (1-0) and (2-1) emission lines of CO isotopo- 
logues, HCO + , and C2H These emission lines have different 
opacities and thus sample gas temp erature at various h eights 
above the disk equatorial plane (see IDartois et al.|[2003l) . CO 
lines are particularly suitable for such studies because they 
are easily excited at low densities of 10 4 cm" 3 , while the 
CO isotopologue abundances vary by orders of magnitude. 
Consistentl y with predictions from a passively heated, flaring 
model (e.gJChiang & GoldreichliggTtlD'Alessio et al.ll999b 
iDullemond & Dominikll2004l) . vertical temperature gradients 
have been found. In disks around hot stars midplane tem- 
peratures are quite high. At ~ 100 AU, the disk around the 
10000 K Herbig AO AB Aur has gas temperatures ranging 
from 20 ± 3 K (as measured by C 18 7=1-0) t o 35-40 K 
( 13 CO 1-0, 2-1) and to 68 K ( 12 CO 2-1), see IPietu etall 
(2005). A similar trend has been obtained for a disk sur- 
rounding the cooler, 8500 K Herbig A4 MWC 480 star, with 
CO temperatures growing from 21 ±4 K in the midplane 
to 48 ± 1 K in th e molecular layer at 100 AU from the star 
(IPietu et alj|2007l) . 

In the bright, large 800 AU disk around cold, 3750 K 
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(Ml) DM Tau star very cold CO temperatures of pa 15 K at 
100 AU have been reported bv lDartois et al.l (12003b . based on 
analysis of the PdBI 13 CO (1-0) spectral map. In the later 
study by lPietu et alJ (12007b even lower temperatures of 8 and 
15 K probed by the (1-0) and (2-1) transitions of I3 CO, and 
r kin = 14 ± 2 K probed by the HCO + (1-0) line have been de- 
rived for the outer DM Tau disk. In contrast, in the same disk 
optica lly thick 12 CO lines give kinetic temperature of about 
25 K. iHenning et alJ ( 120101) found very low kinetic temper- 
ature of rs 7 K for C2H in the DM Tau system, using the 
CCH (1-0) and (2-1) PdBI data. Finally, in Chapillon et 
al. (2011, submitted) similarly low values for the CN (1-0) 
temperatures have been inferred. While derived low temper- 
atures are characteristic for midplane regions of protoplan- 
etary disks surrounding cool T Tauri stars, the presence of 
gaseous molecules at such extreme conditions is puzzling. 
The evaporation temperatures for CO, C2H, and CN are about 
20, 30, and 40 K, respectively (see Table [TJ, and in dense, 
dark midplane (T ~ 10-20 K) they are rapidly depleted, 
T acr < 10 3 - 1 4 years (Tables |3US. 

In Paper I (ISemenov et al.l 12006) the problem of the cold 
CO gas res ervoir in the DM Ta u disk has been tackled using 
the 1+lD lD'Alessio et all d 19991) a = 0.01 flaring disk model 
with the 2D radial- vertical turbulent diffusion. We have found 
that the CO column density of about 10 16 cm" 2 is required in 
the disk midplane to explain the puzzling observations of the 
cold CO in DM Tau. Such CO column densities in the disk 
region with T < 25 K have been obtained in the model with 
vertical and radial mixing at all radii, while in the laminar 
model these values only appeared at r > 500 AU. The CO 
molecules are transported from the warm molecular layer to 
the cold disk midplane, maintaining concentration of CO even 
at late times, t > 10 - 10 6 years. 

Aikawa (2007) have used a steady-state a-disk structure 
and a ID-vertical mixing model to show that the warm CO 
gas from the intermediate layer can be transported down to 
the cold midplane at a rate that can be competitive with the 



CO adsorption rate. The efficiency of this process has been 
found particularly pronounced for the disk model when mod- 
erate grain growth is allowed, with grain sizes < 1/im. A sub- 
stantial grain growth leads to a decrease of the averaged dust 
surface area and thus less active freeze- out (see also results of 
iFogel et alJl20TTtlVasyunin et alj201 ll) . 

Later, iHersant et al.l (|2009) have utilized a power-law, a- 
disk model with an isothermal midplane, strong UV scatter- 
ing, and ID-vertical mixing. They have found that the UV- 
desorption of CO and other ices prevails over their upward 
transport by the vertical mixing to warm regions, leading to 
a large amount of gaseous CO, HCO + , HCN in the midplane 
even in the laminar model. This is due to the adopted high 
UV-desorption yield of - 0. 1 % (e.g. JOberg et al.ll2009bl) and 
assumed high UV penetration efficiency caused by the UV 
scattering on small dust grains in the disk atmosphere. 

Compared to Paper I, we now use the full chemical net- 
work based on the OSU model, higher photodesorption yields 
of 10~ 3 , and the slow surface hy drogenation rates based 
on the results of iKatz et alJ (Il999h . with 5% of products 
to be released back to the gas upon surface recombination 
dGarrod et al.H2007T) . In Figure [l3]we show the column den- 
sities of CO, HCO + , CN, HCN, C 2 H, and CS at 5 Myr in the 
disk regions with T < 20 K (odd rows) and T > 20 K (even 
rows) calculated with the laminar, the slow 2D-, and the fast 
2D-mixing models. Note that in the inner region at r < 30 AU 
kinetic temperatures are higher than 20 K even in the mid- 
plane (Fig.Q] 1st panel). 

As can be clearly seen, turbulent transport does enhance 
molecular abundances and column densities in the DM Tau 
disk midplane for the all species apart from CS. Compared 
to the results of Paper I, the CO column densities in the mid- 
plane region {T < 20 K) are increased by the diffusive mix- 
ing by 2 orders of magnitude, to values between 10 15 and 
10 16 cm" 2 , while in the laminar model these are still rela- 
tively high, > 10 14 cm -2 . As in the study of Hersa ntet alJ 
(2009) powerful CRP-induced and UV-desorption keeps some 
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molecules in the gas phase even in the midplane. Transport 
from the disk regions with T > 20 K, like the molecular layer 
and the inner, accretion-heated midplane at r < 30 AU, where 
CO can be produced on grains, replenishes gaseous CO in 
the cold midplane. The column densities of "warm" CO are 
much higher, sa 10 18 cm -2 and are not affected by the disk 
dynamics. This is because chemical evolution of CO is rela- 
tively fast compared to the transport speeds, while CO is also 
widespread through the disk, canceling o ut th e chemical gra- 
dient needed for its diffusion, see Section l4~3l The amount of 
cold CO in the mixing model is hardly sufficient to explain the 
puzzling 13 CO PdBI observations of DM Tau. On the other 
hand, there is enough warm CO to explain the presence of the 
T f» 26 K CO gas derived from the optically thick 12 CO data 
tracing the disk upper region. 

HCO + shows an increase of midplane column densities by 
up to a factor of 100 in the mixing case, while in the molecular 
layer its column densities are enhanced by a factor of < 30. 
This is due to the sensitivity of H3 ion, required for HCO + 
pro ducti on, to the turbulent transport (see discussion in Sec- 
tion l4.U . In the mixing model, H3 is particularly abundant at 
elevated disk heights, z > 1 —2H r , and so is HCO + . The ra- 
tio of column densities of warm and cold HCO + is large, and 
nearly the same for both the laminar and 2D-mixing models, 
a factor of ps 30 and 10- 100, respectively. With a typical 
column density of 10 13 - 10 14 cm" 2 at r > 100 AU, the HCO + 
(1-0) line remains optically thin in the outer disk, sampling 
both the disk midplane and the warm molecular layer. 

On the other hand, the critical density required to excite a 
molecular line from a linear molecule can be estimated as: 

««*=—, (19) 
av 

where A u \ is the Einstein coefficient for spontaneous emis- 
sion from the level u to /, and a and v are the cross- 
section and velocity of a collisional partner, respectively. 
For the H2-dominated disk midplane with T = 10-20 K 
av rs 3 10" 11 em's" 1 (e.g.. IWilson et al.l 120091). We adopt 
A 10 = 4.25 10" 5 s" 1 from the LAMDA database (ISchoier et al.l 
12005b . and find that the critical density for the HCO + 7=1-0 
excitation is « 1.5 10 6 cm" 3 . The region where densities ex- 
ceed 10 6 cm" 3 extends vertically up to < 3 and 1 H r in the in- 
ner and outer disk, respectively, and covers the midplane and 
the molecular layer (see Fig. [TJ. Thus, HCO + (1-0) emission 
would originate from the T > 20 K disk region, and the re- 
sulting temperature derived from this transition would likely 
exceed 15 K obt ained from interferometric observations by 
iPietu et all (f2007b . 

The behavior of C2H, CN, and HCN is sim ilar as their 
chemical evolution is tightly linked (Section |4.4| >. Their ver- 
tical column densities are not much altered by the turbulent 
mixing, having a typical value of ~ 310 I3 -10 14 cm" 2 at 
100 AU. Since absolute concentrations of C 2 H, CN, and HCN 
are highest in the molecular layer where chemical timescales 
are relatively short, their T > 20 K column densities are 
not strongly increased by diffusion, less than by a factor of 
3. Contrary, the column densities in the DM Tau midplane 
(T < 20 K) are raised by the transport by an order of magni- 
tude for CN and HCN, and by 2 orders of magnitude for C2H. 
This is caused by the long evolutionary timescales in the mid- 
plane associated with gas-surface kinetics (e.g., gaseous HCN 
is produced in the midplane upon recombination of the H and 
CN ices). The resulting column densities through the mid- 



plane are > 10 13 cm" 2 . Thus, in the fast 2D-mixing model 
column densities of warm and cold C2H, CN, and HCN differ 
by only a factor of several, which further decreases outward. 
Using expression (fT9l and the LAMDA database, we estimate 
that the optically thin 1-0 rotational lines of C2H, CN, and 
HCN are excited at densities of about 4-8 10 5 cm" 3 . Conse- 
quently, in the presence of strong turbulent mixing low-lying 
transitions of C2H, CN, and HCN likely trace the cold mid- 
plane of DM Tau, supporting the observational evidence. 

Finally, turbulent transport decreases CS column densities 
in the midplane region by a factor of < 4, simultaneously in- 
creasing it in the T > 20 K zone by up to 2 orders of magni- 
tude. The chemistry of sulfur-bearing molecules is among the 
most altered by the mixing due to slow surfa ce pr ocesses as- 
sociated with these heavy species (see Section l4~5l l. However, 
even in the laminar model column densities of cold (T < 20K) 
CS are lower than that of warm (T > 20K) CS by an order of 
magnitude. The overall CS column densities in the outer disk 
are » 3 10 11 -3 10 12 cm" 2 , and the low-/ CS lines are opti- 
cally thin. The critical densities of the excitation for the (2-1) 
and (3-2) CS transitions at 98 and 147 GHz are - 5 10 6 and 
210 6 cm" 3 , respectively. Therefore, if our model is correct, 
the kinetic temperatures derived from the low-lying CS emis- 
sion lines should be above 20 K in the DM Tau disk. 

Clearly, our modeling shows a potential of turbulent trans- 
port as a cause for the presence of molecules in the cold mid- 
planes in T Tauri disks. The non-thermal broadening of the 
CO lines b y < 150 ms" 1 in the DM Tau system has been 
reported in Piet u et al] (|2007), which corresponds to the tur- 
bulent velocities of < 10% of the sound speed. However, 
the increase in concentrations of the cold gases in our 2D- 
mixing model is not strong enough to rule out other explana- 
tions. First of all, more accurate modeling of the UV scat- 
tering toward the midplane would result in faster photoevap- 
oratio n of ices in the lower molecular layer and upp er mid- 
plane (Ivan Zadelhoff et al l 120031; iHersant et alj|2009h . Sec- 
ond, as mentioned in lAikawal d2007ir moderate grain growth 
beyond 1 [im in disk central regions lengthens the depletion 
time for molecules such that the gas-phase CO abundances 
may remain high even afte r 1 Myr of the evo lution. Indeed, 
IVasyunin et alJ (l201 lb and iFogel et all (1201 lh have modeled 
chemical evolution in disks taking grain evolution into ac- 
count and found that the grain growth substantially decreases 
the depletion zones of molecules. Moreover, from analy- 
sis of the SED slopes in millimeter and centimeter wave- 
lengths large grain sizes of 1 mm have been inferred for 
many young systems in various star-f orming regions, includ- 
ing Taurus-Auriga association (e.g.. [ Rodmann et al.l 120061; 
iLommen etai]l2010t iRicci et al. 2010). The DM Tau disk 
has an inner hole of ~ 4 - 20 AU with grains as large as 
1 mm, and hence is in a tr ansitional phase (e.g., ICalvet et al.1 
12001 ISargent et all 120091; lAndrews et ail UM]\> . From ad- 
vanced theoretical models of grain coagulation, fragmenta- 
tion, sedimentation, and turbulent stirring the grain evolution 
should proceed everywhere in the disk, though the mean grain 
sizes remain smalle r in the outer disk compared to the in - 
ner disk r egion (e.g.JBrauer et al. 2008; Birn stiel et al.ll20 lOh . 
Recently, Guilloteau et alJ (1201 lh have used high-resolution, 
multi-frequency interferometric PdBI observations to discern 
more accurately the dust emissivity slopes at millimeter wave- 
lengths in a sample of young stars. Their analysis has indeed 
shown that the outer disk of DM Tau contains large grains 
with sizes > 1 /im. We advocate for a combined action of 
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turbulent transport and grain growth as a mechanism to main- 12 CO(3-2) 
tain a sizable reservoir of cold gases in the DM Tau system. 

5.3. Comparison with observations 



Table 14 

Observed and modeled column densities in DM Tau at 250 AU 



Species 


Observed 


Refs 


Modeled 

Laminar 2D-mixing (Sc = 100) 


2D-mixing (Sc =1) 


12 CO* 


3.0(17) 


(1) 


3.0(17) 


3.0(17) 


3.0(17) 


HCO+ 


1.7 (13) 


(1,2) 


8.0(12) 


1.2(13) 


2.5 (13) 


H 2 CO 


1.0-2.0(13) 


(3) 


6.2 (12) 


6.7(12) 


4.0(12) 


N 2 H+ 


4.0(11) 


(2) 


3.4(11) 


4.7(11) 


1.2(12) 


CS 


4.0 (12) 


(3) 


8.1 (10) 


1.1(11) 


5.6(11) 


CN 


4.0(13) 


(4) 


1.4(13) 


1.5(13) 


1.8(13) 


HCN 


8.0(12) 


(4,5) 


1.2(13) 


1.2(13) 


2.7(13) 


HNC 


3.0(12) 


(4,5) 


1.0(13) 


1.0(13) 


2.4(13) 


C 2 H 


3.0(13) 


(6) 


1.1 (13) 


1.3 (13) 


1.4(13) 


Agreement** 






7/8 


7/8 


6/8 



References. — (1) Pietu et al. (2007); (2) Dutrey et al. (2007); (3) Dutrey 
et al. (201 1, submitted); (4) Chapillon et al. (201 1, submitted); (5) Schreyer 
et al. (2008); (6) Henning et al. (2010) 

* 

The calculated column densities are scaled down to match the observed 
values for l2 CO. The renormalization factor is applied to the column densi- 
ties of other species. 

** 

Agreement is achieved when the observed and calculated column densi- 
ties for a molecule do not differ by more than a factor of 4). 

Table 15 

Observed and modeled column densities of ices in DM Tau at r = 10 — 30 AU 



Species 


Relative 
abundance 


Ref 


10 AU 


20 AU 


30 AU 


H 2 


100 




100 


100 


100 


CO 


5-30 


(1,2,3,4) 


<0.01 


0.7 


16 


co 2 


3-20 


(1,4) 


10 


20 


0.6 


H 2 CO 


0.04-1 


(1,2,3,4) 


0.02 


0.015 


<0.01 


CH 3 OH 


2 


(2,3,4) 





<0.01 


<0.01 


HCOOH 


0.06-0.09 


(3,4) 


0.02 








HCOOCHt 


0.06-0.08 


(3,4) 











CH 4 


0.6-5 


(1,3,4) 





40 


30 


C 2 H 2 


0.1-0.5 


(3,4) 











C 2 H 6 


0.3-0.4 


(3,4) 


0.07 








NH 3 


0.1-2 


(1,3,4) 


33 


20 


15 


HCN 


< 0.25 


(1,2,3,4) 


12 


3.6 


0.5 


HNC 


0.01-0.04 


(4,4) 


0.5 


0.02 





HNCO 


0.06-0.1 


(3,4) 











N 2 


~0.02 


(1) 











CH 3 CN 


0.01-0.02 


(3,4) 











H 3 CN 


0.02 


(3,4) 


0.003 








NH 2 CHO 


0.01-0.02 


(3,4) 











H 2 S 


0.75-1.5 


(1,2,3,4) 


0.15 


0.09 


0.06 


CS 


~0.1 


(2) 











SO 


0.2-0.8 


(1,3,4) 











so. 


0.1-0.2 


(1,3,4) 











ocs 


0.1-0.4 


(1,3,4) 











H 2 CS 


0.02 


(3,4) 


0.006 








s 2 


0.005 


(4) 












References. — (1) Aikawa et al. (1999); (2) Biver et al. (1999); (3) 
Bockelee-Morvan et al. (2000); (4) Crovisier & Bockelee-Morvan (1999) 

(Molecules in protoplanetary disks) Some previous stud- 
ies of the laminar disk chemistry have reported reason- 
able quantitative agr eement with observationally-inferred col- 
umn densities (e.g.jAikawa et al.ll2 002; Sem enov et al.ll2005t 
iDutrev et al.ll2007bt ISchrever et al.ll2008l) . On the other hand, 
the radial profiles of the column densities derived from fitting 



high-resolution interferometric data have not been fully repro- 
duced by co nventional disk mode ls. In the ID-vertical mix- 
ing study of I Willacy et al.l (120061) the modeling results have 
been compared with single-dish and interferometric observa- 
tions of several molecules in the disks around DM Tau, LkCa 
15, and TW Hya. The good agreement between theoretical 
and observed column densities for the fast mixing model has 
been inferred from the single-dish data, whereas the interfer- 
ometric data have not been reproduced. Calculating infrared 
emission lines for the dis k model with radial a dvective and 
vertical mixing transport. iHeinzeller et al.l (1201 lb have shown 
that it improves agreement with the Spitzer observations of the 
inner disks around A A Tau, DR Tau, and AS 205, compared 
to predictions of the laminar chemical model. 

We compile a table with column densities at 250 AU de- 
rived from the analysis of high-quality PdBI interferometric 
observations and compare them with the results of our lami- 
nar, slow, and fast 2D-mixing models (see Table IT4Y Calcu- 
lated column densities have intrinsic uncertai nties of a factor 
~ 3- 5, caused by reaction rate uncertainties dVasvunin et al.l 
2008), whereas observational data suffer from calibration in- 
accuracies (~ 10-20%), distance uncertainties (10-20%), 
etc. Therefore, we assume that the agreement is good when 
observed and modeled values differ by a factor of < 4 and 
bad otherwise. The disk total surface density cannot be accu- 
rately derived from the continuum data, as these suffer from 
poorly known dust opacities at (su b-)millimeter wavel engths 
that can vary by factors of several (Semenov et al. 2003). In- 
stead, the observed CO column density distribution can be 
used as a proxy to the disk gas density structure. The problem 
here is that the 12 CO rotational lines are optically thick and 
thus probe localized disk region, whereas from optically thin 
CO isotopologue lines the total CO column densities can only 
be recovered if 12 C/ 13 C or 16 0/ 18 isotope ratios are known. 
We use 13 CO and 12 CO column densities obtained for the 
outer DM Tau dis k from the interferometric observations by 
iPietu et ail d2007l) and constrain the total CO column density 
at 250 AU to ss 3 10 17 cm" 2 . The corresponding theoretical 
value is almost the same in the laminar and mixing models, 
N(CO) « 10 18 cm -2 (Fig |9ll. To match this value, modeled 
column d ensities in T abled?] were renormalized accordingly. 

Unlike Willac vet alj (120061) . all the considered models 
match quite well the observed column densities, apart from 
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the CS data. The stellar and disk parameters are well known 
for the nearby DM Tau system, so its physical structure can 
be reliably reconstructed. The detected species are simple and 
their chemical evolution is governed by a limited amount of 
mostly gas-phase processes, with many accurately acquired 
rate constants. However, CS column densities are underpre- 
dicted by factors of ~ 7 - 50, though situation considerably 
improves for the fast 2D-mixing model. Dutrey et al. (2011, 
submitted) have observed SO, H 2 S, and CS in DM Tau, LkCa 
15, and MWC 480 with the IRAM 30-m and PdBI interfer- 
ometer. They have derived upper limits for H2S and SO in the 
DM Tau disk at 300 AU, namely, JV(H 2 S) < 2 10 11 cm" 2 and 
N(SO) < 8 10 11 cm" 2 , whereas CS has been firmly detected 
with the column density of JV(CS) < 3.5 10 12 cm" 2 . Dutrey et 
al. have found that the sulfur content of the three disks cannot 
be explained by modern chemical models as these result in 
higher abundances of H2S and SO compared to the CS abun- 
dances. This is also the case in our modeling. While having 
too less CS, our disk models overproduce column densities of 
sulfur monoxide and hydrogen sulfide. T he re ason for such a 
disagreement, as we discussed in Section |431 is that the sul- 
fur chemistry is hampered by poorly known reaction data and 
may lack key reactions, and the depletion of elemental sulfur 
from the gas is also crudely constrained. 

Generally, the fast 2D-mixing model produces more 
molecules per CO compared to the laminar disk chemistry, 
and overproduces HNC with respect to the observations. The 
utilized chemical network leads to the synthesis of hydro- 
gen cyanide and isocyanide with almost equal probability in 
the all disk models. On the other hand, the observed col- 
umn density of HCN is higher than that of HNC by a fac- 
tor of - 2 - 3. The HCN/HNC ratios of < 1 are typical 
for da rk, dense cores dTenneke s et al. 2 0061 iHily-Blant et al.l 
2010), while in irradia ted environments like PDR regions this 
ra tio exceeds 1 (e.g.jLoenen et al.l l2008). Indeed, in the study 
of lFogel et al.l (1201 ll) with accurate description of the UV RT 
including the L Q radiation and the moderate grain growth the 
resulting ratio of HCN/HNC is about 15. The HCN/H NC ra- 
tio of 3 has been obtained by IVasyunin et al.l (1201 ll) for the 
disk model with consistently calculated grain evolution and 
the UV penetration. An increase in the HCN/HNC ratio to 
3-5 due to lowering the disk mass or including d ust se dimen- 
tation has been reported by Aikawa &Herbst (1999). Fur- 
thermore, ISarrasin et al.l (120101) have accurately calculate d the 
collisional data for HCN and HNC with He and found that 
the use of the HCN rates to interpret HNC observations in 
dark clouds may lead to the underestimation of the HCN/HNC 
ratio. Apparently, the use of ISM-like dust grains or the 
crude approach to calculate UV penetration utilized in our 
disk model make the UV opacities too high in the molecular 
layer to result in HCN/HNC > 1 . 

Column densities or upper limits have also been re- 
ported for o t her sp ecies. Using the IRAM 30-m antenna, 
iDutrev et al.l (119971) have not been able to detect molecular 
lines of SiO, SiS, HC 3 N, C 3 H 2 , CH 3 OH, CO + , S0 2 , HNCS, 
HCOOCH3 in DM Tau. Among these species, only SiO, 
C3H2, and SO2 have large modeled column densities between 
10 12 and 10° cm" 2 at 250 AU. The S0 2 lacks permanent 
dipole moment and has no pure rotational spectra, unless oxy- 
gen atoms are isotopically different. Its ro-vibrational lines 
are excited at T > 50- 100 K, so that the emission arises 
only from an inner disk region (r < 5 - 10 A U). The C3H7 
molecule can be in ortho- or para-state (e.g., [Madden et al.l 



[19891: iMorisawa et aIll2006T) . though in our model these de- 
tails are not taken into account. The linear isomer H 2 CCC 
(propadienylidene) has a large dipole moment of about 4 D 
and rotational spectrum starting at cm wavelengths. However, 
the IRAM 30-m observ ations of TMC1 and IRC +10216 by 
iCernicharo et al.l (119911) have revealed that the abundance of 
the linear isomer is only ~ 1% of that of the cyclopropenyli- 
dene. The cyclic isomer has a rich ro-vibrational spectrum, 
with emission lines at (sub-)millimeter wavelengths that are 
excited at n„ ~ 10 6 cm" 3 (LAMDA database; iGreen et aTl 
119871: ISchoier et aLi r2005). Due to energy partition the indi- 
vidual lines are not as strong as in the case of HCO + , making 
them hard to detect. Finally, the outflow tracer, SiO, has a 
dipole moment of 3.1 D and strong rotational lines at (sub- 
)mm, excited at ~ 10 5 cm" 3 (Schoieret al. 2005). The calcu- 
lated column density of SiO at 250 AU is 4 10 12 cm" 2 in the 
laminar model and 10 13 cm" 2 in the fast 2D-mixing model. 
These values are sensitive to the abundance of the elemen- 
tal silicon remained in the gas and may not be representative 
of DM Tau in th e "low metals" set of initial abundances of 
iLee et all (Il998al) adopted in our modeling. 

The recent non-detection/tentative detection of cold water 
vapor in th e DM Tau disk by t he Herschel/HIFI have been re- 
ported by Bergi n et"aT1 ((2010). They have interpreted the ob- 
servational data with an advanced chemical disk and the line 
radiative transfer model, and inferred the disk-averaged wa- 
ter column densities within ~ 5 10 12 - 3 10 13 cm" 2 . Our cal- 
culated disk-averaged H 2 column densities are 5 10 14 cm" 2 
and 210 I5 cm" 2 in the laminar and fast mixing models, re- 
spectively. Therefore, all our models overestimate the abun- 
dance of gaseous water in DM Tau by a factor of at least 15. 
Similarly large values have been obta ined in the studies of 
disk chemistry with gra in evolution by Vasv unin et al.l (1201 ll) 
and iFogel et all (120111) . Bergin et al. have concluded that 
most water ice is trapped in large dust grain aggregates that 
sedimented toward the DM Tau midplane and thus cannot 
be easily photodesorbed. However, if this is the case, other 
molecules like CO, CN, etc. have to be depleted from the 
gas phase much more severely, and should not be present in 
the cold midplane. The presence of large reservoir of cold 
CO, CN, HCN, C 2 H observed in the DM Tau midplane con- 
tradicts with such hypothesis, unless other molecules begin 
accreting onto dust grains later than water, when grains are 
already large, > 1 /im. 

(Molecules in comets of the Solar system) We cal- 
culate column densities of ices at 10, 20, and 30 AU 
and compar e them with molecu les observed in comets 
as done in lAikawa et all (1 1 9991) . see Table [131 Ob- 
served va lues are ta ken from | Aikawa et all (1999); 
l Biveretal.1 dl999t>; iBockelee-Morvan et al.l (I2000b : 
Crovisier & Bockelee-Morvan (119991 and references therein). 
The difference in calculated ice column densities between the 
laminar and mixin g m odels is smaller than their radial vari- 
ations, so in Table TT31 only results of the laminar model are 
presented. The comets have formed in the early Solar nebula 
at distance s of > 10-20 AU and later have been expelled 
outwards (lLissauer| [T987). Comets are reservoirs of pristine 
compounds, and their chemical composition is indicative of 
evolutionary history of the outer Solar system during the 
first several million years (e.g., lEhrenfreund & Charnleyl 
2000). However, molecules detected in their comae can be 
partly photodissociated or are products of immediate coma 
chemistry and very long irradiation of bulk cometary ices 
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by cosmic ray particles, so direct comparison with model 
predictions is not straightforward. The presence of both 
reduced and oxidize d ices in comets have been explained by 
lAikawa et alJ {T999) as due to CRP-ionization of the disk, 
whereas in cosmochemical models the Fischer- Tropsch catal- 
ysis o r cloud-nebula evolutionary scenario hav e been invoked 
(e.g.. lGreenberglll982t IPrinn & Feglevlll989l) . We find that 
our disk model incorporating the high-energy radiation is 
able to reproduce relatively high abundances of water, CO, 
and CO2 ices, as well as abundances of CH4 and C2H6 ices. 
Also, abundances of H 2 CO, HCOOH, HCN, HNC, H 2 S, and 
H2CS are reproduced. However, the model overproduces 
NH3 ice by a factor of > 7, and severely underpredicts 
abundances of most of the sulfur-bearing ices (CS, SO, 
S0 2 , OCS, S 2 ) and complex organics (CH 3 OH, HCOOCH3, 
HNCO, NH 2 CHO). The chemical model with fast turbulent 
transport gives better results for heavy sulfur-bearing and 
organic molecules, but is still far below the observed values. 
As we discussed above, the sulfur chemistry is rather dubious 
in modern astrochemical databases. The adopted chemical 
network has a limited number of surface reactions leading 
to complex organic molecules and not many endothermic 
neutral-neutral reactions. Calculated ice abundances are also 
sensitive to the adopted binding energies, many of which are 
not accurately derived. 

Apart from that, we conclude, that our disk chemical mod- 
els with and without transport processes are equally and rea- 
sonably well agree to the high-quality interferometric obser- 
vations of DM Tau and the chemical composition of comets 
in the Solar system. 

5.4. Observable molecular tracers of dynamical processes 
Table 16 

Detectable tracers of turbulent mixing 



Steadfast Hypersensitive 

CO Heavy hydrocarbons (e.g.,CgHj) 

H,Oice CoS 

C 3 S 

CO, 

2 

so 
so 2 

OCN 

Complex organics (e.g., HCOOH) 



In Table [16] we show most promising tracers of transport 
processes in protoplanetary disks as found with our model- 
ing. The quantity least biased by observational and model- 
ing uncertainties is the ratio of the observed column densi- 
ties of an abundant steadfast species to that of a hypersensi- 
tive molecule. Among the insensitive species (Table [5]) the 
most promising are CO and the water ice as these termi- 
nal species incorporate substantial fractions of the elemen- 
tal carbon and oxygen in disks, and are easy to observe. 
Cold molecular hydrogen is not observable, whereas con- 
centrations of C + , light hydrocarbons, CN, HCN, and HNC 
are more sensitive not to mixing, but to the stellar X-ray 
and UV radiation (e.g., iFogel et all [201 It lAresuetal.ll20lU 
iKampet aT]|2011l) . Other molecules and ices unresponsive to 
mixing have low abundances or also model-dependent (e.g., 
S-bearing ices). Both pure rotational CO lines at (sub-)mm 



and ro-vibrational CO lines at IR wavelengt hs have been de- 
tected in inner and outer disk region s (e.g.. iPietu et al.l 120071 : 
lOberg et all 120 lit ISalvk et al .1120111) . Water ice absorption 
featur e at 3 nm has also been detected in protoplanetary disks 
(e.g.. lTerada et alj|2007l) and e nvelopes around young proto- 
stars (e.g., Boogert et al. 2008|). 

Among species sensitive and hypersensitive to the turbu- 
lent transport (Tables HH7} most promising are gaseous and 
solid heavy hydrocarbons, C 2 S, C 3 S, SO, S0 2 , C0 2 , 2 , 
and complex organic molecules (e.g., HCOOH). Basically, 
any non-terminal, abundant molecule produced mostly via 
grain-surface kinetics can be used as a tracer of dynamical 
transport in protoplanetary disks. The OH + and H2O" 1 " ions 
are sensitive to the disk ionization structure and thus can- 
not be reliable tracers of turbulent mixing. Due to atmo- 
spheric opa city, molecular oxy gen can only be observed from 
space (e.g.. lLarsson et al.ll2007l) . Gas-phase mol ecules actu- 
ally d etected in disks are warm C0 2 at IR (e.g.. ISalvk et al] 
12011b and SO at millimeter wavelengths dFuente et alll2010l) . 
Rotational lines of C2S, C3S and various hydrocarbons and 
hydroc arbon anions have been detected in cold dense cores 
(e.g. iDickens et al.l l200ll: iKalenskii etaD 12001: ISakai et al.l 
2010). The emission lines of SO, SO2, and several organic 
molecules have been ident ified in submillimeter spectra of 
young stellar objects (e.g.. |j0rgensen et al] 120051) . CO2 ice 
feature at 1 5 .2 pm has been detected in molecular clouds (e.g., 
iKim et aT]|201 ll) and toward em bedded young low-mass stars 
(e.g.,[Pontoppidan et al. 2008b). In addition, plethora of other 
ices like CH 4 , S0 2 , HCOOH, H 2 C O, and CH 3 OH have been 
observed in Class I/II objects (e.g., Zasowski et al.ll2009l) . 

The column density ratios of these sensitive molecules to 
the steadfast CO and H 2 ice vary between our laminar and 
fast 2D-mixing models by 2-4 orders of magnitude (see CDR 
values in Table [7])- In the not so far distant future forthcom- 
ing observational facilities like ALMA, extended VLA, and 
James Webb Space Telescope will allow us to observe many 
of the mixing tracers along with CO and water ice in nearby 
protoplanetary disks of various ages, masses, and sizes, and 
testify our predictions. 

6. SUMMARY AND CONCLUSIONS 

We study the influence of dynamical processes on the 
chemical evolution of protoplanetary disks. Our analysis is 
based on the 2D flared a-model of a ~ 5 Myr DM Tau disk 
coupled to the large-scale gas-grain chemical code. To ac- 
count for production of complex molecules, the chemical 
network is supplied with a set of surface reactions (up to 
HCOOH, CH3OH, CH3OCH3, etc.) and photoprocessing of 
ices. For the first time our disk model covers a wide range 
of radii, 10-800 AU, and includes warm planet-forming zone 
and cold outer region. Turbulent transport of gases and ices 
is modeled using the mixing-length approximation in full 2D 
(based on the a-prescription for viscosity). Since turbulent 
transport efficiency of molecules in disks is not well known, 
we consider two dynamical models with the Schmidt number 
of 1 and 100. We come up with a simple analysis for laminar 
chemical models that allows to highlight the potential sensi- 
tivity of a molecule to turbulent transport. It is shown that 
the higher the ratio of the characteristic chemical timescale 
to the turbulent transport timescale for a given molecule, the 
higher the probability that its column density will be affected 
by dynamical processes. With our chemo-dynamical models, 
we find that the turbulent transport influences abundances and 
column densities of many gas-phase species and especially 
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ices. The results of the chemical model with reduced turbulent 
diffusion are much closer to those from the laminar model, but 
not completely. Mixing is important in disks since a chemi- 
cal steady-state is not reached for many species due to long 
timescales associated with surface chemical processes and 
slow evaporation of heavy molecules (f > 10 5 years). When 
a grain with an icy mantle is transported from a cold disk 
midplane into a warm intermediate/inner region, the warm- 
up makes heavy radicals mobile on dust surface, enriching 
the mantle with complex ices, which can be released into the 
gas in appropriate temperature regions. In contrast, simple 
radicals and molecular ions, which chemical evolution pro- 
ceed solely in the gas phase, are not much affected by dy- 
namics. We divide all molecules into 3 distinct groups with 
respect to the sensitivity of their column densities to the dif- 
fusive mixing. The molecules that are unresponsive to dy- 
namical transport include such observed and potentially de- 
tectable molecules as C 2 H, C + , CH 4 , CN, CO, HCN, HNC, 
H2CO, OH, as well as water and ammonia ices. Their col- 
umn densities computed with the laminar and fast 2D-mixing 
model do not differ by more than the factor of 2 - 5 ("stead- 
fast" species). The molecules which vertical column densi- 
ties in the laminar and dynamical models differ by no more 
than 2 order of magnitude include, e.g., C2H2, some car- 
bon chains, CS, H 2 CS, H 2 0, HCO\ HCOOH, HNCO, N 2 H + , 
NH3, CO ice, H 2 CO ice, CH3OH ice, and electrons ("sensi- 
tive" species). Molecules which column densities are modi- 
fied by diffusion by more than 2 orders of magnitude include, 
e.g., C 2 S, C 3 S, C 6 H 6 , C0 2 , 2 , SiO, SO, S0 2 , long carbon 
chain ices, CH3CHO ice, HCOOH ice, 2 ice, and OCN 
ice ("hypersensitive" species). The sulfur-bearing molecules, 
along with polyatomic (organic) molecules frozen onto the 
dust grains, are among the most sensitive species to the tur- 
bulent mixing. The chemical evolution of assorted molecules 
in the laminar and turbulent models is thouroughly analyzed 
and compared with previous studies. We find that the ob- 
served column densities in the DM Tau disk are well repro- 
duced by both the laminar and the mixing disk models. The 
observed abundances of reduced and oxidized cometary ices 
are also successfully reproduced by the both models. A com- 
bination of efficient UV photodesorption, grain growth, and 
turbulent mixing leads to non-negligible amount of molecular 
gases in the cold disk midplane. We propose several observ- 
able or potentially detectable tracers of dynamical processes 
in protoplanetary disks, e.g. ratios of the C0 2 , 2 , SO, S0 2 , 
C 2 S, C3S column densities to that of CO and the water ice. 
Some of these tracers have been observed in disks by the cur- 
rent radiointerferometers and infrared facilities (e.g. PdBI, 
SMA, Spitzer, Keck, VLT) and some will be targeted by the 
Herschel telescope. The detection of complex species (e.g., 
dimethyl ether, formic acid, methyl formate, etc.) in proto- 
planetary disks with ALMA and JWST will be a strong indi- 
cation that chemical evolution of these objects is influenced 
by transport processes. 
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